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The  reason  for  dividing  this  report  into  three  volumes  is 
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are  not  required  except  as  a  record  for  those  who  would  care  to 
investigate  the  subject  more  extensively. 
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CHAPTER  VI 

SYNOPTIC  SCALE  VARIABILITY  OF  WINDS  WITH 
TIME  AND  HEIGHT  DIFFERENCES  AS  PARAMETERS 

Chapter  VI  is  devoted  to  extending  the  analysis  of  wind 
variability  from  the  case  in  which  the  two  points  concerned 
are  located  on  an  isobaric  surface  and  winds  are  observed 
simultaneously  to  the  more  general  case  in  which  the  winds  are 
observed  at  two  points  that  are  on  different  isobaric  surfaces 
and  may  be  measured  at  significantly  different  times. 

The  wind  variability  in  such  cases  has  been  studied  exten¬ 
sively  for  the  particular  cases  of  the  variability  in  time  at  a 
particular  point.  Another  important  special  case  that  has  been 
studied  extensively  is  that  of  the  variability  of  wind  with 
height  at  a  given  location  (essentially  simultaneous  measurements 
of  winds  over  a  fixed  point  at  different  levels).  The  references 
available  on  these  two  special  cases  are  extensive  and  are  well 
summarised  by  Major  Hugh  W .  Ellsaesser. ^  The  correlation  of 
winds  obtained  in  each  of  these  two  special  cases  show  peculi¬ 
arities  which  require  some  explanation.  This  explanation  is 
far  from  obvious  when  these  special  cases  are  considered  sepa¬ 
rately,  but  when  the  whole  picture  is  considered  the  major 
variation  of  the  observed  correlations  fit  a  reasonably  straight¬ 
forward  analysis.  In  other  words,  to  explain  the  phenomena 
observed  in  the  variability  of  winds  with  time  and  with  height 
requires  a  complete  four -dimensional  analysis  of  the  variability 
of  wind,  three  space  co-ordinates  and  a  time  co-ordinate. 
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There  is  a  tendency  to  "equate"  time  variability  of  wind 

f  1-2 1 

and  variability  of  wind  with  distance.'  ‘  Though  such  a 
rough  and  ready  relation  Is  handy  for  some  purposes,  it  will  be 
seen  that  the  phenomena  concerned  is  somewhat  too  complex  and 
that  a  more  detailed  separation  of  the  influences  is  required. 
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A.  VERBAL.  DESCRIPTION  OP  STRUCTURE  QF  CORRELATION  TENSOR 

C'QMP'ONENTS^  ™  *  *  “ 

The  basic  item  to  consider  in  treating  the  two-point  vari=» 
ability  of  wind  in  space  and  time  is  the  structure  of  the  com¬ 
ponents  of  the  wind  correlation  tensor.  The  form  of  the  com¬ 
ponents  of  this  tensor  was  discussed  extensively  in  Chapter  III 
and  the  relation  to  the  height  correlation  functions  was  given 
in  equation  (7).  An  ideal  description  of  the  components  of  the 
wind  correlation  tens or  was  presented  in’ Fig.  III-3  for  the 
case  of  the  two  points  on  an  isobarie  surface. 

1 .  Variation,  in  the  Vertical 

To.  Introduce  a  third  dimension  into  this  description 
one  starts  by. depicting  a  family  of  illustrations  like  Fig. 

Ill -3  stacked  on  top  of  each  other.  The  point  P  is  located  at 
the  origin,  of  co-ordinates  (in  the  case  of  each  individual  com¬ 
ponent)  on  the  "center"  sheet  of  tne  stack.  The  point  P'  is 
located  on  some  other  sheet  of  the  stack. 

So  far,  all  sheets  of  the  stack  are  Identical.  Now  modify 
somewhat  all  sheets  (except  the  center  sheet  containing  P)  work¬ 
ing  up  and  down  from  the  center  sheet  in  a  roughly  symmetrical 
manner.  The  first  modification  is  to  "degrade"  the  correlation 
coefficient  values,  i.e.  reduce  the  numerical  value.  This  is  • 
done  progressively  outward  from  the  center  sheet.  In  our 
original  stack,  the  surfaces  of  constant  correlation  coefficients 
were  vertical  cylinders.  When  the  correlation  coefficient  values 
are  "degraded",  these  cylinders  become  closed  over  at  "top"  and 
"bottom"  to  become  somewhat  like  ellipsoids. 


The  second  modification  that  should  be  made  is  t©  tilt  the 
stack  or  to  introduce  a  shear  in  the  sheets.  This  corresponds 
.to  the  well-known  tilt  of  axis  of  highs  and  lows  in  the  atmos¬ 
phere.'  .'The.  direction  of  shear  is  predominantly  in  the  east- 
west  direction-',  (though  not  necessarily  exactly  so). 

At  this  point,  one  may  begin  to  see  one  of  the  character¬ 
istics  of  the  behavior 'of  wind '  component  correlations  in  the 
vertical.  In  modifying  the  components  of  the  wind  correlation 
tensor,  each  .component  was  modified  in  the  same  manner.  But 
in  Pig.  HI-3  the  correlation  contours  for  rxl  and  for  raa  (for 
the  east  pointing  and  the  north  pointing  components  at  both  P 
and  P'  respectively)  are  rotated  at  right  angles  to  each  other. 
The  introduction  of  some  tilt  or  shear  in  the  east-west  direction 
affects  these  two  components  of  the  correlation  tensor  quite 
differently.  In  the  case  of  rlx  the  shear  Is  In  the  same  di¬ 
rection  as  the  longest  direction  of  the  ovals  of  equal  correla¬ 
tion  value  while  in  the  case  of  r8a  the  direction  of  shear  is 
in  the  direction  of  the  shortest  dimension  of  the  ovals  of 
equal  correlation  value.  As  a  result,  as  one  progresses  upward 
(or  downward)  from  the  central  sheet  the  value  of  raa  will  decay 
slightly  more  rapidly  than  the  value  of  rxi.  (This  is  more 
pronounced  in  the  lower  troposphere  where  the  tilt  of  highs  and 
lows  is  greatest.) 

2 .  Variation  with  Time 

The  addition  of  time  as  a  variable  (to  the  case  of  P 
and  ?'  on  the  same  isobaric  surface)  may  be  depicted  in  an 
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analogous  manner.  First  one  considers  a  transparent  overlay 
for  Fig.  III-3  which  contains  the  four  sets  of  co-ordinate  axes. 
The  point  P  is  the  origin  of  co-ordinate  for  each  of  these  sets 
of  axes,  one  for  each  component  of  the  wind  correlation  tensor. 

The  remaining  steps  are  similar  to  those  described  in  con¬ 
nection  with  the  wind  correlation  in  the  vertical,  degradation 
and  shift.  First,  consider  the  correlation  contour  values  to 
be  reduced  fin  magnitude)  with  the  passage  of  time.  Second, 
consider  the  degraded  contours  to  be  shifted  eastward  (roughly) 
with  the  passage  of  time.  The  correlation  coefficient  of  a  wind 
component  at  P  with  a  wind  component  at  P'  is  indicated  by  the 
degraded  and  shifted  correlation  contour  patterns.  Note  that 
the  geographical  locations  of  P  and  P'  are  fixed  so  that  these 
two  points  are  located  on  the  overlay  that  has  the  co-ordinate 
axes  . 

The  correlation  coefficients  that  would  be  obtained  in 
considering  the  variation  of  wind  at  a  single  point  are  given 
when  P'  lies  at  the  origin  of  co-ordinates  (coincident  with  P). 
In  this  instance,  it  is  readily  seen  that  the  value  of  ria  drops 
off  more  slowly  than  r33  because  the  direction  of  shift  of  the 
pattern  with  the  passage  of  time  is  roughly  in  the  direction  of 
the  elongated  contours  of  rtl  and  transverse  to  the  direction 
of  elongation  of  the  contours  of  r33. 

The  situation  for  observations  made  at  the  points  P'  2k 
hours  ahead  of  that  at  P  (Columbia,  Mo.,  at  the  Star)  is  illus¬ 
trated  in  Fig.  VI-1  in  the  case  of  components  rlx  and  r38 
(Winter,  100  mb). 
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If  one  combines  variation  with  time  and  distance,  then  the 
points  P  and  P'  are  separate.  The  point  P  is  at  the  origin  of 
co-ordinates.  The  correlation  values  for  each  component  are 
read  at  the  location  of  P'  (not  P)  from  the  shifted  and  degraded 
correlation  contours. 

The  situation  becomes  extremely  complex  when  time  varia¬ 
bility  at  a  point  is  considered  and  the  correlations  are  com¬ 
puted  along  and  across  the  wind  direction.  However,  even  in 
such  cases,  for  the  upper  air,  the  wind  direction  is  sufficiently 
consistent  to  show  a  more  rapid  decrease  in  the  cross-wind  com¬ 
ponent  correlations  than  in  that  along  the  wind. 

3 •  Some  Further  Considerations 

In  the  preceding  verbal  description  of  the  construction 
of  the  components  of  the  wind  correlation  tensor  as  functions  of 
differences  in  level  and  time,  only  the  most  elementary  aspects 
of  the  situation  are  described.  It  will  be  noted  that  reference 
is  made  to  Fig.  III-3  and  not  to  Fig.  III-4  where  the  tensor 
components  are  in  terms  of  transverse  and  longitudinal  wind 
components.  When  the  height  and  time  variables  are  introduced, 
the  concept  of  longitudinal  and  transverse  wind  components  be¬ 
come  unnecessarily  complex,  more  complex  than  reverting  back  to 
the  fact  that  the  point  of  maximum  correlation  departs  appreciably 
from  the  location  of  the  point  P.  In  order  to  maintain  the  idea 
of  longitudinal  and  transverse  wind  components,  it  would  be  nec¬ 
essary  to  measure  these  along  and  across  the  direction  joining 
the  point  of  maximum  correlation  to  the  point  P'.  This  direction 


is.  not  known-  bsfone  the  fact,  Purthsr,  oven. whan  iho  point  of  ' 
maximum  correlation  is  known,  the  fact  that  it  is  not  at  P  makes 
the  practical  usefulness  of  the  concept  less  convenient  than  the 
reference  ordinary  geographical  co-ordinates. 

In  the’  case  of  P  and  P.'  on  the  same  isobaric  surface,  the' 
fact  that  the  standard  deviation  of  the  wind  components  is  a 
function  of  geographical  location  introduced  only  some  rather 
minor  terms  into  expressions  for  the  two  components  of  the  wind 
correlation  tensor.  The  variation  of  the  standard  deviation  of 
the  wind  component  in  the-  vertical  is  several  times  larger  than 
in  the  horizontal'.  Consequently,  this  non -homogeneous  structure 
of  the  atmospheric  wind .  field .  .brings  about  some  important  modi¬ 
fications  v;hen  P  and  P'  are  .Separated  In  the  vertical. 
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B.  YALIDITy  OF  PRSVTQUS  KESUHPS  M  QmSTKQSmQ  M  SA3B 

The  basic  equation  relating  the  wind  covariance  tensor  com¬ 
ponents  to  the  height  correlation  function  is  of  the  same  form 
as  that  of  (7)  in  Chapter  III.  The  justification  for  this  same 
form  is  easily  seen  by  considering  its  derivation  in  Appendix 
A -I;  We  repeat  here  only  the  modifications  that  would  have  to 
be  made  in  the  arguments  of  Appendix  A-I  to  show  the  relatively 
trivial  changes  involved. 

Equation  (7)  of  Chapter  III  is  repeated  for  convenience 

here . 


The  significance  of  the  quantities  appearing  here  is  the  same  as 
that  tabulated  on  Page  III-7.  The  geometrical  relations  involved 
are  illustrated  in  Pigs.  Ill -2  and  VI -2 (a). 
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The  basic  equation  from  (8)  of  Appendix  A-I 


ei  3*  ej  a 


3  ^3  (Uk  UlT  ~  S2^  1?  j  VVj 


So  ..  f  SRn  h  5o~ '  5Rh  h 

^  i?r  ’  cte[  wr 


,  3a  3o'ti  , 

+  WT’Wf**  + 


Tso^ .  SRn.i 
[35c7  3xt 


o'£f2jLi_ 

a?33x. 


in  which  el31c  (and  eja£)  -  the  alternating  tensor  and  has  the 
value  1  If  1  =  1,  k  =  2;  -  1  if  i  =  2,  k  =  1;  0  otherwise. 
uk  =  wind  components  at  P,  k  =  1,  2;  =  wind  components  at 

P',  The  symbol  (ukUj)  is  the  mean  value  of  the  product,  and 
the  factors  of  the  product  are  departures  from  their  respective 
means  (these  are  the  components  of  the  wind  covariance  tensor), 
g  =  acceleration  of  gravity;  a,  a '  are  standard  deviations  of 

« 

height  at  P  and  P'  respectively;  Rhh  =  height  correlation  func¬ 
tion,  the  correlation  coefficient  relating  height  of  an  isobaric 
surface  at  P  to  that  at  P';  xt  =  co-ordinates  at  P  and  Xj  =  co¬ 
ordinates  at  P' . 

The  basic  idea  at  this  point  was  (in  Appendix  A-I),  that 
Rhh  is  a  function  of  the  relation  of  P'  to  P  and  only  secondarily 
a  function  of  the  location  of  P  itself. 

Symbolically 


-  ^-h  h  (xi  »xi  >x3  )  —  Uti  h  (  ?i  >  ?s  *  xl  >  xs  )  ~  ^iih(?i>?a) 


where 


?  _ 
n 


“  xx  -  x, 


-  Xa  -  X3  . 

Some  assumptions  about  the  particular  relation  Rhh  (|1;§ a)  were 
then  made.  It  was  found  that  if  r2  =  ?!  2  +  ?22  =  square  of  dis¬ 
tance  from  P  to  P'  and  Rbb  =  Rbb(r),  the  expressions  for  (uku^) 
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were  particularly  simple.  (See  Fig.  HI-3  'and  equation  (13) 
of  Chapter  III.) 

Four*  new  variables  have  been  Introduced  when  height  and 
time  separation  are  involved.  These  are  z,  z' ,  t,  t';  height 
and  time  for  P  and  P' .  We  will  not  Introduce  these  specifi¬ 
cally  in  the  symbolism  since  differentiation  with  respect  to 
these  are  not  involved  in  (l).  The  presence  of  these  variables 
implies  that  the  largest  correlation  between  height  at  P  and 
height  at  P/  is  no  longer  reached  if  P'  s  p.  (Formerly 
Phh(0)0)  =  !•)  Let  the  largest  value  of  Rhn  be  attained  at 
some  point  Q  (see  Fig.  VI-2(b)).  Let  the  vector  PQ  be  given 
by  a  =  (aj,aa).  And  let  the  vector  QP'  be  given  by  p  =  (pa,pa). 
Then  (Fig.  VI -2 ) 

=  xt  +  a !  +  px , 

x'  =  x  +  a  +  p,  or 

-  xa  +  a8  +  pa  . 

It  is  now  convenient  to  describe  the  correlation  contours  of 
Rhh  =  constant  in  terms  of  a  co-ordinate  pair  (p,d)  where 
pa  =  pi2  =  pi 2  +  pss  and  a  =  p,Aj 

=  p2K2) 

where  vt  is  a  unit  vector  giving  the  azimuth  of  P'  from  Q; 
is  a  preferred  reference  direction  (an  axis  of  symmetry  of 
Rh  h  =  constant,  for  example).'  (As  =  ( \  ,  7\a  )  has  nothing  to  do 
with  A3 ,  X3 ,  values  of  the  Coriolis  parameter;  a  here  is  tem¬ 
porarily  used  as  a  co-ordinate,  projection  of  QP'  on  the  ref¬ 
erence  direction,  and  has  nothing  to  do  with  standard  deviation 
of  height.) 
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h,  =  ^SiUu.  So  +  . 

3p  3d  S'ai 


Since 


p2  —  (?i-ai)2  (?a”aa)2  an^  0  —  (?i“ai)^i  +  ( ?a  “aa  )^a 


then 


S£_  _  fii.  =  v 

S§i  p  Vl# 


3CT  _  -v 

afT  ■  x‘ 


where  vt  is  the  vector  giving  the  azimuth  of  P'  from  Q.  Then 

SRjlX  -  SR),  t  ,  5Ph  h  -> 

a?t  ap  1  aa  A» 

which  is  exactly  the  same  as  (9)  of  Appendix  A-I  except  that 
(p,a)  replaced  (r,s).  The  remainder  of  the  transformation  is 
only  tedious  and  involves  no  new  ideas  not  already  covered. 

The  end  result  is  that  to  estimate  the  wind  covariance  compon¬ 
ents,  equation  (l)  is  valid  provided  that  the  center  of  the  co¬ 
ordinates  (r,s)  is  at  Q,  a  "floating  point"  not  specified  exactly 
a  priori,  and  is  that  point  where  Rhh  has  its  maximum  value. 

To  rotate  co-ordinates  along  which  wind  components  are  re¬ 
solved  into  longitudinal  and  transverse  form  it  is  required  to 
multiply  (l)  by  the  proper  unit  vectors.  In  (l)  the  unit  vector 
Vj  gives  the  azimuth  of  P'  with  respect  to  Q,  the  center  of 
maximum  Rhh.  If  one  multiplies  (i)  by  VjVj  and  sums  over  the 
subscripts  i  and  j  (since  they  are  now  repeated),  the  same  re¬ 
sults  are  obtained  as  before  as  far  as  their  form  is  concerned. 
This  gives  the  covariance  form  (u.uj).  The  geometrical  inter¬ 


pretation  is  a  bit  awkward,  however.  As  shown  in  Fig.  VI-3, 
the  direction  of  the  longitudinal  component  is  QP'  (and  the 


direction  of  the  transverse  component  would  be  perpendicular 
thereto).  This  is  by  no  means  the  direction  FP'. "The  longi¬ 
tudinal  and  transverse  component  form  was  convenient  because, 
among  other  things,  these  could  be  thought  of  as  head -wind  and 
cross-wind  components  along  a  flight  path.  This  would  be  true 
now  only  if  P,  Q,  P'  were  colinear. 

If  winds  are  resolved  along  and  across  PP',  other  diffi¬ 
culties  arise.  Let  be  the  unit  vector  from  P  to  P' .  Mul¬ 
tiplying  (l)  by  \>l  one  obtains  (u^u')  on  the  left-hand  side, 
but  the  direction  factor  in  the  product  on  the  right-hand  side 
in  the  case  of  the  first  term  becomes 

v1vJ(v1vJ)  =  (vjOjHvjOj)  =  cos2  (9  -9* ) 
where  9  is  the  azimuth  of  P'  with  respect  to  Q  and  where  0*  is 
the  azimuth  of  P'  with  respect  to  P  (see  Fig.  VI-3).  In  a  sim¬ 
ilar  manner,  the  term  In  61 j  will  have  the  factor  6tjv,vj  =  coss9 
Thus  on  the  whole  the  terms  on  the  right  of  (l)  become  strongly 
direction  dependent  even  when  Rhh  =  ( p )  is  a  function  of 

radial  distance  only.  As  a  consequence,  the  simplifications 
that  were  present  when  P  and  P'  were  on  the  same  isobaric  sur¬ 
face  are  completely  lost. 


Pig.  VI-3  Geometric  Relations  Between  the  Locations 

of  P,  Q,  and  P1  when  Winds  are  not  Observed 
Simultaneously  or  on  the  Same  Isobaric 
Surface,  or  both 
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C.  HEIGHT  CORRELATION  EmMJMJjiiiS  IN-  IffgBIS.  MD 

Tusur  .  '  "  —  - 


In  order  to  use  (l)  to  estimate  the  form  of  the  components 
of  the  wind  correlation  tensor  (covariance)  components  or  the 
corresponding  wind  component  correlations,  it  is  necessary  to 
estimate  the  form  of  the  height  correlation  function  Rhh.  This 
is  done  in  only  a  very  simplified  case. 

When  P  and  P'  lie  on  the  same  isobaric  surface,  a  simple 
reasonable  estimate  of  Rhh  appeared  to  be 

R„  =  ke^1'2/'2 

where  0  <  k  s  1  and  L  =  a-1  =  constant  =  size  parameter.  This 
particularly  simple  form  is  now  modified  into  a  form  like 


R 


h  h 


=  k, 


i1  -  ^Xl  - 


zf_\  .  2  r  e/2  p  -z  2/2  D£ 

C3/ 


(3) 


aar2  =  (x^.-x^t-b^)5  +  (*a/-xa-aa_t-aaz)g 

[L(t,z)]a 

The  many  factors  in  this  function  form  serve  only  the  purpose 
of  introducing  effects  that  are  likely  to  be  encountered.  It 
is  perhaps  too  much  to  ask  that  the  details  of  such  an  expres¬ 
sion  be  confirmed. 

The  factor  k0,  0  <  k0  s  1  is  introduced  to  take  into  account 
errors  of  measurement  and  effects  of  small  scale  motions  that  are 
too  small  to  be  adequately  accounted  for  in  the  network  of  data 
points  that  are  available. 

The  factor  (l  -  ta/T2)  is  indicative  of  the  degradation  of 
correlation  at  its  maximum  point  on  an  isobaric  surface  with  the 
passage  of  time. 
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The  factor,  e  '  Is  indicative  of  the  degradation  of 
the  maximum  correlation  value  when  P  and  P'  are  separated  by  a 
height  difference  z. 

The  factor  (l  -  za/£3)  is  to  account  for  the  fact  that  the 
correlations  are  not  degraded  symmetrically  above  and  below  the 
level  z  =  0  (level  of  P) .  This  is  due  to  the  fact  that  the 
standard  deviation  of  height  increases  to  a  maximum  and  then 
decreases.  The  rate  of  decrease  and  increase  are  large  enough 
to  materially  affect  the  shape  of  the  correlation  function. 

(See  Appendix  A-III  where  the  subject  is  discussed  in  detail.) 

The  terms  in  a^  +  btz  allow  for  the  drift  of  the  center  of 
correlation  maximum  as  a  function  of  time  and  also  the  tilt  as  a 
function  of  the  difference  between  the  levels  of  P  and  P' . 

The  factor  L(t, z)  allows  for  spreading  of  the  correlation 
pattern  as  a  function  of  the  time  lags  and  height  lags  (differ¬ 
ences). 

None  of  the  differentiations  in  (l)  Is  with  respect  to  z 
or  t  so  that  all  of  the  factors  containing  these  quantities  may 
be  treated  as  constants .  This  means  that  one  may  actually  ignore 
the  conjectures  made  in  the  specific  formulation  of  (3)  and  write 
it  In  the  form 


Rft  n  =  Ke^V2 

a2r2  =  [  (xx7  -xx  -Aj  ) 2  +  (x'  -xs  -A s  )  2]/L2 


(3a) 


in  which  K,  L,  Ax ,  Aa  are  all  functions  of  t  and  z.  In  particular, 
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one  would  expect  that 


~  <1  -  SX1  -  SKzV2D' 


Aj  ~  aj_  t  +  "bj  z 


A2  ~  aat  +  baz 


L  ~  L(t,z) 

The  form  of  (l)  Is  much  too  complex  for  direct  substitution 
of  (3a)  to  yield  reasonably  tractable  results.  Only  the  terms 
of  primary  importance  are  retained.  These,  as  was  pointed  out 
in  Chapter  III,  consist  of  the  terms  underlined  thrice  in  (l). 

We  also  retain  the  terms  underlined  twice  in  (l)  but  with 
the  assumption  that  the  gradient  of  the  standard  deviation  of 
height  is  small  and  uniform,  i.e.  that  a  as  o',  |Vc|  =  |Vcr'|,  p.,  *  m' . 
Then  substituting  the  values  of  1  and  2  for  the  subscripts  one 


obtains 

tf-K < )  =  - 


+  d^|,AGlngej  +  |7c|2Rhhsin2i|;,  (4a) 


ua 

g 


(uaua)  =  -  aa'  in 2 8  +  -  R|acosg9  +  |Va  |  gR„  h  cos2iji ,  (4b) 

x  ax 


^-(UiUg)  -  -  ua'  sln9cos0  +  l^rl 2rh *sin*cost 

0  mm 


-5— (u.  Uj  )  =  -  aa1  5  R|n  -  X.  a  n  sinecose  +  |^|2Rh.sin(icosi|[ 

r  4  04  U 
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where  0  is  the  azimuth  of  P'  with  respect  to  Q  and  \ji  is  the 
azimuth  of  the  gradient  of  the  standard  deviation  of  height, 

Va. 

Note  that  when  (4a)  and  (4b)  are  added  and  subtracted,  one 
ob  ta  ins 


XX' 


g  8ao' 


>,%')  +  KV)]  - 


and 


(5a) 

(5b) 


AX' 

g2acr 


r  . 

-  ( n  nM 

1  _  .  / 1. 

&Rh  h 

.  3% 

\ui  ui ) 

\  / 

J  \  r 

dr 

dr2  ) 

|i-)cos2e  +  |~|2Rhhcos2\|i 


It  is  worthy  of  note  that  (5a)  is  very  closely  related  to 

(21 

the  stretch  vector  correlation  coefficient v  '  and  is  seen  to  be 
independent  of  the  azimuth  of  P  with  respect  to  Q  when  Rbh  is 
also  independent  of  9.  On  the  other  hand,  even  when  Rhb  is  in¬ 
dependent  of  9,  the  factor  eoc29  makes  (5b)  highly  dependent  on 
the  azimuth  of  P'  with  respect  to  Q. 

The  sum  and  difference  of  (4c)  and  (4d)  likewise  have  dis¬ 
tinct  characteristics.  Thus, 


AX' 


g2acr 


7^ ( Uj_ Uj  )  +  (ug u^  ) ]  -  -(—V1-  "  F’^t^)sin20  "  l^r|aRtmSin2i!r  (5c) 


and 


XX ' 
g2acr 


{(ui«s)  -  (ua )]  =  2  |~|-^k;in(e-ili)  . 


(5a) 


In  (5c)  it  is  readily  seen  that  the  first  term  on  the  right  is 
highly  direction  dependent  while  the  second  term  on  the  right 
is  independent  of  direction  and  is  small  if  |7c/e|  is  small. 

The  expression  (5d)  is  dependent  on  |Vo/cr|  and  would  become  zero 
if  the  gradient  of  the  standard  deviation  of  height  were  com¬ 
pletely  negligible. 
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If  we  use  (3a)  for  Rih  and  convert  to  correlation  coeffi¬ 
cients  rather  than  to  using  covariances,  one  obtains  (on  re¬ 
taining  only  the  most  important  terms*) 


s*  = 

+  raa  =  2K(l-fa/2)e"^/2, 

(6a) 

D*  -  rlx 

-  raa  =  K(fla-fs2)e-f'S/2, 

(6b) 

S  =  r18 

+  r8l  -  2K(f1fa+A)e"?2/2, 

(6c) 

D  =  rl  3 

-  r81  =  2K(B?a+Cf1  )e-^2//2, 

(6d) 

where 

=  (xi'  -Xj  -Aj  )/L,  ?a  =  (x8  -x8  -A 8  )/L,  f2  =  2  +  f82. 

Thus,  the  four  relations  (6a),  ---,  (6d)  contain  seven  un¬ 
determined  parameters.  Pour  of  these  arise  from  the  assumptions 
regarding  the  form  of  Rhh:  K,  L,  Ax ,  Aa;  and  three  depend  on 
the  fact  that  the  horizontal  gradient  of  height  is  not  entirely 
negligible:  A,  B,  C.  Thus, 

2KA  a  L2  |~jSsin2t|i  (7a) 

2KB  s  -  L|~|cos^  (7b) 

2KC  =  L|2£|sin*  (7c) 

One  should  note  that  in  assigning  meaning  to  A,  B,  and  C  in  the 
above  one  has  only  retained  here  the  effect  of  the  nonvanishing 


*The  arithmetic  of  the  reduction  is  tedious.  In  carrving  out 

the  reduction  it  is  necessary  to  go  to  the  forms  (4a),  - , 

(4d)  and  to  account  for  the  component  standard  deviations  by 
using  the  limit  relations  for  P'-P,  Q-Patt=z-0. 
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of  the  gradient  of  the  standard  deviation  of  height*  This 
serves  merely  to  justify  the  fact  that  such  terms  should  he 
retained  in  S  and  D.  Actually,  parameters  A,  B,  C  pick  up 
contributions  due  to  the  fact  that  Rhh  may  he  of  an  elliptical 
shape,  that  the  standard  deviation  of  height  may  he  non-uniform, 
etc.  Consequently,  one  is  inclined  to  look  on  A,  B,  and  0  as 
three  empirical  factors  that  should  be  retained  and  to  de- 
emphasize  any  meaning  that  might  be  ascribed  to  them  in  the 
initial  analysis. 

The  individual  correlation  coefficient  has  the  forms 
r:i  a:  K(1  -  fa2)e"?  /2 

raa  -  K(  1  -  V)e“f2/2 

ria  =s  +  A  +  Brx  +  '2 

r*a i  2=  K[Va  +  A  ~ 

The  first  pair  is  by  far  the  most  important.  Both  ria  and 
rai  are  small  and  their  expressions  are  somewhat  complex  for 


practical  purposes. 
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d.  ■  Baaggggas. . fig,  the: m 

The  values  of  K,  L,  kl ,  and  Aa  were  estimated  from  (6a) 

1 

and  (6b)  by  minimizing  the  sum  of  the  squares  of  the  differences  i 

rf 

4 

between  S*  and  D*  as  computed  and  as  determined  from  the  corre-  j 

lations  actually  observed.  The  method  used  was  that  of  Gauss  for 
determining  parameters  that  appear  in  a  non-linear  form  (see  ( 

Appendix  A-IV  for  a  similar  example). 

The  values  of  A,  B,  and  C  appear  linearly  in  (6b)  and  (6d). 

They  were  estimated  by  standard  least  squares  methods  from  these 
expressions  after  the  values  of  K,  L,  Ax ,  and  Aa  from  (6a)  and 
(6b)  had  been  obtained. 

The  correlation  data  used  throughout  were  those  computed  by 
the  Sandia  Corporation  using  serially  complete  wind  data  prepared 
by  the  U.S.  Weather  Bureau  for  the  Federal  Civil  Defense  Agency 
covering  a  period  of  five  years.'-" 
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E.  TEE  EMPIRICAL  PAHAffBTBBS 

Pigs.  VI-4  through  VI -45  Illustrate  the  results  of  the 
numerical  estimates  of  the  parameters  K,  L,  Ax ,  Aa ,  A,  B,  and 
C.  These  illustrations  are  divided  into  seven  groups,  one  for 
each  parameter.  These  groups  each  contain  6  figures  and  are 
divided  into  three  illustrations  of  summer  conditions  and  three 
for  winter  conditions.  There  are  three  "master  levels"  or 
levels  for  the  point  P,  500  mb,  300  mb,  and  100  mb. 

The  500  mb  level  illustrates  conditions  in  the  lower  trop- 
opause  where  the  gradient  of  standard  deviation  is  increasing 
rapidly  with  height,  the  300  mb  level  represents  conditions  in 
the  upper  troposphere  where  the  standard  deviations  are  near 
their  maximum,  the  100  mb  level  represents  conditions  in  the 
lower  stratosphere  where  the  standard  deviations  are  decreasing 
with  heights . 

Each  illustration  is  divided  into  three  parts,  one  repre¬ 
senting  conditions  with  Washington,  D.C.  as  the  point  P  or 
"master  station",  the  next  with  Omaha,  Nebraska,  and  the  third 
with  Oakland,  California,  serving  in  this  role.  The  identity  of 
the  master  station  has  been  preserved  not  only  to  show  the  dis¬ 
tinction  between  conditions  centered  at  these  locations,  but 
also  to  emphasize  the  great  similarities  that  occur. 

The  characteristics  of  each  parameter  as  functions  of  z, 
height  difference,  and  t,  time  difference,  are  discussed  in  turn 
with  the  exception  of  A1  and  A3  (which  represent  a  displacement 
vector)  and  of  B  and  C  (which  also  represent  vector  components). 
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In  all  cases,  the  functions  of  z  and  t  have  been  determined 
at  25  points  each,  5  time  lags  of  -2,  -1,  0,  Tl,  and  +2  days 
and  5  height  lags  (which  amounts  to  roughly  -4,  -2,  0,  +2,  +4 
kilometers,  but  which  differ  from  one  master  level  to  another 
since  winds  on  standard  isobaric  surfaces  are  used.  .* 

1 .  Characteristics  of  Parameter  k 

The  variation  of  the  parameter  k  with  height  and  time 
lag  is  illustrated  in  Fig.  VI-4  through  VI-9.  There  are  common 
characteristics  present  in  all  of  these  figures.  The  most  im¬ 
portant  is  the  decay  of  the  value  of  k  with  both  height  and  time 
lags.  The  largest  value  is,  of  course,  that  for  zero  height  and 
time  lag  located  at  the  "center"  of  the  5x5  array  of  points. 

(a )  100  mb  Level 

The  most  important  property  that  is  typical  of 
the  100  mb  level  is  the  somewhat  more  rapid  decay  of  k  for  P' 
located  above  P  as  compared  with  the  location  of  P'  below  P. 

This  characteristic  is  prominent  at  all  locations  and  seasons, 
though  somewhat  more  notieable  in  summer  (Fig.  VI-4)  than  in 
winter  (Fig.  VI -5).  This  is  one  of  the  effects  brought  about 
by  the  non-homogeneity  of  the  wind  field  In  the  vertical;  speci¬ 
fically,  the  fact  that  at  the  100  mb  level  In  both  seasons  the 


standard  deviation  of  wind  is  decreasing  with  height 


*The  "slave  levels"  for  each  master  level  are: 


Master  Level  (P) 


Slave  Level  (P' 


500  mb 
300  mb 
100  mb 


850,  700,  500,  400,  300  mb 
500,  400,  300,  250,  200  mb 
200,  180,  100,  80,  50  mb 
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(b)  3QQ  mb  Level 

The  300  mb  level,  located  reasonably  close  to 
the  level  of  maximum  wind  and  of  the  maximum  standard  deviation 
of  wind,  displays  the  characteristic  that  the  decrease  of  k 
with  increasing  numerical  value  of  time  lag  is  nearly  the  same 
for  both  positive  and  negative  values  (P'  above  or  below  P) . 

The  values  of  k  for  both  summer  (Fig.  VI -6)  and  winter  (Fig. 
VI-7)  appear  to  behave  in  much  the  same  way  so  that  seasonal 
changes  appear  to  be  minor. 

(c)  500  mb  Level 

At  the  300  mb  level  in  both  summer  (Fig.  VI-8) 
and  winter  (Fig.  VI-9)  the  change  of  k  with  height  lag  is  just 
the  reverse  of  that  illustrated  at  the  100  mb  level.  At  the 
500  mb  level  the  value  of  k  drops  off  more  rapidly  for  negative 
height  lags  (P'  below  P)  than  for  positive  height  lags.  This 
is,  of  course,  due  to  the  fact  that  at  300  mb  the  standard  de¬ 
viation  of  wind  is  increasing  with  height  In  both  seasons. 

2 ,  Characteristics  of  Parameter  L 

The  major  over-all  characteristic  of  the  parameter  L, 
the  size  parameter  in  the  horizontal,  that  appears  to  some  ex¬ 
tent  in  Figs.  VI-10  through  VI-13  is  the  apparent  uniformity. 
There  appears  to  be  a  trend  for  L  to  be  larger  as  the  magnitude 
of  the  time  difference  between  observations  increases,  but  this 
is  not  invariably  observed.  The  behavior  of  L  in  detail  is 
characterized  primarily  by  the  level  of  the  point  P.  The  three 
particular  cases  are  discussed  separately. 
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One  characteristic  that  appears  at  each  level  to  a  greater 
or  less  degree,  and  which  warrants  special  mention  is  the  tend¬ 
ency  of  the  size  parameter  to  increase  somewhat  more  rapidly 
with  negative  time  lag  at  Oakland  than  it  does  when  time  lags 
are  positive.  In  contrast  with  this  at  Washington  the  value 
of  L  tends  to  increase  more  rapidly  with  positive  time  lags 
than  with  negative.  This  unsymmetrical  change  of  L  with  time 
lag  at  these  locations  may  be  a  phenomena  introduced  by  the 
geometrical  relation  of  the  family  of  points  P'  v/ith  respect 
to  the  point  P.  The  points  P'  are  in  each  instance  radiosonde 
(rawinsonde)  locations  in  North  America.  As  discussed  in  the 
next  section,  the  location  of  the  correlation  "center"  (the 
point  of  maximum  correlation  for  rlx  and  r8a)  lies  to  the  we st 
of  the  point  P  for  negative  time  lags  and  to  the  east  of  P  for 
positive  lags.  The  distance  Is  roughly  proportional  to  the  lag 
magnitude.  This  means  that,  with  negative  lags  at  Oakland  and 
positive  lags  at  Washington,  the  "center  of  correlation"  lies 
to  sea  from  the  data  used  to  estimate  the  parameters.  This 
indicates  that  instead  of  estimating  L  from  data  points  surround¬ 
ing  the  "center  of  correlation",  all  of  the  available  data  lie 
to  one  side  of  the  center.  It  may  well  be  that  this  situation 
has  had  a  tendency  to  produce  values  of  L  which  are  larger  than 
would  be  really  present  if  data  had  been  available  surrounding 
the  center  of  correlation  maximum  symmetrically.  (The  same  un¬ 
certainty  Is  present  in  the  calculation  of  values  for  Al  and  Aa , 
the  co-ordinates  of  the  center  of  correlation  maximum  for  rn 
and  ra  8  . ) 
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(a)  100  mb  level  -  / 

The  variation  of  the  parameter  L  with  height  and 
time  lags  is  illustrated  in  Figs.  VI-10  and  VI-11  for  summer 
and  winter  respectively.  In  both  summer  and  winter  the  size 
parameter  (L)  tends  to  increase  with  increasing  height  lag, 
the  rapid  rate  increase  taking  place  at  the  positive  height 
lags.  The  situation  is  much  more  pronounced  in  summer  (Fig. 
VI-10)  than  in  winter  (Fig.  VI-ll).  It  is  interesting  to  note 
that  such  a  shift  of  the  size  parameter  implies  that  the  funda¬ 
mental  dynamics  of  the  air  motions  are  undergoing  a  fundamental 
change  (or  are  about  to  undergo  such  a  change)  at  the  levels 
where  such  a  change  in  size  parameter  appears. 

The  strong  tendency  for  unusually  large  L  with  large 
negative  time  lags  at  Oakland  and  large  positive  time  lags  at 
Washington  shows  up  strongly  in  Fig.  VI-11.  It  is  not  without 
Interest  to  note  that  the  situation  at  Omaha  Is  at  the  same 
time  unsymmetrical  and  corresponds  to  the  situation  at  Washing¬ 
ton.  This  lack  of  symmetry  at  Omaha  cannot  be  explained  by  the 
lack  of  data  eastward  of  Omaha  because  this  is  the  area  In 
which  the  data  density  is  greatest.  One  might  conjecture  that 
the  effect  of  the  mountain  region  lying  between  Oakland  on  the 
west  and  Omaha  and  Washington  on  the  east  may  contribute  to 
this  effect,  but  the  detailed  mechanism  to  account  for  this  is 
by  no  means  clear. 

(b )  300  mb  Level 

The  variation  of  L  with  time  and  height  lags  for 
summer  (Fig.  VI-12)  and  winter  (Fig.  VI-13)  at  the  300  mb  level 
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are  the  least  interesting  of  the  three  levels.  The  value  oT  X 
increases  slowly  with  the  magnitude  of  the  time  lags  and  shows 
little  variation  with  height  lag. 

(c)  500  mb  Level 

Both  the  summer  and  winter  situations  of  Pigs. 
VI-14  and  VI-15  respectively  indicate  a  tendency  for  larger 
values  of  L  at  the  negative  height  lags  (near  the  surface). 

This  is  just  the  opposite  to  the  situation  encountered  at  the 
100  mb  level.  This,  of  course,  corresponds  to  the  difference 
in  structure  of  atmospheric  disturbances  in  the  surface  layers 
and  that  of  the  middle  troposphere  where  surface  effects  have 
lost  a  great  deal  of  this  effectiveness. 

It  should  be  pointed  out  that  the  variation  of  correlations 
rlt  and  raa  in  the  vertical  include  an  "interaction"  between  the 
variation  of  both  the  parameters  k  and  L  as  functions  of  height. 

3 .  Characteristics  of  Parameters  A,  and  AP 

The  parameters  A1  and  Aa  are  discussed  together  since 
they  are  the  components  of  the  displacement  for  the  center  of 
maximum  correlation  rlx  and  rsa  from  the  point  P.  The  major 
feature  of  the  displacement  is  the  movement  from  west  to  east 
through  P  with  increasing  time  lag.  There  is  a  pronounced 
tendency  for  displacement  from  east  to  west  v/ith  increasing 
height  lag.  This  is  particularly  prominent  at  the  500  mb  level. 

The  situation  of  the  "center  of  correlation"  and  the  data 
area,  mentioned  in  connection  with  the  parameter  L,  should  be 
kept  in  mind.  Generally  for  negative  lags  the  center  of  corre¬ 
lation  at  Oakland  is  well  west  of  the  area  of  data  while  at 
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Washington  with  positive  lags  the  center  is  well  east  of  the 
area  of  data.  In  these  cases  it  is  felt  that  the  displacement 
values  Ax,  Aa  (and  all  other  parameters,  for  that  matter)  may 
he  unreliable. 

In  the  case  of  parameter  Aa  ,  many  cases  occur  where  the 
values  are  quite  small.  In  such  cases  the  reality  of  any  trend 
with  height  or  time  lag  is  seriously  questioned.  Such  trends 
may  be  entirely  due  to  the  errors  in  the  basic  data  and  may  not 
be  real.  An  estimate  of  the  error  in  the  values  of  Ax  and  Aa 
may  be  obtained  from  the  values  listed  in  the  figures  at  lag 
zero  in  time  and  height.  For  zero  height  and  time  lags  the 
center  of  maximum  correlation  in  rxx  and  raa  should  be  at  P. 
i.e.  Ax  -  Aa  =0.0.  The  computing  program,  however,  did  not 
supply  this  information.  Consequently,  the  computing  procedure 
computed  what  Ax ,  Aa  should  be  even  in  this  case  from  the  cor¬ 
relation  values.  The  root  mean  square  values  of  Ax  and  Aa  to¬ 
gether  Is  34  n  mi  and  is  considered  as  measure  of  the  error 
included  In  these  parameters. 

(a)  100  mb  Level 

The  summer  displacements  at  100  mb,  Figs.  VI-16 
and  VI-17,  tend  to  be  somewhat  anomalous.  The  trend  for  east¬ 
ward  displacement  (Ax )  of  about  100  n  mi/day  appears  generally, 
but  there  are  pronounced  exceptions.  For  positive  height  lags 
at  Oakland  the  displacements  are  of  the  wrong  sign  for  positive 
time  lags.  These  particular  values  are  also  associated  with 
unusually  large  values  of  L  (see  Section  2a). 
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The  values  of  As  In  this  instance  are  also  anomalous > 
tending  to  he  negative  and  large  too.  The  values  of  A2 ,  ex¬ 
cept  at  Washington,  tend  to  be  increasingly  negative  for  both 
positive  and  negative  time  lags.  This  would  lead  to  a  curved 
path  for  the  center  of  correlation.  Only  at  Washington  is  the 
displacement  of  the  center  of  correlation  definitely  from  the 
northwest  to  the  southeast.  (This  case  is  exceptional  at 
height  lag  of  +2 . ) 

The  winter  picture  in  Pigs.  VI-18  and  VI-19  is  much  more 
definite.  The  value  of  A1  shows  rapid  eastward  displacement 
of  the  center  of  correlation.  The  rate  appears  to  be  more 
rapid  at  Washington  and  least  at  Omaha  (near  350  and  170  n  mi/day 
respectively  with  a  value  near  220  n  mi/day  at  Oakland).  There 
also  appears  to  be  a  tendency  for  more  rapid  movement  in  the 
interval  from  lag  -1  day  to  lag  +1  day  than  in  the  outlying  lag 
intervals . 

All  cf  the  values  of  As  (Fig.  VI-19)  are  positive  (as  con¬ 
trasted  to  negative  values  in  the  summer)  with  a  strong  tendency 
for  the  larger  values  at  the  larger  numerical  time  lags.  This 
'would  indicate  that  as  the  center  of  correlation  moves  from 
west  to  east  it  approaches  from  the  north  of  west  (negative 
time  lags)  and  departs  (positive  time  lags)  toward  the  north  of 
east . 

(b )  300  mb  Level 

The  trend  of  the  center  of  correlation  at  the  300 
mb  level  is  much  more  definite  in  both  summer  and  winter.  Trie 
summer  situation,  Pigs.  VI -20  and.  VI -21,  shows  a  moderate  west 
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to  east  motion  of  near  150  n  ml/d  ay  (again  more  rapid  near  time 
lag  zero  than  for  larger  numerical  values).  The  vertical  dis¬ 


placement  at  Oakland  appears  to  be  slightly  toward  the  east 
while  at  Washington  and  Omaha  it  appears  to  be  tilted  toward  the 
west  (expected).  The  values  of  Aa  (Fig.  VI-21)  indicate  a  trend 
from  north  to  south  with  time  lag  and  from  south  to  north  with 
height  lag.  This  trend  at  Oakland  is  particularly  small.  As  a 
consequence,  the  center  of  correlation  moves  from  somewhat  north 
of  west  to  somewhat  south  of  east.  The  systems  are  tilted  from 
southeast  to  northwest. 

In  the  v.'inter  (Pigs.  VT-22  and  VI-23)  the  same  situation 
prevails  except  with  a  greatly  increased  speed  of  movement  in 
the  eastward  direction.  The  decreased  rate  of  movement  from  -2 
to  -1  and  +1  to  +2  is  particularly  pronounced  (Pig.  VI-22).  The 
north-south  displacement  with  time  lag  is  also  smaller  than  in 
summer  except  in  the  case  of  Washington  when  the  sign  of  the 
time  displacement  is  reversed,  being  from  south  to  north. 

( c )  300  mb  Level 

The  situation  at  500  mb  Is  similar  to  that  at  300 
mb  in  its  general  aspect.  The  summer  situation  (Figs.  VI -24  and 
VI-25)  shows  the  eastward  movement  of  the  center  of  correlation 
with  a  tendency  to  move  also  from  somewhat  north  of  west  to 
definitely  south  of  east.  The  tilt  is  reasonably  strong  from 
southeast  to  northwest.  (The  value  of  kl  at  height  lag  -1  and 
time  lag  +2  at  Omaha  is  definitely  anomalous.  The  same  appears 
to  be  true  of  the  values  of  At  at  Oakland  for  height  lag  of  -1, 
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time  lag  +2  and  height  lag  -2  and  time  lags  -1  and  -2.  These 
values  were  carefully  checked,  but  the  source  of  the  error  was 
not  located.)  In  the  case  of  Oakland,  the  north-south  compon¬ 
ent  of  motion  is  much  smaller  than  at  the  other  locations . 

The  winter  situation.  Pigs.  ¥1-26  and  VI-27,  appears  to 
be  analogous  to  that  at  300  mb  in  general.  At  Oakland  the  move¬ 
ment  of  the  center  of  correlation  is  from  west  to  east  at  about 
200  n  ml/day  with  a  tilt  from  southeast  to  northwest.  The  same 
situation  prevails  at  Omaha  with  a  much  more  pronounced  tilt. 

At  Washington  the  motion  is  much  more  rapid,  near  400  n  mi /day, 
and  with  a  trend  toward  a  curved  path  from  near  west-east  curving 
to  the  northward  for  positive  time  lags.  At  Washington  the  tilt 
tends  to  be  from  south  to  north  with  little  west -east  component 
except  for  the  largest  negative  height  lags. 

In  Table  I  the  values  of  the  parameters  a  and  b  are  broken 
down  on  the  basis  of  a  least  squares  fit  to  the  form 

Ai  =  ai  t  +  bx  z 
Aa  -  aat  +  baz 

i.e.  strictly  linear  translation  of  the  center  of  the  correla¬ 
tion  pattern  with  time  anu  linear  tilt  with  height. 

(ax,a8)  =  Components  of  horizontal  translation  with 
time  (n  mi  per  day) 

(bj,  ,ba  )  =  Component  of  tilt  with  height  (n  rni  per 
kilometer) 

The  fit  was  made  in  two  ways,  first  for  all  five  time  lags, 
-2  to  +2  days,  then  with  only  the  three  central  time  lags,  -1 
to  +1  days.  The  speed  of  translation  is  larger  in  this  last 


Vi-45 


co  meo 

0v=t c- 

4-  OVGV 

vo  (T\cn 

cococo 

t-o  o 

tit 

rH  COO 

4*  104* 

1  1  1 

CM  t>~V£j> 

O  rH  o 

I  i  i 

Loot" 

1  rH  1 

mate  in 

1 

immo 

oco  ts- 

m  i — i  ov 

O  C4 

mvo  m 

4-  o  m 

!!=>  -P 

cm  imco 

imt- 4- 

i — 1  LOCO 
rH  OJ  rH 

oo  4-  cd 
com  co 

OJ  rH  OV 

H  (OH 

O  O'.  4 

CO  COOJ 

OJVO  co 

H  OJ  rH 

nrlo 

Error 

rH  OJ  1—1 

H  OJVO 

1 — 1 

OVrl  OV 

OJ  CO  rH 

l—VO  .H 
OOOJ  OJ 

m  oj  co 
COOJ  CO 

OH(T\ 

H  OJ  H 
rH 

t-VD4- 

m  m 

LTV  CO  VO 

4-  LO  C — 

CO  OJ  CO 

CVI  o  cv 

mco  oj 

rH  00  rH 

H 

*h 

W 

rC 

o  loco 

i — 1  rH  ! 

f 

CTstVOJ 
rH  OJ 

mmco 

1  H  OJ 

mov.4 

CCJ 

HO\m 

Ol 

cot— vo 

i — 1  i — 1 

VI 

Eh 

N 

CO 

C —  rH  in 

! 

VO  f— CO 

ov  oco 

1  rH  1 
| 

OJC0  OJ 

i  com 

i  i 

rH  mo 

OJ  4-  CO 

i  i  i 

O  004- 

oj  com 

i  i  i 

GW  VO 
rH  i — 1 

1 

to 

*  ! 

3 

vt 

RMS 

Error 

comae 
(7\- 1  t~ 

COOVLTi 

HHfO 

rH 

-=t  in  cj 

rH  rH  VO 

Lf\rH  O 
4-  OJ  4- 

CSV  t— CO 
rH  rH  CO 

101 

8 

62 

1 

Ol  C —  rH 

i — l  4"  c — 

OJ  t —  ov 

VO  rH  00 

oj  4-  m 

Cvi  Ol  t— 

H 

& 

SOON 

1  014 

1  1 

OJ  VO  4- 
CO  1  Ol 

Ol  oco 

( — 1  rH  CO 

1  1  t 

COOWD 

:  i  oi 

1 

comoj 

rH  1  i — 1 

1 

COVO  rH 
(O  H 

1 

r- 

03 

coco  vo 

C  OJ  o 
4"  4-  lA 

i — 1  OJ  ov 
OJVOVO 

OJ  rH  rH 

0VVO4- 
VOVO  O 
rH  OJ  CO 

co  4-  m 
t— mvo 

i — 1  rH 

4"  004- 
t'-C'-CO 
rH  L\J  OJ 

VD  CVI-i 

rH  rH  C7\ 

i— 1 

RMS 

Error 

CO  OJ  VO 
Vi)  VO  H 

I- 1  I- 1 

OV4-  VO 
4-  COOJ 

VO  rH  OJ 
V04-  m 

mO  rH 

m4-  co 

CO  OJ  co 
m4-co 

rH 

VO  VO  rH 
41  H  O 
rH  rH 

aj 

+ 

0 

rO 

-14.3 
12.8 
-12  .2 

co  mo 

•  •  • 

CO  CO 

rH  OJ 

lt»4-  m 

•  •  • 

vo  OV4- 

I  rH  CM 

O-COrH 

•  •  • 

mo  o 

rH  CO 

O  coco 

oco’  CO 
CO 

mvo  co 

o  r-ov 

CVJ  rH 

VI 

Eh 

OJ 

!  mooco 

CO 

IAH4- 

t~a\co 

:  i  i 

-3"  OJ  CO 
r— 1  rH 

I  | 

co  vo  co 

rH  COCO 

1  1  I 

COt—  CO 
rH  1  CO 

1 

OJ  OJ  rH 
rH  rH  5 

t 

be 

3 

VI 

RMS 

Error 

VO  C^-  OJ 

mcvco 

rH  rH 

4-  OJVO 

vo  me- 

l — 1 

4"  H  ov 
-4-  E-  OV 

72 

48 

102 

c— vo  o 

4- CO  CO 

^co^t 
Ov  CVI  1—1 

i — i  i — i 

cy 

r 

P 

t —  OJ  — 

00  OVD 

1  OJ  CO 

I  1 

co  [■—  c — 

•  •  • 

CO  i  1  VO 
CO  1  1 

1 

OVO  CO 

•  •  • 

COCO  OJ 
rH  1  CO 

1  1 

movo 

•  ♦  * 

covo  o 

1  cy 

m  rH  co 

•  »  • 

to  mco 

rH  !  t 

I 

LPv4-  rH 

•  •  O 

CTNC-O 
CO  r\J 

a 

1  mow 

J  COCO  CT\ 

1  COCOCO 

cooo  t- 

O  i — i  t — l 

OJ  rH  rH 

OVrH  OV 

i  1  C—  o 

rH  i  1  OJ 

t>-  ctut\ 

VO  H  H 
rH  rH 

CO  rH  c — 
VO  t— 4- 

i — 1  rH  i — 1 

COOJ  C— 

:  co  r- 

t — i 

a;— 
>  £ 

S-E 

o  o  o 
o  o  o 

rH  COLO 

o  o  o 
o  o  o 
t-H  com 

o  o  o 
o  o  o 
I-H  com 

o  o  o 
o  o  o 

rH  lOim 

o  o  o 
o  o  o 

rH  lO  LI  V 

o  o  o 
o  o  o 

rH  00  m 

1 

aJ  f 
<D  C 
CO  V 

3: 

CO 

t3 

co 

3: 

co 

f- 

a 

4* 

% 

Wash . 

Omaha 

Oakland 

VI-46 


ease  and  more  representative  of  a  speed  value,  Hhe  root  mean 
square  error  is  also  tabulated  in  each  case.  The  mean  wind 
components  are  shown  in  the  last  two  columns. 


The  association  between  the  mean  wind  and  the  values  of  a! 

and  aa  may  be  estimated  by  means  of  the  following  relations: 

at  the  100  mb  level  displacement  =  0.30  x  24-hour 
average  wind  movement 

at  the  300  mb  level  displacement  =  0.23  x  24-hour 
average  wind  movement 

at  the  500  mb  level  displacement  =  0.40  x  24-hour 
average  wind  movement  with  RMS  errors  of  68,  36,  and 
26  n  mi  respectively 

(if  the  displacement  is  related  to  24-hour  wind  movement  by  a 
two-term  regression  of  the  form  a  +  j3(ITx24),  if  x  24  =  movement 
of  wind  in  24  hours,  then 


a 

RMS 

Error 

-8  n  ml 

0.31 

68 

n  mi 

-23 

0.25 

31 

-18 

0.42 

21 

and  it  is  seen  that  the  additional  refinement  in  a  /  0  is  scarcely 
worth  while.) 

4.  Characteristics  of  Parameter  A 

Pigs.  VI -28  through  VI-33  show  the  variation  of  the 
parameter  A  with  height  lag  and  time  lag.  The  general  charac¬ 
teristics  of  this  variation  is  that  the  parameter  A  has  its 
maximum  value  near  height  and  time  lag  zero  and  decreases  as 
the  height  and  time  lags  depart  from  zero. 


I  irp  ■-^iirt!tT'»ny;p--|,'<|y'i"' 
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From .equation  (7a)  it  is  seen  that  the  value  of  A. depends 
on  i)  the  direction  and  orientation  of  the  gradient  of  the 
standard  deviation  of  height,  ii)  the  square  of  the  size  param¬ 
eter  L. 

The  direction  and  orientation  of  the  gradient  of  the  stand¬ 
ard  deviation  of  height  is  constant  with  respect  to  time  lag  so 
that  the  time  variation  of  A  should  be  a  reflection  of  that  of 
L.  Since  L  tends  to  increase  with  time  lag  the  reflected  varia¬ 
tion  with  respect  to  time  as  it  appears  in  A  would  be  such  as  to 
make  A  increase  with  time  lag.  This  is  just  the  reverse  of  the 
variation  actually  experienced.  One  can  only  conclude  that  the 
simplified  assumptions  that  enter  into  the  deviation  of  the 
estimate  of  the  form  of  equation  (6c)  are  inadequate. 

The  decrease  in  the  parameter  k  with  time  lag  is  ignored 
since  it-  would  be  expected  that  the  correlations  entering  the 
equations  carry  v/lth  them  such  a  factor  which  will  account  for 
the  degradation  of  correlation  due  to  the  various  factors  in¬ 
volved  in  the  parameter  k. 

The  form  of  the  expression  (6c)  was  derived  on  the  basis 
of  the  assumption  that  the  major  contribution  to  the  departure 
from  a  strictly  ideal  situation  consists  of  the  non-zero  gradient 
of  the  standard  deviation  of  height.  It  has  been  pointed  out 
that  the  effects  of  the  non -circularity  of  the  height  correlation 
contours  may  be  detected  in  the  structure  of  the  components  of 
the  correlation  tensor.  This  and  other  effects  such  as  other- 
departures  from  homogeneity  in  the  correlations  themselves  (the 


faot  that  parameters  in  the  description  of  the  correlation  func¬ 
tion  also  have  gradients)  may  well  contribute  to  the  anomalous 
behavior  of  parameter  A. 

5 .  Characteristics  of  Parameters  B  and  Q 

Figs.  Y 1-3^  through  VI-45  show  the  behavior  of  the  pa¬ 
rameters  B  and  C  as  functions  of  height  lag  and  time  lag.  The 
illustrations  in  this  sequence  are  ordered  like  those  for  param • 
eters  Al  and  A3  since  B  and  C  are  related  to  vector  components.; 
i.e.  illustration  for  B  and  C  are  in  successive  pairs  for  each 
season  and  elevation. 

The  parameters  B  and  C  are  related  (see  equations  (7a)  and 
(7b))  to  the  components  of  the  gradient  of  the  standard  devia¬ 
tion  of  height  and  as  such  they  should  be  independent  of  the 
time  lag.  Descriptions  of  Figs.  VI-34  through  VI-45  clearly 
indicate  a  time  lag  dependence  in  these  parameters.  On  the 
other  hand,  equations  (7a)  and  (7b)  contain  the  time  lag  depend 
ent  parameter  L  which  increases  with  time  lag;  but,  the  nature 
of  the  variation  of  B  and  C  with  time  could  hardly  be  attribute 
to  a  reflection  of  the  variation  of  the  parameter  L 

The  situation  with  regard  to  parameters  B  and  0  appears  to 
be  similar  to  that  encountered  in  the  case  of  parameter  A;  ther 
appear  to  have  been  so  many  simplifying  assumptions  made  in  the 
derivation  of  the  form  of  equation  (6d)  that  contact  with  the 
actual  situation  has  been  largely  lost. 


mr 


VI -5  8 


ABB-LICATXQN  TO  YAHIABILI1Y.  QP 


1.  Preliminary  Considerations 

The  variability  of  wind  in  the  vertical  may. be  described 
by  using  t  =  xx  =  xs  =  x{  =  x3  =  0  in  equations  (8)  and  (3b). 


Then 


-  k(x  -  f)('  -  ^)«p  [- 
r*»  -  k(i  -  p-X1  -  i^2)exp  [-  4( 


^2+ba2  + 


5*)]’ 

f)]  ’ 


-  =  K1  -  +  *£*'  +  (Bb1+Cba)z]exp  [-  +  jjr)J  > 

■  k(b  -  ^-)[a  *  -  (Bbl+Cb,)z]exp  [-  +  -L)]. 


x3  1 

Prom  the  first  pair,  it  is  evident  that  the  change  of  rlx  and  raa 
with  respect  to  height  difference,  z,  are  different  from  each 
other.  This  is  because  the  one  involves  factor  (l  -  b82z2/L2) 
while  the  other  involves  the  factor  (l  -  bx2z2/L2).  These  may 
be  simplified  to: 


rxx  =  k(l-Az3  )  ( 1  -  E32za)exp(-C2z2/2), 

r38  =  k(l-Az3)(l  -  B1az2)exp(-C2z2/2), 

r,g  =  k(l-Az3)(Ae  +  BxB3z2  +  Ez)exp(  -C2z2/2 ) , 

r3i  =  h(l-Az3)(A0  +  BxBsz2  -  Ez)exp(-C2z2/2)  . 


The  first  few  terms  of  the  series  expansions  of  these  correla¬ 
tions  become 


r 


i  i 


r 


as 


Az3  + 


Az3  + 
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If  only  the  part  -due  to  the  differentiable  part  ®f  the  wind  com¬ 
ponents  (the  part  In  parenthesis  on  the  right)  is  considered, 
then  from  Appendix  A-III  it  follows  that 


au  =  Standard  deviation  of  u  component 
ct'  =  Rate  of  change  of  cu  with  height 
ou ,  =  Standard  deviation  of  the  component  of  wind  shear. 


This  Implies,  since  the  right  sides  are  unequal,  that  there  are 
basic  differences  in  the  behavior  of  either  (i)  the  standard  de¬ 
viations  of  the  u  and  v  wind  components  as  a  function  of  height, 
(ii)  the  standard  deviations  of  the  wind  shear  components,  or 
(ill)  both.  The  inequality  of  the  right-hand  terms  is,  of  e curse 
brought  about  by  the  tilt  of  the  wind  correlation  patterns  in  the 
horizontal  (isobaric  surface)  as  a  function  of  height  dt  ff- rer.cf  s 
It  is  not  necessary  to  begin  with  the  specific  empirical 
expression  employed  in  Chapter  III  to  obtain  these  results. 
Consider  a  three-dimensional  differentiable  function.  Then  lei 
r(  5  J  hi  C  )  =  ao*C/  0  (  C  )  (  ?  “  biC)  +  2 &u(C)(?  “  blO(t|  -  t8 


+  ao a ( C ) (h  -  baC)2]  +  — 

where  the  coefficients  are  functions  of  Q,  the  height  difference. 
The  parameter  £  is  introduced  in  such  a  way  that  the  translation 
of  the  center  of  correlation  with  height  difference  is  explicit 
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and  Is  located  at  |  -  'biC  ~  Q>  n-baC  B  0*  Then  rearranging 
terms  to  appear  as  ascending  powers  of  £,  expanding 

a0(c)  =  1  +  a: Cs/2 :  +  — 

and  setting  §  =  r\  =  0 

r(0,0;C)  =  1  +  [a«  +  (a2obi2  +  2anbiba  +  aoaba  2  )]c2/2 1  + 

The  quantities  a",  a20,  a0 a ,  (and  perhaps  alx)  are  negative. 

The  coefficient  of  Q2  is  made  up  of  two  parts.  The  first,  a^' , 
is  the  way  in  which  the  correlation  coefficient  maximum  decays 
with  Q.  The  second,  the  parenthetical  terms,  is  the  part  due  to 
the  shift  of  the  correlation  maximum  with  height  difference. 

Consider  the  empirical  model  used  previously,  with  ai'  = 
-l/(2D2).  When  the  correlation  concerned  is  rxl,  aao  =  l/(2Ls), 
an  =  0,  a.oa  =  3/(2 L2)  while  for  raa,  aa0  =  3/(2Ls),  an  =  0, 
aos  =  l/(2La).  The  reason  for  the  difference  in  behavior  lies 
not  so  much  in  these  specific  values  as  in  the  fact  that  the 
values  are  different  from  rlx  to  ra8 . 

2 .  Practical  Analysis  of  Wind  Correlations  in  the  Vertical 
The  preceding  outline  is  unsulted  for  the  detailed  a- 
r.alysis  of  wind  component  correlations  in  the  vertical  since  it 
requires  many  levels  above  and  below  the  "master  level"  to  de¬ 
termine  the  parameters  and  since  some  of  the  parameters  indicated 
as  constants  herein  turn  out  to  be  functions  of  the  level  dif¬ 
ference,  (see  the  earlier  part  of  this  Chapter  in  which  these 
variations  were  discussed  extensively).  It  is,  however,  possible 
to  use  a  series  approximation  for  the  correlation  in  the  vertical 
over  a  limited  range  of  level  differences  and  to  obtain  specific 


rr^rf^'iri'^nir p  rr^.n 
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information  on  the  coefficients  of  the  first  and  second  terras. 

Thus,  let '  — 

r1]L  =  a0  +  a8z2  +  a3z3  +  a4z4  + - 

* 

r83  “  b»  +  t>g z2  +  b3za  +  b4z4  + 

Then  if  the  series  forms  are  terminated  after  the  fourth  degree 
term,  reasonable  estimates  of  a0 ,  a8 ,  b0,  ba  may  be  obtained 
(also  a3  and  b3 ,  but  the  meaning  of  these  terms  is  not  as  useful 
as  the  first  pairs). 

This  procedure  was  carried  out  for  the  summer  and  winter- 

wind  component  correlations  at  Washington,  D.C.,  and  Wiesbaden, 

(M 

Germany.'  '  The  results  of  the  determination  of  the  coefficients 
are  given  in  Table  II.  (in  this  table  the  values  pertaining  to 
the  north-south  wind  component  are  tabulated  under  column  headings 
that  may  make  it  appear  that  the  east-west  component  is  intended. 

The  component  intended  is  given  on  the  line  indicating  station 
and  season.)  In  the  table  the  height  of  the  "master  level"  is 
indicated  in  kilometers.  The  correlation  values  at  the  master 
level  (z  =  o)  are,  of  course,  1  and  are  not  used.  The  data  used 
were  the  correlations  for  height  lags  of  +  2,  +4,  and  +  €  km, 
the  three  values  above  and  below  the  master  level.  Values  of  the 
four  coefficients  and  the  root  mean  square  value  of  the  residuals, 
■£*,  are  tabulated.  The  coefficients  are  in  terms  of  a  2  km 
height  interval  so  that  the  units  are  (2  km)-s ,  (2  km)"3 ,  and 
(2  km)-4  for  aa ,  a3 ,  and  a4  .  The  values  of  a0  and  JP*  are  dimen¬ 
sionless  . 

The  value  of  a0  is  expressed  in  terms  of  the  standard  d-.via- 
tion  of  wind  component  and  the  standard  deviation  of  the  observation 
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errors  and  the  component  due  to  small  scale  addles «  Thus. 


k  = 


J  *  -  cr« 


where?  o*  is  the  observed  standard  deviation  of  wind  component, 
a9  is  the  standard  deviation  of  the  part  due  to  observation 
errors  and  small  scale  eddies,  ct„  is  the  standard  deviation  of 
the  wind  component  in  the  absence  of  these  effects 

°*2  =  ao2  + 

The  coefficient  of  the  second  term  of  the  series  is  (see 
Appendix  A-ITl)  given  by 


(&l)S  f  n  2(1^ 

\az  /  -°  \  a0  / 

^  dz  / 

where  is  the  standard  deviation  of  the  wind  shear  component. 

In  this  instance,  the  wind  shear  concerned  is  that  due  to  the 
differentiable  part  of  the  wind. 

To  discuss  meaningfully  the  wind  shear  vector,  it  is  assumed 
that  the  wind  shear  exists  in  the  derivative  sense,  du/dz  and 
3v/dz .  This  implies  that  the  wind  component  is  a  differentiable 
function  of  height.  The  correlation  coefficient  of  such  a  func¬ 
tion  must  be  quadratic  near  zero  height  separation,  and  hence 
the  form  of  the  series  approximation  to  the  correlation  coeffi¬ 
cients  must  omit  the  linear  term. 

The  values  of  o0  and  ae  were  determined  from  o*  and  a„  and 
are  tabulated  In  Table  II  and  Illustrated  in  Fig.  VI-46.  The 
shape  of  the  curves  of  ov.  (and  c0 )  are  well  known  and  deserve 
no  particular  comment.  The  shape  of  the  curves  of  c?e  pose  some 
rather  interesting  problems. 


Fig.  VI-46  Variation  of  c*,  aQ ,  and  ae  in  the  Vertical 
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In  all  four  eases  there  Is  a  slight  maximum  in  a6  just 
above  the  level  of  the  maximum  a0  or  d*.  This  is  even  true 
in  the  case  of  Washington,  winter,  where  the  behavior  of  ce 
at.  the  1.6  to  20  km  levels  is  anomalous.  The  value  of  <?e  is 
made  up  of  two  parts, 

a93  =  af  +  of 
i  3 

where  ot  may  be  thought  of  as  the  part  due  to  the  error  of 

i 

observation  and  oB  may  be  thought  of  as  that  part  due  to  small 

a 

scale  eddies.  It  would  seem  reasonable  to  assume  that  the  part 
due  to  observation  error  is  either  constant  or  slightly  increas 
Ing  in  a  monotonic  sense  with  altitude.  The  value  of  such  a 
function  would  appear  to  be  bounded  by  the  smaller  values  of 
cB  just  above  the  maximum  of  n0  to  be  explained  mostly  by  small 
scale  eddies. 

If  a  reasonable  value  of  o,  is  taken  as  about  5  kts,  then 

l 

to  explain  values  of  cr,  (total)  of  near  10  kts  requires  that 

the  eddy  part,  ae  ,  have  a  value  in  the  neighborhood  of  8  kts. 

a 

The  correlation  values  near  the  20  km  level  for  Washington 
are  extremely  low  even  for  a  separation  of  only  2  km  between 
levels.  This  shows  up  in  the  values  of  a0  which  decrease  to 
near  0.4  at  this  altitude.  The  implication  is  that  at  these 
levels  the  larger  part  of  the  observed  variability  of  wind  com¬ 
ponent  is  made  up  of  observation  error  and  the  effects  of  the 
small  scale  eddy  components.  At  these  levels,  the  correlation 
data  showed  some  pronounced  irregularities  so  that  the  values 
of  as  turned  out  to  be  positive.  It  is  presumed  that  such 
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v-a-l-uss -are  due-  to  the  statistical-  ^luctuatdApns-uf— tha  -data - 

with  decreasing  number  of  observations  at  these  altitudes. 

The  value  of  the  standard  deviation  of  the  wind  shear 
component  was  calculated  using  the  relation  above  in  terms  of 
the  standard  deviation  of  wind  component  and  its  vertical  gra¬ 
dient  and  the  value  of  aa  (or  b8 )  from  an  empirical  fit  of  the 
correlation  coefficient  values.  The  calculated  values  of  the 
standard  deviation  of  wind  shear  are  shown  in  Fig.  VI-47  and  in 
Table  II.  The  numerical  values  are  remarkably  consistent  show¬ 
ing  values  near  5  kts/km  from  8  to  14  km  and  decreasing  some¬ 
what  above  the  14  km  level.  (The  full  curves  in  Fig.  VI-47 
represent  values  for  the  east  pointing  component  while  the 
dashed  curves  are  for  the  north  pointing  component.) 

The  relative  contribution  of  the  two  terms  that  constitute 
the  expression  for  the  standard  deviation  of  the  wind  shear 
varies  markedly  with  the  level  concerned.  The  term  in  the  ver¬ 
tical  gradient  of  the  standard  deviation  of  height  contributes 
the  greater  part  at  levels  above  and  below  the  level  of  maximum 
standard  deviation  of  wind  component.  Near  the  level  of  maxi¬ 
mum  standard  deviation  of  wind  component  this  term  becomes  small, 
but  through  this  level  the  value  of  a8  (or  b8 )  remains  nearly 
constant  while  its  multiplier  (a02)  is  a  maximum  so  that  as  a 
result  the  value  of  cx  remains  nearly  constant. 

At  the  18  and  20  km  levels  the  values  of  are  somewhat 
uncertain.  The  correlation  coefficient  values  near  these  levels 
were  low  and  showed  some  peculiar  variations  which  caused  posi¬ 
tive  values  of  a8  (or  b8 )  to  appear  in  some  instances. 
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Fig.  VI-47  Variation  of  a  in  the  Vertical 


The  preceding  analysis  is  carried  out  or.,  the  assumption 
that  the  vertical  structure  of  the  wind  field  is  made  up  of  a 
part  that  is  differentiable  on  the  scale  of  the  spacing  of  the 
correlation  values  (2  km)  and  that  the- remaining  part  is  divided 
into  independent  random  variations,  one  of  these  due  to  observa¬ 
tion  error  and  another  due  to  small  scale  eddies  with  their 
velocity  components  considered  as  random  if  separated  by  as  much 
as  1  km. 

The  analysis  may  be  pushed  a  little  farther  on  a  conjectured 
basis  by  an  inspection  of  the  bounds  allowed  to  the  wind  shear 
due  to  the  small  scale  eddies.  If  these  eddies  can  be  considered 
homogeneous  in  the  vertical,  i.e.  the  vertical  gradient  of  the 
eddy  velocity  component  standard  deviation  is  negligible,  then 

°i2  ^  c08(-2a3) 

where  now  the  Cj  and  ae  refer  to  the  eddy  components.  If  the 
eddy  velocities  are  differentiable  in  the  vertical  (and  we  must 
stick  to  this  assumption  or  else  give  up  the  concept  of  a  v/ind 
component  shear  with  a  finite  second  moment),  then  the  eddy 
component  correlation  in  the  vertical  must  be  of  the  form 

r  =  1  -  +  — - 

2  A2 

where  A  is  a  depth  parameter  and  aa  =  l/(2A2).  Then 


CTj  (eddy)  =2  c0  (eddy)/A(eddy) , 

Ci  2:  o0/A. 

Under  our  previous  assumptions,  it  was  found  that  cr0(eddy)  —  8 
kts.  Since  it  was  assumed  that  the  eddy  velocity  components 
were  essentially  independent  at  2  km  height  separation,  the 


correlation  coefficient,  for  a  =  .2  . must  1  g  reas.Qnatly-.t-mal^  __ 
This  would  indicate  a  value  of  7\  of  not  more  than  about  1  km. 

The  resulting  value  of  at  (eddy)  is  then  approximately  8  kts/km. 
This  would  lead  to  an  estimate  of  the  total  standard  deviation 
of  the  wind  component  shear  of  near  10  kts/km. 

The  preceding  estimates  are,  of  course,  conjectured.  Until 
such  time  as  data  is  available,  one  cannot  truly  say  what  the 
small  scale  eddy  structure  is  in  the  vertical.  It  is  not  beyond 
possibility  that  it  is  indeed  not  differentiable  and  that  a 
standard  deviation  of  the  wind  shear  simply  does  not  exist. 
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*  *  *  FART  THREE  *  *  * 

CHAPTER  VII 

WIND  VARIABILITY  AT  SHORTER  DISTANCES 


A.  INTRODUCTION 

The  data  from  rawinsonde  stations  in  the  synoptic  network 
provides  information  on  the  components  of  the  wind  correlation 
tensor  over  a  range  of  distances  from  near  200  n  mi  upward. 

The  effective  cutoff  at  the  long  distance  end  is  near  3000  n 
mi,  an  arbitrary  figure  depending  primarily  on  the  collection 
or  data  available  for  analysis.  The  lower  figure  is  set  by  the 
spacing  of  the  rawinsonde  stations  in  the  synoptic  network. 
These  limitations,  especially  the  lower  distance,  implies  that 
the  structure  of  the  wind  correlation  tensor  components  is  de¬ 
termined  by  the  "synoptic"  scale  disturbances  of  the  wind  field 
Since  data  for  separation  of  less  than  200  n  mi  has  not  been 
used,  the  way  in  which  the  correlation  coefficient  rll  and  r88 


or  r^  and  rtt  approach  1  as  separation  approaches  zero  is  left 
to  conjecture  due  to  the  absence  of  data.  The  simplest  w ay  out 


of  this  situation  is  to  assume  that  the  winds  are  subject  to 


observational  error  and  small  scale  disturbances  that  are  un- 


correlatea  with  the  true  value  of  the  wind  components.  This 
assumption  leads  to  a  form  of  the  correlation  coefficient  which 
is  discontinuous  at  zero  distance.  It  is  undoubtedly  true  that 
the  rawind  observations  are  subject  to  observational  error  so 
that  the  correlations  computed  from  them  should  be  diseontinuou 


at  zero  distance.  The  smaller  scale  disturbances,  however* 
also  have  a  definite  correlation  structure  of  their  own.  ff 
there  were  data  available  on  winds  at  distances  down  to  zero, 
the  principle  components  of  the  wind  correlation  tensor  (nor¬ 
malised  form)  should  display  a  transition  zone  from  the  struc¬ 
ture  of  the  synoptic  scale  disturbances  to  a  value  nearer  to 
unity  at  zero  without  so  large  a  discontinuity.  The  situation 
is  illustrated  in  Pig.  VII-1.  If  winds  at  short  distance  sepa¬ 
rations  are  determined  by  an  electronic  system  such  as  a  doppler 
navigator  which  has  an  appreciable  lag  constant,  a  false  approach 
curve  for  small  distances  may  be  obtained  due  to  an  entirely 
different  kind  of  effect  that  might  be  classed  as  another  kind 
of  observation  error. 

The  data  presented  in  this  Chapter  is  divided  into  three 
major  parts.  In  Section  B  data  on  wind  correlation  tensor  com¬ 
ponents  is  presented  that  were  obtained  from  closely  spaced  ob¬ 
servations  of  the  conventional  synoptic  type.  In  Section  C 
some  data  is  presented  which  were  obtained  from  flights  on  which 
doppler  navigator  winds  were  read  at  five-minute  intervals.  In 
Section  D  the  wind  data  was  obtained  from  doppler  navigator 
readings  of  Project  Jet  Stream  flights.  Section  E  is  devoted 
to  a  discussion  of  the  combined  picture  of  data  from  all  sources. 
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Some  possible  correlation  structures  near  zero  diato'-.ce . 

A:  The  extrapolation  assuming  lumped  errors  end  smaller 
scale  eddies j  E:  Tire  extrapolation  assuming  random  errors 
but  accounting  for  the  structure  of  smaller  scale  eddies; 

C:  A  possible  situation  if  there  are  no  random  errors; 

D:  A  possible  case  in  which  correlations  are  computed  from 
wind  data  determined  by  a  linear  system  with  an  appreciable 
lag  constant  (over  smoothing  of  winds). 
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The  location  of  the  station  pairs,  their  separation,  and 
the  years  in  which  the  observations  were  made  .are  given  in  the 
following  table: 


TABLE  I 

STATION  PAIRS,  SEPARATION,  AND  DATA  PERIOD 


Station  Pair 

Distance 

Period  of  Data 

Port land -Brunswick 

21  n  mi 

1953-1955 

Washington-Dahlgren 

30 

1953-1957 

Rome -A lbany 

Y6 

1950-1951 

Merced -Oakland 

81 

1955-1957 

Norfolk-Hatteras 

103 

1955-1958 

One  of  each  station  pair  is  obviously  one  of  the  regular  synoptic 
network.  The  other  station  made  upper  wind  observations  for 
other  purposes.  As  will  be  noted  from  Table  II,  the  number  of 
observation  pairs  available,  the  number  of  simultaneous  wind 
observations  is  far  from  the  maximum  possible.  It  is  only  in 
the  case  of  the  last  two  station  pairs  that  there  is  reasonable 
assurance  that  the  winds  at  both  locations  are  determined  by 
electronic  tracking  methods. 

The  correlation  coefficients  and  supplementary  data  calcu¬ 
lated  from  the  data  are  listed  in  Table  II.  Figs.  VII -2  through 
VII -5  show  the  values  of  the  longitudinal  and  transverse  corre¬ 
lations  as  a  function  of  observation  point  separation.  The  data 
points  shown  in  black  are  those  from  Table  II.  The  remaining 
data  points  are  from  the  data  with  Washington  as  master  station 
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..  ...  _  - -  - -TAB-LB  II 

WIND  CORRELATIONS  -DOR 
CLOSELY  SPACED  STATION'S  AND  RELATED  DATA 


Station 

Brunswick 

Brunswick 

Washington 

Washington 

Rome 

Pair 

Portland 

Portland 

Dahlgren 

Dahlgren 

Albany 

Season 

Summer 

Winter 

Summer 

Winter 

Summer 

Level 

700  mb 

700  mb 

700  mb 

700  mb 

700  mb 

lie  (mps) 

4.3 

0.4 

3.0 

0.5 

7.5 

ut 

-6.6 

-10.8 

7.0 

14.5 

0.8 

fi; 

3.0 

1.5 

2.7 

0.9 

5.8 

V 

-7.5 

-10.1 

6.4 

14.6 

0.0 

e 

6.3 

8.4 

5.8 

9.5 

4.8 

4.7 

8.5 

4.6 

7.5 

5.4 

Vl\ 

5.8 

7.6 

5.6 

9.6 

4.3 

V 

5.2 

7.5 

4.2 

6.7 

4.8 

ru 

0.911 

0.899 

0 . 863 

0.867 

0,751 

rn 

-0.267 

0.284 

0.069 

0.001 

0.123 

rtx 

-0.048 

0.140 

-O.I38 

-0.055 

-0.121 

rtt 

0.755 

O.885 

0.855 

0.817 

0.763 

rxt  (at  p)N 

r»(9/  (of  P  *  ^ 

-0 . 092 

0.287 

-O.O37 

0.060 

-0.014 

1  it  U  r  / 

-0. 166 

0.139 

0.075 

0.018 

0.017 

N 

50 

36 

44 

46 

73 

Station 

Rome 

Oakland 

Oakland 

Oakland 

Oakland 

Pair 

Albany 

Merced 

Merced 

Merced 

Merced 

Seas  on 

Winter 

Summer 

Winter 

Summer 

Winter 

Level 

700  mb 

700  mb 

700  mb 

500  mb 

500  mb 

He  (mps) 

9.8 

-0.4 

-5.3 

-5.2 

-14.1 

Ut, 

1.9 

-2.3 

-0.8 

-4.2 

-1.2 

Ue 

9-3 

-2.2 

-6.8 

-5.4 

-13.6 

fit' 

-1.0 

-1.0 

-1.5 

-3.7 

-1.0 

6.8 

3.6 

6.6 

6.6 

11.2 

11.2 

3-7 

9.0 

6.0 

13.1 

8.8 

4.0 

6.9 

6.7 

11.2 

•V 

10.3 

4.7 

10.5 

7.0 

13-4 

r.e;e 

0.596 

0.802 

0.84? 

. .  0,900 

O.9IT 

rit 

0.010 

-0.337 

0.107 

-0.195 

0.166 

rtA 

-0.180 

-O.271 

0.077 

-0.134 

0.150 

rtt 

0.821 

0,686 

0.882 

0.872 

0.907 

r^t  (at  P)v 

0.095 

-0.4  07 

0.008 

-0.173 

0.187 

r,t  (at  P') 

-0.075 

-O.36O 

0.172 

-0.164 

0.158 

N 

21 

91 

173 

89 

167 

VII- -6 


TABLE  H 

(continued) 


Station 

Norfolk 

Norfolk 

Norfolk 

Norfolk 

Pair 

Hatteras 

Hatteras 

Hatteras 

Hattejras 

Season 

Summer 

Winter 

Summer 

Winter 

Level 

700  mb 

700  mb 

500  mb 

500  mb 

U 4  (mps ) 

3.5 

4.8 

5.4 

7.4 

ut 

4.1 

14 .0 

6.7 

19-5 

St' 

2.3 

4.2 

3.7 

6.8 

V 

2.9 

14.5 

4.9 

19.0 

tfe  (mPs) 

5.2 

8.5 

6.2 

13.3 

5.7 

7.9 

6.9 

16 . 8 

4.3 

7.7 

5.7 

12.6 

<v 

5.1 

7.5 

6.2 

15.1 

vu 

0.771 

0.842 

0 . 860 

0.926 

rjt 

-0.077 

-0.203 

-0.051 

-O.323 

rt  l 

-0.139 

-0.258 

-0.151 

-0.373 

rtt 

0.755 

0.723 

0.819 

0.920 

ru  (at  p). 

-0.054 

-0.282 

-0.076 

-0.348 

rn  (at  P') 

-0.143 

-0.250 

-0.163 

-0.370 

N 

176 

250 

174 

208 

and  at  the  500  mb  level.  The  data  points  from  Table  II  are 
labeled  with  the  first  digit  of  the  level  to  which  the  values 
pertain. 

The  spread  of  the  data  points  precludes  any  serious  analysis 
of  the  data  to  determine  the  shape  of  the  curve  near  zero  sepa¬ 
ration.  This  is  primarily  due  to  the  unreliability  of  the  cor¬ 
relation  values  which  occurs  for  at  least  two  reasons:  First, 
the  number  of  cases  available  was  small  in  several  instances; 
second,  the  observations  were  made  using  several  kinds  of  equip¬ 
ment  and  hence  the  effect  of  observation  errors  differ  from  one 
station  pair  to  another.  There  is  no  readily  available  North 
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Fig.  VII-3  r, .  as  a  Function  of  Distance,  Summer  Data 


Fie.  VII-4  r  as  a  Function  of  Distance,  Winter  Data 
0  6  6 
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American  data  available  at  700  mb  with  which  -to  -compare  t-he  five 
values  of  the  component  correlation  coefficients  at  this  level . 
The  values  at  both  500  mb  and  700  mb  are  shown  in  Figs,  VII -2 
through  VII -5  compared  with  the  data  at  500  mb  with  Washington 
as  the  master  station. 


To  show  at  least  part  of  the  difficulties  in  a  somewhat 
different  way,  all  of  the  correlations  of  Figs.  VII -2  through 
VII-5  are  shown  in  Fig.  VII -6.  In  this  case  the  values  of 
log(l-r  )  and  log(l-r^.^.)  are  used  as  ordinate  and  the  value 
of  log  (r)  (r  =  distance)  is  the  abscissa.  This  method  of  plot¬ 
ting  the  data  greatly  enlarges  the  region  near  zero  distance. 

The  data  points  lying  to  the  left  of  110  n  mi  distance  are  those 
tabulated  in  Table  II.  Superimposed  on  the  data  points  are  the 


curves  for 


v =  k  exp(-r2/2L2) 

and 

r  =  k(l  -  ra/L2)exp(-r2/2L2) 

1/  \j 

where  L  =  525  n  mi  and  values  of  k  =  1.0,  0.9,  and  0.8.  One 
would  certainly  hesitate  to  estimate  the  value  of  k  from  such 
a  spread  of  data  points,  but  the  general  behavior  of  the  points 
with  decreasing  values  of  distance  appears  to  be  of  this  type 
and  certainly  k  <  1.0. 


The  Fourth  Weather  Wing,*  in  connection  with  a  problem 
quite  different  from  that  at  hand,  obtained  data  on  wind  direc¬ 
tion  and  velocity  during  flights  near  25,000  feet  over  the 
northern  United  States,  The  data  was  obtained  by  recording 
the  direction  and  velocity  as  indicated  by  the  doppler  navigator 
equipment.  The  meteorological  observer  tabulated  these  data  at 
five -minute  intervals.  The  flights  concerned  were  as  follows: 

Route  1:  19  Jan  60,  1  Feb  60,  21  March  60 

Route  2 :  10  Jan  60,  2  Feb  60,  23  March  60 

Route  3 :  21  Jan  60,  3  Feb  60,  26  March  60. 

The  routes  flov/n  were: 

Route  1:  From  near  San  Francisco  to  Seattle  and  thence 
across  the  northern  states  to  near  Fargo  and 
thence  to  Dayton. 

Route  2:  From  Dayton  Ohio  to  Kylertown,  Pa.,  Torento  Ont . , 
Quebec,  Que.,  Caribou,  Maine,  then  south  and 
southeast  (passing  southeast  of  Cape  Cod)  to 
Norfolk,  Va . ,  and  return  to  Dayton,  Ohio. 

Route  3‘  Same  as  Route  1  only  flown  In  the  reverse  direc¬ 
tion  . 

The  total  data  recorded  consisted  of  time,  true  air  speed, 
ground  speed,  altitude,  wind  speed  and  direction,  and  latitude 
and  longitude. 

^Abbreviated  in  the  following  as  4WW. 

^We  wish  to  express  our  great  appreciation  of  the  kindness  of  4WW 
personnel  in  providing  us  with  a  copy  of  these  data  and  in  answer¬ 
ing  many  questions  on  this  program. 
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1.  Check  .of  the  Bata  Qualify. 

The  4WW  data  was  manually  recorded  by  the  meteorologi¬ 
cal  observer.  For  this  reason,  the  tabulations  were  Inspected 
to  determine  the  extent  to  which  the  data  were  affected  by  human 
and  other  errors . 

The  first  check  was  on  the  accuracy  of  the  last  digit  of 
the  recorded  wind  speed.  Wind  speed  was  tabulated  to  the  nearest 
whole  knot.  It  seemed  reasonable  to  assume  that  the  last  digit 
of  the  recorded  wind  speed  would  be  randomly  distributed  among 

the  digits  0,  1,  ,  8,  9.  Three  examples  of  the  distribution 

of  the  last  digit  are  shown  in  Fig.  VII-7.  In  the  example  for 
2  Feb  60  the  digit  count  seems  to  be  reasonably  uniform.  From 
example  for  21  Jan  60  one  would  conclude  that  the  observer  had 
a  strong  tendency  to  read  velocity  to  the  nearest  ten  knots. 

On  26  March  60  the  digit  distribution  indicates  a  tendency  to 
read  the  instrument  to  the  nearest  five  knots.  Since  the  digit 
distribution  differed  radically  from  flight  to  flight,  it  is  im¬ 
possible  to  reach  any  firm  conclusion  on  errors  in  instrument 
reading.  It  is  certainly  evident  that  an  accuracy  of  one  knot 
in  wind  speed  is  much  too  low.  An  error  of  at  least  two  or 
three  knots  due  to  instrument  reading  alone  would  seem  more 
appropriate . 

In  "eye -balling"  the  4WW  data,  a  few  cases  were  found  in 
which  the  recorded  data  seemed  inconsistent.  That  is,  the  wind 
speed  and  direction  did  not  tally  with  the  recorded  values  of 
true  air  speed,  ground  speed,  and  distance  and  direction  flown. 
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It  was  felt  that  such  items  of  inconsistent  -4a-ta-  sh-ould-not  he 
used,  To  check  the  data  self -consistency,,  the  true  .air  speed 
was  calculated  for  each  five-minute  Interval  using  the  co-ordi¬ 
nate  at  the  beginning  and  end  of  this  interval  and  the  winds 
recorded.  The  calculated  true  air  speed  was  then  compared  with 
the  recorded  true  air  speed.  The  results  are  shown  in  the 
following  Table  III: 


TABLE  III 

TRUE  AIR  SPEED  DISCREPANCIES 


Day 

Mean  Difference 

Std.  Deviation 

19  Jan  60 

0.02  kts 

5.8  kts 

20  Jan  60 

0.26 

4.0 

21  Jan  60 

0.21 

3.6 

1  Feb  60 

-2.33 

3.3 

2  Feb  60 

1.70 

3  Feb  60 

-4.63 

5.4 

21  March  60 

-0.40 

1.8 

23  March  60 

-1.54 

4.3 

26  March  60 

1.34 

3.3 

There  appeared  to  be  no  connection  between  the  days  on 
which  the  wind  speed  records  appeared  to  be  biased  toward  "o" 
in  the  last  digit  and  the  days  with  large  standard  deviation  of 
the  true  air  speed  discrepancy. 

In  this  form  of  a  check  on  self-consistent  data,  Inconsist¬ 
encies  arise  not  only  from  Incorrect  winds,  but  many  other  items 
of  data  involved.  Specifically,  the  computation  of  true  air 
speed  was  subject  to  the  errors  that  had  been  made,  not  only  In 
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wind  speed  and  direction,  but  also  .in  ground,  speed  and  track .  ... 

No  attempt  to  correct  discrepancies  is  possible.  It  Is  also 
quite  possible  that  the  check  may  reveal  a  self-consistent  situ¬ 
ation  when  the  data  items  are  actually  incorrect.  One  merely 
hopes  that  this  case  occurs  infrequently.  Rejection  of  the  cases 
in  which  the  data  was  found  to  be  inconsistent  does  not  neces¬ 
sarily  indicate  erroneous  wind  data.  A  total  of  33  out  of 
approximately  900  data  items  were  rejected  on  the  basis  of  in¬ 
consistency. 

2 .  Treatment  of  the  Data 

The  data  we re  processed  in  two  ways;  first  on  the  basis 
of  distance  intervals  in  terms  of  ground  distance,  and  second  in 
terms  of  time  lag  intervals  or  essentially  in  terms  of  air  dis¬ 
tance  intervals.  In  the  first  instance,  the  distance  between 
observation  pairs  was  computed  using  the  position  co-ordinates 
(latitude  and  longitude)  of  the  wind  data  pairs.  The  difference 
between  longitudinal  and  transverse  wind  components  at  the  points 
was  accumulated  in  classes  according  to  the  computed  distance. 

The  class  intervals  used  are  shown  :Ln  Table  IV.  It  will  be  noted 
that  at  the  lower  distances  the  class  Intervals  overlap  somewhat. 
In  the  case  of  distances  falling  into  two  class  intervals,  the 
data  were  tabulated  in  both  intervals. 

In  the  case  of  the  distance  lags  on  the  basis  of  time  in¬ 
tervals  between  observations,  the  true  air  speed  is  used  to  con¬ 
vert  time  lag  to  distance  lag.  Tables  V  and  VI  give  estimates 
of  the  conversion  factors  involved. 
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CLASS  INTERVALS 

TABLE  IV 

BASED  ON  GROUND  DISTANCE  SEPARATION 

Nominal 

Separation 

Lower  Limit 

Upper  Limit 

15  n  mi 

12 . 0  n  mi 

18,0  n  mi 

19 

15.4 

22.5 

24 

19.7 

28.6 

30 

25.2 

35.4 

38 

32.2 

44.3 

48 

41.2 

55.4 

61 

52.8 

69.4 

77 

67.5 

87.0 

98 

86.4 

108.9 

124 

110.6 

136.4 

156 

141.6 

170.8 

198 

181.2 

213.8 

250 

232.0 

267.8 

316 

296.9 

335.4 

400 

380.0 

420.0 

TABLE  V 

AVERAGE  TIME  AIR 

SPEED  AND  ITS  STANDARD 

DEVIATION 

Flight 

Average  T.A.S. 

Std .  Dev . 

19  Jan  60 

252  lets 

3  kts 

20  Jan  60 

257 

7 

21  Jan  60 

259 

6 

1  Feb  60 

271 

3 

2  Feb  60 

269 

0 

3  Feb  60 

261 

4 

19  March  60 

270 

7 

23  March  60 

266 

5 

26  March  60 

261 

7 

Composite 

262 

8 
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TABLE  VI 

CONVERSION  OP  TIME  LAG  TO  AIR  DISTANCE  LAG  (COMPOSITE) 


Time  lag 

Distance  Lag 

Std.  Dev. 

5  min. 

22  n  ml 

1  n  mi 

10 

44 

1 

15 

65 

2 

20 

87 

3 

25 

109 

3 

30 

131 

4 

40 

175 

5 

50 

218 

7 

60 

262 

8 

80 

349 

11 

100 

436 

13 

The  actual  ground  distance  between  observation  pairs  was 
computed  In  each  instance  and  the  mean  and  standard  deviation 
of  the  ground  distance  were  tabulated  along  with  the  other  data 
in  the  case  of  point  pairs  based  on  time  lag  separation. 

The  aircraft  was  not  operating  at  exactly  the  same  altitude 
during  any  given  day,  but  altitude  changes  are  relatively  minor. 
For  the  most  part  the  aircraft  flew  near  25,000  feet.  However, 
to  be  certain  that  the  effect  of  altitude  changes  between  data 
points  was  minimized,  only  data  pairs  were  accepted  in  which 
there  was  less  than  1000  feet  altitude  difference. 

In  the  case  of  ground  distance  separation  of  the  data 
points,  the  longitudinal  and  transverse  components  of  the  wind 
were  determined  as  lying  along  and  across  the  line  joining  the 
two  points.  In  this  instance,  it  is  immaterial  whether  or  not 
the  aircraft  actually  flew  this  track  or  not.  By  far  the  majority 
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of  data  points  were  along  straight  flight  fib£h&+.  tut  -oonasionally 
there  are  bends  In  the  track  flown,  this  is  particularly  true,  in 
the  ease  of  Route  No.  2  (above). 

When  the  time  lag  (air  distance)  is  used  to  classify  data 
pairs,  the  straightness  of  the  flight  path  assumes  a  role  of  some 
importance.  Consequently,  to  assure  that  the  track  between  data 
points  was  reasonably  straight,  the  track  was  computed  at  each 
point  of  a  pair.  When  the  tracks  at  the  two  points  differed  by 
more  than  10°,  the  data  pair  was  rejected.  It  is,  of  course, 
possible  that  the  track  at  two  widely  separated  points  could  have 
been  the  same,  but  the  actual  aircraft  path  between  them  may  be 
quite  sinuous.  This  sort  of  situation  did  not  occur  in  any  in¬ 
stance  . 

3.  Results 

The  results  of  the  calculation  of  the  quantities 

V  -  (VUP*  anS  V  -  <Vuhs 

from  the  4WW  data  are  given  in  Table  VII,  for  ground  distance 
lags.  The  values  of  (U.,U. )  are  for  one  point  on  the  track  and 

X  0 

for  (U',U£)  are  for  a  point  ahead.  The  distance  between  the 
points  is  given  in  the  left-hand  columns  and  refers  to  the  nomi¬ 
nal  separation  of  the  points  which  were  grouped  in  the  class  in¬ 
tervals  of  Table  IV.  The  subscripts  i  and  t  refer  to  longitudinal 
and  transverse  components. 

The  results  from  all  nine  flights  are  shown  in  Table  VII 
together  with  a  summary  in  which  all  of  the  data  have  been  com¬ 
bined.  In  each  instance  the  number  of  data  pairs  entering  into 
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"TABLE  TIT 


VALUES  OF  MEAN  SQUARED  WIND  COMPOMT  DIFFEOTOES  FROM  4VfW  BATAi 
FIXED  GROUND  DISTANCE  IADS.  DISTANCE  LAG  IN  NAUTICAL  MILES,* 
MEAN  SQUARE  VELOCITY  DIFFERENCES  IN  KNOTS  SQUARED. 


FLIGHT: 

19 

January  60 

FLIGHT: 

20  January  60 

Dist. 

a  2 

A  2 

Dist. 

A  2 

A  2 

.  Lag 

N 

A,e 

At 

Lag 

N 

A* 

At 

15 

7 

7.7 

5-3 

15 

17 

40.3 

51.2 

19 

49 

15.3 

14.9 

19 

50 

48.5 

62.5 

24 

58 

17.8 

27.7 

24 

48 

67.7 

87.0 

30 

14 

33.4 

53.9 

30 

18 

52.7 

33.7 

38 

46 

25.5 

21.3 

38 

58 

89.6 

99.7 

48 

58 

53.5 

42 .1 

48 

55 

99.3 

165.3 

6l 

58 

53.3 

26.6 

61 

70 

127.1 

228.0 

77 

71 

99-5 

66.9 

77 

78 

162.8 

306.3 

98 

84 

129.8 

76.9 

98 

97 

i47.7 

444.1 

124 

95 

227.8 

90.8 

124 

113 

208.0 

684.1 

156 

107 

260.4 

106.7 

156 

129 

251.8 

955.4 

198 

115 

375.0 

129.3 

198 

147 

323.7 

1250.8 

250 

122 

622.2 

180.6 

250 

171 

422.7 

1614.8 

316 

122 

928.7 

205.6 

316 

102 

869.0 

2841.3 

4  00 

90 

1206.2 

301.9 

400 

56 

1482.5 

2859.6 

FLIGHT: 

21 

January 

60 

FLIGHT : 

1 

February  60 

15 

51 

76.2 

43a 

15 

8 

2.8 

13.7 

19 

73 

63.1 

34.3 

19 

44 

20.6 

14.6 

24 

33 

52.5 

19.3 

24 

55 

20.7 

16.6 

30 

47 

134.9 

84.9 

30 

6 

28.7 

17.6 

38 

72 

112.6 

71.8 

38 

41 

49.0 

35.1 

48 

79 

146.2 

99.1 

48 

43 

33.4 

30.5 

61 

91 

120.2 

142 .6 

6l 

55 

73.6 

52.3 

77 

96 

218.8 

156.4 

77 

61 

83.4 

59.8 

98 

122 

278.1 

202  .6 

98 

67 

111.2 

67.0 

124 

119 

395.2 

284.5 

124 

78 

181.1 

96.6 

156 

145 

544.3 

398.4 

156 

75 

201.0 

103.9 

198 

151 

837.6 

524.6 

198 

93 

265.2 

120.2 

250 

152 

1153.0 

581.2 

250 

97 

389.4 

138.3 

nl  Y! 

116 

1833.8 

803.7 

316 

98 

526.2 

57.9 

400 

43 

1215,0 

190.9 

400 

88 

696.6 

117.0 
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PLIGHT  t 

2  February  60 

FLIGHT; 

3  February  60 

Dist. 

A  2 

A  2 

Dist. 

A  2 

A  2 

lag 

N 

At 

Lag 

N 

At 

15 

2 

26.1 

7.5 

15 

16 

43.8 

59.7 

19 

10 

33.4 

48.7 

19 

57 

25.2 

43.1 

24 

28 

16.0 

16.9 

2  4 

45 

18.7 

51.4 

30 

25 

30.0 

22.5 

30 

*i 

C.  X 

121.1 

193.8 

38 

10 

76.6 

55.9 

38 

56 

ol.l 

74.8 

48 

31 

89.5 

38.9 

48 

48 

118.3 

193.6 

61 

26 

80.4 

75.2 

61 

70 

152.3 

184.7 

77 

38 

120.6 

136.7 

77 

68 

203.0 

269.8 

98 

30 

225.0 

182.0 

98 

83 

209.4 

352.9 

124 

33 

216.2 

192.5 

12  4 

89 

259.4 

482.2 

156 

33 

442  .2 

410.7 

156 

107 

300.5 

710.5 

198 

42 

564.2 

709.6 

198 

115 

378.6 

1119.1 

250 

53 

1055.0 

737.6 

250 

129 

500.4 

1595.0 

316 

43 

1244.3 

682.1 

316 

136 

684.8 

2287.7 

400 

31 

739.0 

331.5 

400 

68 

538.0 

1541.0 

FLIGHT: 

21 

March  60 

FLIGHT: 

23 

March  60 

15 

1 

— 

15 

27 

24.2 

20.4 

19 

14 

23.5 

26.7 

19 

49 

14.2 

21.8 

24 

63 

11.3 

17.2 

24 

34 

25.3 

30.0 

30 

33 

18.2 

39.6 

30 

45 

35.1 

78.0 

38 

8 

26.4 

81.1 

38 

54 

33.9 

73.0 

48 

61 

32.0 

38.6 

48 

55 

36.8 

144.4 

61 

31 

60.3 

56.9 

61 

69 

56.1 

193.5 

77 

65 

44.2 

45.8 

77 

79 

63.3 

289.4 

98 

64 

52.1 

57.1 

98 

86 

113.2 

376.3 

124 

72 

70.3 

79.1 

124 

99 

171.3 

548.8 

156 

83 

75.2 

83.0 

156 

102 

302.6 

649.0 

198 

95 

86.6 

102.6 

198 

120 

400.3 

665.2 

250 

98 

97.7 

131.5 

250 

129 

678.7 

867.4 

316 

101 

124.7 

175.3 

316 

74 

1076.3 

1070.6 

400 

66 

223.9 

123.2 

400 

11 

422.1 

871.2 
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TABLE  “VII 
(CONTINUED) 


FLIGHT 

;  26 

March  60 

FLIGHT: 

SUMMARY 

Dlst- . 

— y 

A  2 

Dist. 

A  2 

V 

Lag 

N 

At 

Lag 

N 

15 

69 

13.8 

27.9 

15 

198 

37.1 

34.4 

19 

82 

14.1 

22.0 

19 

428 

30.2 

31.7 

24 

20 

10.0 

50.5 

24 

384 

27.4 

36.7 

30 

71 

31.5 

53.3 

30 

280 

56.1 

69.2 

38 

87 

30.1 

42.9 

38 

432 

60.5 

63.9 

48 

90 

44.7 

57.6 

48 

520 

75.0 

94.0 

61 

90 

47.9 

60.7 

61 

560 

97.1 

127.5 

77 

114 

54.8 

60.9 

77 

670 

123.3 

171.6 

98 

142 

78.5 

59.1 

98 

775 

157.8 

227.3 

12  4 

137 

76.8 

73-5 

124 

835 

217.4 

332.7 

156 

162 

93.3 

80.6 

156 

943 

296.8 

478.8 

198 

167 

82.8 

87.5 

198 

1045 

412.2 

674.4 

250 

168 

100.6 

88.7 

250 

1119 

609.6 

906.0 

318 

148 

92.0 

91.3 

316 

940 

893.8 

1193.7 

*100 

8 

35.0 

23.2 

400 

461 

1140.6 

1175.8 

the 

calculation  of 

T*  and 

Ats  are  shown 

in  the 

column 

headed  N 

The  value  of  N  differs  widely  from  one  distance  lag  to  another 


in  Individual  flights.  This  is  due  to  the  fact  that  the  data 
were  read  at  five-minute  intervals,  but  the  aircraft,  though 
flown  at  a  reasonably  constant  true  air  speed,  was  influenced  by 
head  and  tail  winds.  Thus  under  different  wind  conditions  the 
data  values  accumulate  in  different  class  intervals.  This  effect 
is  especially  pronounced  when  data  point  separation  is  less  than 
40  n  mi.  Here  we  are  dealing  with  data  points  separated  in  time 
by  one  and  two  five -minute  intervals.  The  class  intervals  are 
also  narrow  in  this  range  so  that  the  effect  of  head  or  tail 
winds  becomes  pronounced. 
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There  is  a  tendency  for  the  number  of  points  to  dacrease  - 
abruptly  at  the  400  n  mi  class  interval.  This  is  due  to  a  com¬ 
bination  of  many  factors  which  include  the  large  time  separation 
and  the  effect  of  the  rejection  criteria  such  as  data  at  same 
altitude  for  each  of  a  pair. 

The  data  for  any  one  flight  represent  a  kind  of  space  aver¬ 
age  of  the  mean  square  value  of  the  component  differences  for 
each  day.  There  Is  a  great  variation  in  the  behavior  of  the  mean 
square  component  differences  from  one  day  to  another.  This  is 
brought  about  by  several  factors  of  which  some  are: 

i)  The  natural  variability  associated  with  the  rela¬ 
tively  few  samples  entering  into  each  individual  day's  distance 
lag  computations . 

ii)  The  variability  of  the  Inherent  wind  variability 
from  place  to  place  from  day  to  day  In  the  area  covered  by  the 
flight . 

It  is  felt  that  the  data  for  individual  days  is  without 
significance  and  that  only  the  summary  data  can  be  Interpreted 
seriously. 

The  data  for  the  same  quantities  for  air  distance  lags  is 
tabulated  in  Table  VIII.  In  this  table,  the  rows  are  indicated 
by  the  time  separation  of  the  wind  component  pairs.  The  corres¬ 
ponding  distance  separation  is  given  in  Table  VI.  The  average 
ground  distance  separation  and  Its  standard  deviation,  d  and  a4 , 
are  also  tabulated.  In  addition  to  A^2  and  A^2,  the  values  of 
and  S^T  are  also  shown.  (This  last  pair  was  not  tabulated  in 
the  case  of  the  ground  distance  lags.  In  Table  VII,  the  data 
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TOOTTTOX' 


VALUES  OF  MEAN  SQUARE  WIND  COMPONENT  - - 

DIFFERENCES  FROM  4  WW  DATA  WITH  AIR  DISTANCE  (TIME)  LACS 


Time 

7T 

/T" 

a  2 

A- s 

Lag 
(min . ) 

N 

d 

(n  ml) 

°4 

(n  ml) 

(knots)  (knots) 

X 

(knots®) 

X/ 

(knots2) 

FLIGHT: 

19  January 

60 

5 

78 

23 

6 

0.4 

-0.4 

20.1 

26.1 

10 

68 

44 

8 

0.4 

-0.9 

50.5 

39.0 

15 

63 

66 

8 

0.7 

-1.7 

96.4 

62.0 

20 

61 

87 

10 

1.1 

-2.3 

152 .8 

81.8 

°5 

57 

109 

12 

1.0 

-2.8 

228.2 

107.1 

30 

53 

131 

14 

1.3 

-3.8 

304.1 

133.9 

40 

48 

174 

18 

2  .0 

-5.7 

446.7 

180.5 

50 

44 

218 

19 

3.7 

-6 . 1 

523.7 

251.1 

60 

41 

258 

24 

3.7 

-6.5 

592.3 

238.6 

80 

34 

338 

30 

3.9 

-6.1 

618.4 

326.4 

100 

32 

415 

32 

6.4 

-7.5 

1068.0 

475.5 

FLIGHT: 

20  January 

60 

5 

84 

22 

5 

-1.1 

-1.9 

61.9 

73.8 

10 

68 

44 

8 

-1.9 

-4.1 

U9.3 

148.9 

15 

61 

64 

11 

-3.5 

-4.9 

162.3 

278.0 

20 

52 

86 

13 

-4.5 

-6.0 

175.1 

351.4 

25 

44 

107 

15 

-5.5 

-6.9 

267.5 

479.8 

30 

38 

127 

18 

-7.6 

-10.3 

336.0 

714.6 

4o 

27 

167 

20 

-11.9 

-16.8 

616.9 

1437.9 

50 

21 

216 

26 

-13.5 

-17.8 

1009.4 

1846.8 

60 

14 

248 

29 

-21.5 

-31.3 

1398.7 

3138.8 

80 

8 

306 

14 

-51.3 

-79-2 

2781.4 

7  632.4 

100 

4 

384 

5 

-50.4 

107.3 

2750.2 

11676.5 

FLIGHT: 

21  January 

■  60 

5 

97 

18 

3 

-0.6 

0.6 

63.7 

34.8 

10 

84 

38 

5 

-0.5 

0.6 

98.9 

65.8 

15 

79 

54 

7 

-1.2 

1.4 

149.2 

116.5 

20 

78 

71 

10 

-1.8 

1.2 

209.2 

152.7 

25 

75 

89 

11 

-2.1 

2  .2 

266.2 

223.0 

30 

70 

106 

14 

-1.4 

2.6 

319.1 

270.7 

40 

60 

141 

19 

-0.6 

3.3 

437.0 

374.6 

50 

52 

176 

25 

-2.7 

6.2 

628.4 

482.6 

60 

52 

211 

29 

-2  .0 

4.1 

786.5 

550.0 

8o 

43 

279 

38 

-2.4 

3-8 

1234.6 

662.7 

100 

35 

346 

46 

-2.7 

7.0 

1827.5 

888.3 
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TABLE  ¥1X1 
(CONTimJIB) 


Time 

Lag 

a- 

04 

At 

V 

Ata 

(min. ) 

JL 

(njnl). 

(rL.mi) 

(toots  5 

(toots) 

(toots2) 

(toots2) 

PLIGHT: 

1  February  60 

5 

69 

22 

4 

0.8 

0.7 

20.1 

18.1 

10 

62 

43 

5 

1.2 

1.0 

45.5 

32.8 

15 

56 

64 

6 

1.7 

1.2 

77.8 

51.0 

20 

55 

84 

8 

1.9 

1.0 

129.6 

68.7 

25 

53 

105 

9 

1.9 

0.9 

172.8 

73.1 

3° 

52 

126 

10 

2.4 

1.5 

216.6 

85.8 

4o 

48 

168 

11 

3.9 

3.0 

294.4 

111.4 

50 

45 

210 

12 

3-9 

1.1 

403.6 

92.3 

00 

40 

252 

14 

2.0 

-0.3 

541.6 

59.7 

80 

35 

334 

17 

-0.1 

3.6 

736.1 

108.3 

100 

31 

414 

18 

-3-3 

6.1 

801.2 

154.0 

FLIGHT: 

2  February  60 

5 

46 

26 

4 

1.3 

-2.1 

27.9 

25.1 

10 

32 

52 

5 

4.5 

-5.5 

76.5 

59.7 

15 

26 

78 

6 

9.4 

-7.8 

161.1 

113.9 

20 

21 

105 

7 

11.3 

-12.5 

267.9 

267.4 

25 

18 

131 

7 

11.4 

-18.3 

404.1 

481.3 

30 

16 

156 

8 

7.7 

-26.9 

484.0 

932.8 

40 

11 

211 

6 

-1.5 

-43.3 

567.0 

2065.7 

50 

10 

260 

5 

-3.3 

-55.9 

801.1 

3286.5 

60 

8 

311 

7 

-9.3 

-60.0 

539.8 

4015.0 

80 

2 

407 

0 

-10.4 

-72.5 

115.0 

5268.4 

100 

- 

- 

- 

- 

- 

- 

- 

PLIGHT: 

3  February  60 

5 

76 

20 

4 

-0.5 

0.7 

3°.3 

54.0 

10 

71 

41 

6 

-1.5 

1.3 

81.7 

135.9 

15 

65 

61 

8 

-2.7 

2.4 

133.5 

205.7 

20 

61 

81 

10 

-3.0 

3.1 

199.2 

317.4 

25 

56 

102 

12 

-6 . 0 

3.9 

248.0 

438.1 

30 

49 

122 

13 

-8.6 

8.1 

312.6 

589.4 

40 

45 

162 

17 

-13.0 

1.2 

440.5 

897.2 

50 

38 

202 

19 

-17.4 

-3.8 

580.2 

1176.9 

60 

31 

242 

21 

-22,6 

-10.7 

904.7 

1461.5 

8o 

28 

290 

86 

-25.O 

-15.4 

1291.1 

192  0.3 

100 

26 

333 

105 

-23.7 

-17.2 

1566.5 

2385.2 
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TABU!  VIII 
(CONTINUED) 


Time 

Lag 

d 

At 

0s 

V 

( min . ) 

N 

(n  mi) 

(n  mi) 

(knots) 

(knots ) 

(knots2) 

(knots2) 

PLIGHT: 

21  March 

60 

5 

75 

25 

3 

0.3 

0.8 

17.2 

23.8 

10 

66 

50 

5 

0.8 

2.1 

36.2 

41.8 

15 

61 

75 

6 

1.2 

3.6 

52.1 

51.2 

20 

53 

99 

7 

1.2 

5.0 

58.9 

77.1 

25 

47 

123 

8 

1.3 

6.3 

72.1 

121.4 

3° 

45 

147 

9 

2.1 

7.9 

83.1 

150.3 

4o 

40 

193 

14 

3-0 

9.6 

93. 1 

191.7 

50 

35 

241 

14 

4.2 

13.0 

IO8.5 

285.5 

60 

31 

288 

11 

7.0 

17.1 

122.5 

390.5 

80 

27 

383 

15 

10.3 

21.2 

255.5 

563.4 

100 

27 

476 

10 

14.5 

27.7 

353.0 

870.4 

PLIGHT: 

23  March 

60 

5 

83 

20 

5 

-0.2 

0.3 

18.9 

21.9 

10 

63 

40 

9 

-1.0 

0.5 

30.6 

67.0 

15 

54 

58 

11 

-1.0 

-0.1 

49.8 

162 .1 

20 

45 

75 

13 

-2.5 

-0.4 

64.2 

281.3 

25 

42 

95 

17 

-1.9 

-1.1 

79.6 

423.6 

30 

37 

110 

18 

-3.0 

-2.5 

126.2 

597.8 

40 

29 

151 

26 

-0.1 

-6.1 

263.2 

894.8 

50 

24 

190 

26 

1.8 

-4.9 

494  .5 

975.2 

60 

20 

225 

25 

3.4 

-7.4 

680.2 

1155.1 

80 

11 

282 

26 

-9.5 

-10,0 

883.4 

825.8 

100 

6 

320 

36 

-26.4 

-5.4 

948.0 

124.7 

PLIGHT: 

26  March 

60 

5 

105 

17 

4 

-0.1 

0.5 

13.6 

24.1 

10 

91 

34 

4 

0.0 

1.4 

30.7 

48.0 

15 

81 

52 

4 

0.4 

2.5 

41.9 

54.4 

20 

72 

69 

5 

0.3 

4.2 

70.7 

75.6 

25 

63 

86 

6 

-0.1 

5.3 

85.8 

85.0 

3° 

59 

103 

6 

-0.1 

5.2 

89.7 

91.6 

4o 

56 

139 

8 

l.l 

7.7 

96.5 

133.8 

50 

47 

173 

9 

-0.3 

11.5 

H5.5 

216.5 

60 

43 

209 

10 

-0.6 

12.5 

113.7 

250.4 

80 

39 

276 

14 

1.1 

15.4 

87.6 

339.1 

100 

38 

345 

16 

0.7 

18.0 

122  .7 

367.8 
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Time 

Lag 

d 

TABLE  VIII 

(CONTINUED) 

*7 

V 

V 

(min.) 

N 

(n  mi) 

(n  mi)  (knots) 

(knots) 

(knots*) 

(knots*. 

5 

713 

21 

FLIGHT :  SUMMARY 

5  0.0  0.0 

31.1 

34.0 

10 

605 

42 

8 

-0.1 

0.0 

62.9 

71.5 

15 

546 

62 

11 

-0.1 

0.2 

99.3 

120.0 

20 

498 

82 

14 

-0.5 

0.4 

141.8 

172.8 

25 

455 

102 

17 

-0.9 

0.5 

190.9 

241.3 

3° 

419 

122 

20 

-1.2 

-0.1 

235.3 

324.3 

40 

364 

162 

25 

-i  .4 

-0.5 

331.5 

492.9 

50 

316 

203 

31 

-1.9 

-0.3 

459.2 

634.0 

60 

280 

241 

35 

-2.5 

-1.4 

570.8 

745.4 

80 

227 

312 

54 

-3.9 

-0.3 

771.0 

908.8 

100 

199 

383 

68 

-2.8 

3.8 

969.6 

989.7 

points  were  paired  together  even  though  the  flight  path  had  sharp 
corners  so  that  a  mean  value  of  the  longitudinal  or  transverse 
wind  component  difference  would  be  without  meaning.  In  Table  VIII, 
pairs  were  rejected  if  the  flight  path  between  points  was  not 
reasonably  straight,) 

The  number  of  data  pairs  for  each  flight  decreases  with  in¬ 
creasing  time  lag.  This  is  due  to  several  items,  mostly  the  ef¬ 
fect  of  rejection  criteria  of  which  the  condition  for  a  reasonably 
straight  flight  path  is  prominent.  The  effect  is  easily  seen  on 
20  Jan  60,  2  Feb  60,  and  23  March  60  when  the  flight  path  formed  a 
loop  from  Dayton  around  the  northeastern  states  and  return  (Route 
2). 

The  values  of  A^2  and  A^2  from  Table  VIII  for  each  flight 
show  the  same  variation  in  wind  behavior  as  a  function  of  lag  as 
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those  In.  the  previous  table  and  for  the  same  reasons,.  Only  the 
summary  table  in  which  values  from  all  nine  flights  are  combined 
is  considered  of  any  real  significance. 

The  composite  values  of  A^a  and  A^2  as  functions  of  distance 
lag  are  shown  graphically  in  Pig.  VII -8.  The  horizontal  scale 
is  ground  distance  in  the  case  of  the  values  determined  from  such 
class  intervals  and  is  air  distance  lag  for  the  data  determined 
from  observation  time  lags.  The  two  pairs  of  data  points  appear 
to  agree  reasonably  well.  The  agreement  would  have  been  even 
better  if  the  average  ground  distance  had  been  used  for  the  ab¬ 
scissa  for  data  grouped  on  the  basis  of  time  lags.  (Compare  dis¬ 
tances  of  Table  VI  with  the  average  distance  of  Table  VIII. 

The  central  parts  of  the  curves  (140  n  mi  to  200  n  mi)  seem 
to  be  well  behaved  as  far  as  the  mean  square  of  component  differ¬ 
ence  vs  distance  relative  is  concerned.  At  the  short  distance 
end,  the  points  for  five-minute  time  lag  (22  n  mi)  appear  reason¬ 
ably  placed  as  an  extension  of  the  remaining  points.  The  point 
pairs  from  ground  distance  lag  classes  seem  to  jump  around  un¬ 
necessarily.  It  is  easily  seen  from  Table  VII  that  the  first 
five  points  select  samples  from  the  nine  flights  in  a  very  ir¬ 
regular  manner  whereas  the  sampling  from  the  middle  group  of 
points  is  reasonably  uniform. 

At  large  distances,  the  same  nonrepresentative  sampling  may 
account  for  the  departure  of  the  data  points  from  the  trend  of 
the  middle  points. 


Square  Component  Differences  (knots*) 


Distance  (n.  mi.) 


Fig.  VII-8  Comparison  of  Value  of  Mean  Squared  Component 
Differences  Based  on  Ground  Distance  Classes 
and  Air  Distances  (time  lag)  Classes 
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D.  MEAN  SG.UABE  WIND  COMPONENT  .B-IEBiBlMCES  WM  Ea&EBfliE.  JBffi 

SmBAM  "BAf?  ~  “  ~~ 

1 .  General  Remarks 

Early  in  this  study  Mr,  I,  I.  Gringorten  sailed  our 
attention  to  the  wind  data  available  from  Project  Jet  Stream. ^ 
In  connection  with  the  Project  Jet  Stream  program,  a  series  of 
volumes  were  prepared  entitled,  "Flight  Data  Pertaining  to  Jet 
Streams,  Jet  No.  XXX,  Flown  XX  Month  Year"  where  the  X's,  month, 
and  year  are  appropriately  chosen.  These  volumes  contain  graphs 
of  the  quantities  measured  and  a  printout  of  the  data  from  which 
the  graph?1-  were  prepared,  together  with  maps  showing  the  general 
v;eather  situation,  the  flight  path  and  some  wind  data  along  it, 
and  other  pertinent  information.  The  data  were  obtained  from 
photographic  recordings  of  the  aircraft  Instrumentation  and 
were  reduced  at  the  Data  Reduction  Laboratory  of  the  University 
of  Dayton,  Dayton,  Ohio.  For  our  purposes,  we  are  interested 
in  the  wind  data  from  the  doppler  navigator  equipment. 

The  frame  speed  of  the  camera  used  to  photograph  the  in¬ 
strument  dials  was  near  two  frames  per  minute.  Though  the  air¬ 
craft  had  a  true  air  speed  of  near  423  knots  compared  with  260 
knots  for  the  4WW  data,  the  fact  that  the  data  was  recorded  ten 
times  more  frequently  makes  possible  wind  comparisons  with  a 
distance  separation  of  order  of  3  n  mi  instead  of  22  n  mi. 

Had  the  flights  all  been  made  entirely  in  the  jet  stream 
region,  there  would  be  considerable  question  of  the  validity  of 
the  use  of  this  data  to  make  inferences  concerning  the  turbulent 
structure  of  the  atmosphere  In  general  at  or  near  the  jet  stream 
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level.  Fortunately,  there  were  several  flights  that  were  re¬ 
quired  for  miscellaneous  purposes  on  which  the  instrument  arid 
recorder  were  operating  and  which  were  either  wholly  or  at  least 
extensively  outside  the  jet  stream  as  such.  Only  such  flights 
or  flight  legs  were  selected  for  inclusion  in  this  study. 

2 .  The  Data 

The  data  used  was  selected  from  the  flights  listed  in 
Table  IX.  The  fact  that  nearly  all  the  flight  legs  In  Table  IX 
were  on  a  westerly  track  was  quite  accidental.  The  subdivision 
Into  several  subsections  is  based  on  minor  changes  of  track. 

The  data  from  Project  Jet  Stream  were  handled  in  a  more 
straightforward  manner  than  the  4WW  data.  Instead  of  dealing 
with  a  flight  along  a  crooked  track,  Project  Jet  Stream  data 
were  selected  so  as  to  lie  along  essentially  straight  flight  path 
segments.  Altitude  in  each  segment  was  carefully  maintained  and 
each  segment  was  checked  for  such  constancy  of  altitude  before 
Inclusion  In  the  list  of  Table  IX. 

The  air  distances  corresponding  to  the  time  lags  used  are 
listed  in  Table  X  for  each  JETNO  (our  numbers  1  through  6). 

The  lag  number  Increases  by  one  unit  for  simple  indexing  purposes. 
The  number  of  time  lag  intervals  corresponding  to  each  lag  number 
is  shown  in  the  second  column  of  Table  X.  These  were  selected 
to  roughly  approximate  an  exponential  form  (i.e.,  that  the  loga¬ 
rithm  of  the  time  lags  would  be  roughly  uniformly  spaced).  The 
air  distances  corresponding  to  each  time  lag  were  computed  from 
the  mean  true  air  speed  and  mean  interval  of  time  lag  for  each 
JETNO.  These  distances  are  subject  to  some  dispersion  that 
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arises.  from  two  sour g as,  the  fact  .that  the  truoair  speed- .-is  net 
really  uniform. and  the  fact  that  the  tinier  actuating  the  camera 
shutter  did  not  run  exactly  uniformly.  The  largest  factor  is 
the  true  air  speed  variation,  which  (Table  IX)  amounts  to  approxi 
mately  1.0^  (standard  deviation),  (in  JETNO-2  it  amounts  to  h°/0 
and  JETNO-1  and  JETNO-4  to  about  1.7$.  The  case  of  JETNO-2  is 
the  only  situation  in  which  flight  was  in  the  jet  stream  core. 
This  is  presumed  to  account  for  the  higher  variability  of  the 
true  air  speed.) 


TABLE  X 

LAG  NUMBERS,  TIME  INTERVALS,  AND 
AIR  DISTANCES  FOR  JET  STREAM  DATA 


Lag 

No.  of 

Distance 

(n  mi) 

No. 

Intervals 

1 

2 

'  2 

4 

5 

b 

1 

1 

3.9 

3.8 

3.9 

3.1 

2.9 

0.8 

2 

2 

7.7 

7.7 

7.7 

6.2 

5.9 

1.6 

3 

3 

11.6 

11.5 

11.6 

9-3 

8.8 

2.4 

4 

4 

15.4 

15.4 

15.4 

12.4 

11.7 

3.2 

5 

6 

23.2 

23.I 

23.I 

18.6 

17.6 

4.9 

6 

8 

30.9 

30.8 

30.8 

24.7 

23.4 

6.5 

7 

10 

38.6 

38.5 

38.5 

30.9 

29.3 

8.1 

8 

12 

46.3 

46.2 

46.2 

37.1 

35.1 

9.7 

9 

16 

61.8 

61.5 

61.6 

49.5 

46.8 

12.9 

10 

20 

77.2 

76.9 

77.0 

61.9 

58.5 

16.2 

11 

24 

92.7 

92.3 

92.4 

74.2 

70.2 

19.4 

12 

28 

108 

108 

108 

86.6 

81.9 

22.7 

13 

36 

139 

138 

139 

111 

105 

29.1 

14 

44 

169 

169 

169 

136 

129 

35.6 

15 

52 

201 

200 

200 

161 

152 

42.1 

16 

60 

232 

231 

231 

186 

176 

48.5 

3.  Results  ...  - 

The  values  of  the  mean  square  wind  component  differences 
that  are  obtained  from  the  Project  Jet  Stream  observations  are 
tabulated  in  Tabic  XI.  The  tabulated  data  is  somewhat  more  exten¬ 
sive  than  that  from  the  4WW  observations  since  the  flight  paths 
were  straight  and  some  meaning  could  be  assigned  to  such  items  as 
the  mean  longitudinal  and  transverse  wind  components  and  their 
standard  deviations. 

It  is  felt  that  the  values  of  the  mean  square  w ind  veloc¬ 
ity  component  differences  from  the  individual  JETNO's  are  somewhat 
unreliable  (with  the  possible  exception  of  JETNO-5).  The  amount 
of  data  entering  into  each  is  not  only  limited  in  quantity,  but 
also  the  observations  are  made  in  a  highly  specific  atmospheric 
situation . 

Remarks  on  tne  wind  variability  in  general  are  usually  under¬ 
stood  in  the  sense  of  a  kind  of  average  of  a  wide  variety  of  situ¬ 
ations.  Consequently,  a  composite  of  the  JETNO-1  through  JETNO-5 
was  prepared  and  the  results  tabulated  in  Table  XII.  (Here  the 
mean  component  values  and  their  standard  deviation  are  omitted 
since  they  have  no  significance  in  the  composite  picture.) 

In  spite  of  the  fact  that  we  consider  the  mean  square  com¬ 
ponent  differences  of  the  individual  JETNO's  as  unreliable,  it 
is  interesting  to  compare  the  internal  variations  with  those  of 
the  variations  of  the  nine  4WW  flights.  The  values  of  A^2  and 
A^2  in  Table  VIII  vary  widely  in  this  regard;  in  some  instance 
A^2  <  A^2  holds  and  in  others  the  equality  is  (strongly)  reversed. 


TABLE  XI 

DATA  ON  MEAN  SQUARE  VELOCITY  DIFFERENCES  FROM  PROJECT  JET  STREAM  DATA 
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TABU  XII 

COMPOSITE  OP  DATA  ON  VELOCITY  COMPONENT 
DIPPERENCES  PROM  JETNO-1  THROUGH  JETNO-5 


lag 

No. 

N 

d 

1 

1074 

3-3 

0.4 

2 

1069 

6 . 6 

0.8 

3 

1064 

9.8 

1.3 

4 

1059 

13.1 

1.7 

5 

1049 

19.6 

2.5 

6 

1039 

26.2 

3.3 

7 

1029 

32.7 

4.2 

8 

1019 

39.2 

5.0 

9 

999 

52.2 

6.6 

10 

979 

65.2 

8.3 

11 

959 

78.0 

9.9 

12 

939 

91.1 

10.7 

13 

899 

116.4 

14.6 

14 

859 

141.6 

17.4 

15 

819 

166.8 

20.5 

16 

779 

191.5 

23.2 

A  . 

A. 

Aa2 

A  2 

& 

t 

i 

t 

-0.1 

0.0 

1.8 

4.0 

-0.2 

0.1 

4.1 

8.9 

-0.2 

0.1 

6.6 

14.3 

-0.4 

0.1 

8.9 

19.6 

-0.5 

0.2 

13.4 

30.0 

-0.6 

0.2 

18.5 

39.1 

-0.7 

0.2 

23.2 

47.7 

-0.9 

0.3 

27.0 

56.1 

-1.2 

0.4 

33.9 

68.8 

-1.6 

0.5 

38.9 

79.0 

-2.1 

0.6 

42.7 

92.0 

-2.5 

0.8 

47.3 

106.1 

-3.2 

1.0 

55.0 

123.8 

-4 . 4 

0.9 

58.1 

142.3 

-5.5 

0.4 

55.8 

149.7 

-6.7 

-0.2 

58.2 

150.9 

This  inequality  is  much  more  consistently  valid  in  the  data  of 
Table  XI.  It  is  not  known  whether  this  Is  a  result  of  the  dif¬ 
ference  due  to  different  eddy  structures  at  25,000  and  35>000 
feet. 

Columns  for  the  mean  air  distance  between  points  (d)  and 
its  standard  deviation  (crt )  are  shown  in  Table  XII.  The  values 
of  mean  air  distance  are  the  weighted  mean  values  for  the  indi¬ 
vidual  JETNO's.  The  standard  deviation  tabulated  should  be  in¬ 
terpreted  with  care.  It  is  a  value  computed  on  the  basis  of  the 
five  distance  values  of  Table  X  with  weights  equal  to  the  number 
of  data  pairs.  Since  the  number  of  data  pairs  is  highly  variable 
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the  "distribution"  of  distances  is  far  from- "normal".  Also  no 
account  is  taken  of  the  variability  of  true  air  speed  within 
JETNO's,  shutter  time  variations,  etc. 

The  detailed  discussion  of  the  significance  of  the  data 
from  the  Project  Jet  Stream  flights  is  deferred  to  Section  F. 
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E.  EEIA-T UM  BETWEEM  MEAN  S  QUAEE_  CXMPQNENT  JJXEEEEEN  GES  AND  THE 

COMPONENTS  GF  IKE  WIND  COHRBOTIW  TE3JS0R  '  """" 

The  values  of  the  mean  square  velocity  component  differences 
are  related  to  the  components  of  the  wind  correlation  tensor 
through  the  following  elementary  equations.  A  mean  square  dif¬ 
ference  may  be  expressed  as 

A2  -  (u-u')e  (1) 

=  [(u-u)  -  (u'-u')  +  (U-U')]2 

=  |  u-u'  +  (U-U' ) J2 
=  u2  -  2uu '  +  (u')2  +  2u(U-U') 

(2) 

-2 u  '  ( U -U ' )  +  (U-TJ')2 

where  the  subscripts  are  suppressed  for  the  moment  and  where 

u  =  U  -  U,  u'  =  u'  -  U' 

If  the  points  P  and  p'  were  fixed  (where  U  and  U',  respec¬ 
tively,  are  measured),  then  u  =  ^  0  (these  are  the  means  of 

the  departure  from  the  mean)  and  if  the  mean  velocity  field  were 
uniform,  then  U  =  U'  so  that  (U-U')2  =  0. 

If  the  mean  velocity  field  is  nearly  uniform,  this  last 
term  may  be  considered  negligible.  If,  in  addition,  the  wind 
variability  at  P  and  P'  is  nearly  the  same,  u2  s;  (u')2,  one  may 
write 

A2  ==  2  (u2-uu  '  )  =  2ua(l-ruu)  (3) 

where  ruil  is  the  correlation  coefficient  relating  U  at  P  to  U' 
at  P'.  In  the  case  in  which  P  and  P'  are  fixed,  the  means  are 
based  on  a  period  of  observations  in  time. 
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In  the  case  of  the  4WW  and  Project  Jet  Stream  observations, 
the  points  P  and  F'  are  not  fixed  in  space  but  are  fixed'  with 
respect  to  each  other  and  the  means  are  taken  with  respect  to 
space  (i.e.,  means  over  a  flight  path).  This  raises  several 
questions  as  to  the  validity  and  interpretation  of  the  results. 

When  the  flight  path  is  straight,  as  in  the  case  of  the 
individual  Jet  Stream  Flights,  the  situation  is  nearly  the  same 
as  that  for  fixed  P  and  P'.  In  this  case  P  and  P-  are  fixed 
with  respect  to  the  orientation  of  P'  with  respect  to  P  (i.e., 
in  the  direction  of  the  heading).  The  distance  between  P  and  P' 
is  subject  to  some  variation  (when  fixed  time  lags  are  used)  due 
to  variations  in  the  Integrated  true  air  speed.  These  variations 
are  relatively  minor.  It  was  seen  that  for  the  Jet  Stream  data, 
the  value  of  the  standard  deviation  of  the  true  air  speed  was 
only  1.0  to  4.0$.  The  value  of  the  integrated  true  air  speed 
over  a  flight  path  of  some  length  would  be  less  the  longer  the 
flight  path.  Some  other  factors  also  Influence  the  air  spacing 
of  P  and  P'  such  as  the  uniformity  of  the  camera  shutter  trigger, 
the  uniformity  of  the  timing  clock,  etc.  It  is  assumed  that 
these  are  relatively  minor  factors. 

The  evaluation  of  the  correlation  coefficient,  ruu,  from 
the  values  of  A2  raises  the  question  of  determining  the  value 
of  u2,  the  basic  wind  component  variance.  In  the  case  of  flight 
data,  the  flight  is  made  at  a  specific  level  on  a  specific  day 
and  hence  in  a  very  particular  synoptic  situation.  The  value 
of  u2  from  such  data  would  be  a  reflection  only  of  the  line  of 
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flight  in  that  situation.  Consequently,  the  value  of  u2  would 
be  most  unreliable.  To  obtain  a  reliable  estimate  of  us  would 
require  many  flights  (more  or  less  identical)  so  tint  the  data 
on  u2  from  each  would  be  reasonably  representative.  (Another 
alternative,  flying  many  aircraft  over  the  world  simultaneously 
to  sample  the  entire  field  of  motion  seems  (very)  unfeasible 
and  not  free  from  objections.)  Consequently,  it  seems  unwise 
to  use  data  on  A2  to  actually  compute  correlation  coefficient 
values.  The  values  of  A2,  from  the  few  situations  at  hand  begin 
to  approach  a  "representative  sample"  (even  though  poorly)  so 
that  one  has  much  more  confidence  in  the  composite  results  of 
the  several  flights  than  in  those  for  a  single  flight. 

When  the  observations  are  taken  from  flight  paths  on  dif¬ 
ferent  headings,  the  question  of  the  structure  of  the  components 
of  the  wind  correlation  tensor  arises.  It  has  been  shown  that 
the  way  in  which  the  wind  component  correlation  coefficients 
vary  with  distance  is  a  function  of  the  orientation  of  P'  with 
respect  to  P  when  the  components  in  a  fixed  direction  are  used. 
If  we  were  to  use  for  U  and  U'  the  east  pointing  wind  component, 
then  to  adequately  study  the  variation  of  A2,  it  would  be  nec¬ 
essary  to  determine  Its  value  not  only  as  a  function  of  the 
time  lags  (distance)  but  also  as  a  function  of  flight  heading. 
This  increases  the  magnitude  of  the  problem  greatly.  On  the 
other  hand,  it  has  been  seen  that  the  wind  vector  is  resolved 
along  and  across  the  flight  path,  the  components  of  the  wind 
correlation  tensor  become  independent  if  the  orientation  of  P' 
with  respect  to  P  (in  certain  idealized  cases  which  seem  to  be 
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verified  for  synoptic  scale  disturbances),  Consequently,  this 
orientation  of  the  reference  frame  is  used.  At  least  in  this 
case  we  may  look  for  a  connection  between  the  values  of  A®  from 
synoptic  observations  and  values  from  aircraft  observations  where 
the  point  separation  distances  overlap  and  still  make  use  of 
what  aircraft  data  are  available. 

Then  the  flight  path  is  not  straight,  as  in  the  case  of  the 
4WW  observations,  the  only  interpretation  that  seems  feasible 
for  values  of  A2  Is  in  terms  of  the  longitudinal  and  transverse 
wind  components.  In  such  a  case  the  values  of  A;2  =  (U^-U')2 
and  At2  =  (U^-U')2  are  compiled  in  the  sense  of  an  ensemble 
average  in  which  the  values  over  the  flight  path  are  considered 
as  a  space  average.  In  this  case,  the  values  of  individual 
means,  their  differences,  and  standard  deviations  usually  asso¬ 
ciated  with  a  correlation  coefficient  computation  are  completely 
meaningless . 

The  interpretation  of  the  data  on  A^2  and  A^2  is  made  via 
equation  (3) .  This  requires  that  certain  assumptions  are  made 
concerning  a  meaningful  value  of  u2,  the  wind  component  variance. 
It  is  necessary  to  assume  that  such  a  value  exists  and  that  either 
it  is  Independent  of  the  orientation  of  P'  with  respect  to  P  or 
that  the  ensemble  averages  have  been  taken  with  the  orientation 
sufficiently  randomized  that  Its  value  is  a  suitable  mean  (and, 
of  course,  that  its  variation  with  respect  to  orientation  is  not 
too  large) . 

In  spite  of  the  lengthy  discussion  of  the  conceptual  diffi¬ 
culties  encountered  in  using  flight  observations  for  estimating 


VII-45 


the  behavior  of  A.s  and  A.  s  as  a  function  of  the  separation  of 

&  V  >•  . 

P  and  there  are  certain  outstanding  advantages  in  this  method 
of  attack. 

First,  the  values  of  (l-ruu),  ruu  =  correlation  coefficient, 
are  particularly  suitable  for  estimating  the  behavior  of  ruu  for 
the  separation  of  P'  from  P  near  zero.  The  use  of  the  co-ordi¬ 
nates  (log(d),  log(l-ruu))  or  equivalently  (log(d),  log  A2), 
where  d  is  separation  distance,  permit  a  "magnification"  of  the 
area  near  d  =  0,  ruu  =  1  that  increases  as  this  point  is  approached. 

A  second,  and  more  important,  advantage  lies  in  the  fact 
that  the  effects  of  at  least  certain  types  of  observation  error 
are  minimized.  It  was  pointed  out  in  the  first  section  of  this 
Chapter  that,  when  dealing  with  rawlnsonde  observations  it  is 
expected  that  for  very  closely  spaced  observations  one  is  re¬ 
duced  to  the  analysis  of  observation  error  and  that  for  small 
enough  separation  of  observation  points  correlation  behavior  is 
almost  entirely  determined  by  the  size  of  these  errors.  In 
dealing  with  velocity  component  differences,  one  can  by  no  means 
dismiss  observation  errors,  but  the  way  in  which  they  must  be 
accounted  for  and  the  range  of  their  effects  are  quite  different. 

Errors  that  are  relatively  constant  from  one  observation  of 
a  wind  component  to  the  next  cancel  out  when  the  component  is 
observed  with  the  same  "instrument"  and  component  differences 
are  used.  One  is  not,  then,  faced  with  evaluating  the  total 
instrument  error,  but  only  the  extent  to  which  the  error  differs 
from  one  observation  to  the  next.  This  has  a  tendency  to  reduce 
the  effect  of  errors  by  an  order  of  magnitude  or  so. 
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In  using  the  wind  observations  from  a  device  such  as  the 
dappler  navigator,  it  is  necessary  to  account  for  a  completely 
different  type  of  observational  error,  that  of  the  tendency  of 
the  electro-mechanical  devioe  to  "smooth"  the  observations. 

This  effect  is  most  prominent  when  the  observations  are  closely 
spaced  and  is  of  decreasing  Importance  for  increasing  spacing. 

The  use  of  quantity  (U^-U^) (U^-U^ )  to  estimate  the  struc¬ 
ture  of  the  correlation  tensor  components  u^u^  and  u'u^,  or  the 
correlation  coefficients  r^t  and  rtA,  is  impossible  from  the 
flight  data  at  hand.  This  is  due  to  the  fact  that  these  com¬ 
ponents  (correlation  coefficients)  are  quite  small  in  the  first 
place  and  In  the  second  place  are  direction  sensitive  (i.e., 
they  depend  on  the  orientation  of  P'  with  respect  to  P) .  To 
treat  these  components  would  require  that  the  flight  legs  be 
carefully  classed  in  terms  of  heading.  This  had  not  been  done, 
and  the  data  available  makes  It  Impractical. 
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F.  STRUCTURE  OF  WIND  CORRELATION  TBH3QR  OQMPQN1NT&  PROM  ’FLIGHT 

riflET"  '  ~  '  " 

1 .  Presentation  of  Combined  Flight  and  Rawlnsonde  Data 

The  composite  data  on  A^2  and  Afc®  from  both  the  4WW  and 
Project  Jet  Stream  flights  together  v/lth  the  values  from  JETNO-6 
and  computed  values  of  those  quantities  from  the  winter,  2-00  mb 
level,  SC-USWB-FCPA  data  are  plotted  in  Fig.  VII -9.  The  ordi¬ 
nate  is  the  logarithm  of  the  A^2  or  A^.2  (knots  squared)  while 
the  abscissa  is  the  logarithm  of  the  separation  of  point  P' 
from  P  (nautical  miles). 

The  first  striking  feature  of  this  illustration  lies  in 
the  large  scatter  of  the  points  of  the  3C-U3WB-FCDA  data.  This 
is  due  to  two  factors.  First,  port  of  the  scatter  is  due  to 
the  fact  that  the  values,  of  and  of  rtt  possess  some  scatter 
in  themselves.  The  second,  and  more  important  factor,  is  the 
fact  that  the  values  of  u,,2  and  u'2  (or  u^.2  and  u|2)  that  were 
used  were  those  actually  observed.  These  values  strongly  re¬ 
flect  the  fact  that  the  wind  variability  is  variable  from  place 
to  place  (especially  with  latitude).  Thus  the  values  of  A?2 
and  Afc®  which  had  been  plotted  show  a  much  larger  scatter  than 
those  values  of  1-r ^  and  l-r^. 

The  data  points  for  A^2  and  A^.2  from  the  two  flight  obser¬ 
vation  sources  demand  some  brief  comment.  The  most  striking 
property  that  these  points  display  is  their  tendency  to  lie 
parallel  to  each  other.  This  property  will  he  discussed  in 
more  detail  later.  That  they  do  not  lie  along  the  same  pair 
of  "parallel"  curves  (one  for  A^2  and  the  other  for  a"2')  is  due 
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Pig.  VII -9  Composite  Data  on  A®  and 
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at  least  -in  pant  to  the  effeot  of  the  term  in  u^2  (or  u^2)-  of 
equation  (3),  l,e.,  on  the  basic  wind  component  variance  con¬ 
cerned  . 

The  4WW  observations  are  made  in  middle  and  late  winter 
at  near  25 ,000  feet.  The  fact  that  these  are  winter  observa¬ 
tions  could  urge  a  large  value  for  u2.  The  fact  that  they  were 
at  25.000  feet  would  indicate  a  somewhat  smaller  value  of  u2 
than  would  be  expected  near  the  35,000-foot  level  where  the 
Project  Jet  Stream  flights  operated. 

The  composite  data  from  Project  Jet  Stream  reflects  very 
strongly  the  values  from  J'EThO-5,  a  long  flight  made  in  late 
June  and  reflecting  essentially  summer  (low)  values  of  u2. 

The  flight,  JETNO-0,  shown  individually  because  of  the 
close  spacing  of  P'  to  P,  was  made  in  late  May  and  so  tends  to 
reflect  transition  conditions  from  winter  to  summer.  On  the 
other  hand,  this  flight  w as  made  at  41,000  feet,  very  near  the 
level  or  maximum  u2 . 

To  summarize,  the  vert leal  displacement  of  the  flight  data 
points  is  strongly  dependent  on  the  value  of  u2  that  might  be 
considered  as  prevailing.  This  value  is  affected  by  such  a 
variety  of  factors,  ouch  as  season,  altitude,  heterogeneity  of 
data,  etc.,  that  it  seems  somewhat  unreasonable  to  evaluate 
the  amount  of  the  relative  vertical  displacement  explicitly. 

2 .  Correlation  Structure  at  Shorter  Distances 

The  form  that  the  curves  of  Fig.  VII -9  should  take  is 
subject  to  some  conjecture.  The  empirical  form  of  the  correla¬ 
tion  coefficient  that  was  used  for  the  synoptic  scale  data  in 
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previous  chapters  has  been,  in  simplified  form,  such  as  to  give 

If  -  2^(1  -  e-rS/SL£) 

and  (4) 

=  2u1"[l  -  (1  -  re/LS)e"rS/2L2] 

where  u2  is  of  no  significance  in  detail  but  represents  a  mean 
square  component  standard  deviation.  The  curves  for  this  pair 
are  shown  by  the  dashed  lines  in  Pig.  VII-9,  labeled  k  -  1.0. 
Obviously,  they  drop  much  too  rapidly  with  decreasing  distance. 

As  shown  previously,  if  these  correlations  are  modified  by  in¬ 
troducing  an  "error"  of  observation,  then  they  have  a  horizontal 
asymptote  for  small  distance  values.  These  are  the  dashed  curves 
labeled  k  =  0.9  and  k  =  0.8  where  k  is  a  measure  of  the  ratio  of 
observation  error  to  natural  variability.  Neither  situation  Is 
tolerable . 

The  situation  may  be  described  in  somewhat  different  terms 
by  assuming  that  the  wind  components  are  made  up  of  terms  addi¬ 
tive  ly  combined  and  Independent  of  each  other.  Thus  let 

u  =  u:  +  u8  + - +  u# 

so  that 

uu '  =  u:  +  ua  u8'  + - +  un  un' . 

The  cross  product  terras  are  all  zero  from  the  independence  assump 
tlon .  Then 

r  =  kjrj,  +  k8r8  +  —  +  k,  rB 
kx  +  ka  h -  +  k„  «  1. 

where  the  k's  are  constant  factors  which  depend  on  the  relative 
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contribution  of  each  term  to  the  total  one-point  wind  variability 
and  the  r's  are  the  correlation  coefficients  for  each  term  sepa¬ 
rately. 

For  example,  suppose  that  the  term  in  uz  is  due  to  the 
"synoptic  scale"  wind  variability  as  above  and  that  the  u2  term 
is  due  to  smaller  eddies  with  a  correlation  structure  given  by 

e-r/D 


where  D  is  a  distance  parameter  that  measures  the  size  of  the 
smaller  eddies.  Then  the  previous  expressions  become 

17 . 2^(1  -  *,.-'***  -  kae-r/D) 


and- 


(5) 


At2  =  2u2  (  1  -  !o  (1  -  i»a/L8)e“rS/^L*J  -  k8  ( 1  -  r/D)e“r//D) 


The  fit  of  this  empirical  form  to  the  data  points  leaves 
something  to  be  desired,  but  is  setter  than  the  former  alterna¬ 
tives.  The  main  difficulty  seems  to  lie  that  the  "shift"  from 
the  form  of  (4)  to  the  form  that  would  predominate  at  very  small 
separations  (where  the  slope  of  the  curve  is  near  unity)  appears 
prominently  in  the  empirical  curves  but  does  not  appear  In  the 
data  points. 


*  The  form  of  the  correlations  in  A^2  is  associated  with  those 
in  the  expression  for  A  2  through  the  continuity  equation. 

X/ 

See  Chapter  V,  Section  B. 
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G.  "MSP  SCALE"  EDDY  STRUCTURE  ALOFT 

The  data  on  A^a  and  shown  in  Pig.  VII -9  clearly  indi¬ 
cate  that  the  wind  variability  at  two  points  is  due  to  at  least 
two  factors,  first  a  systematic  eddy  structure  on  the  "synoptic" 
scale,  i.e.,  to  the  traveling  disturbances  seen  on  the  synoptic 
charts  with  an  average  size  of  near  600  n  mi  (as  defined  by  the 
statistical  wind  parameter  L)  and,  second,  a  "tapering  off"  in¬ 
to  eddies  of  smaller  scale  and  intensity.  The  Jet  Stream  data^ 
clearly  indicate  the  presence  of  these  smaller  scale  eddies  which 
are  loosely  called  "meso  scale"  eddies  here.  Such  eddies  are  of 
a  size  that  would  be  undetected  in  the  present  network  of  rawin- 
sonde  observations  except  for  their  effect  on  a  single  station. 
The  presence  of  such  eddies  is  evident  on  the  graphs  of  the 
Project  Jet  Stream  data  which  appear  in  the  data  volumes. 

Exerpts  from  these  graphs  are  shown  in  Figs.  Vll-lOa  and  Vll-lOb. 

In  these  figures,  the  top  curve  of  data  points  represents 
the  magnitude  of  the  parallel  wind  component,  the  middle  curve 
the  perpendicular  component,  and  the  bottom  curve  the  tempera¬ 
ture  (of  no  immediate  concern  here).  The  terms  parallel  and 
perpendicular  are  defined  with  respect  to  the  direction  of  the 
maximum  wind  speed  encountered  during  the  flight.  It  is  evident 
from  the  data  points  that  there  are  small  perturbations  of  the 
wind  velocity  components  with  a  magnitude  of  a  few  knots  and  that 
these  perturbations  are  not  entirely  random. 

To  obtain  a  closer  look  at  the  meso  scale  eddy  structure 
shown  in  these  figures,  the  "synoptic  scale"  wind  was  subtracted 
from  the  total  wind.  This  was  relatively  easy  to  approximate  in 


Pig.  VTI-lOa.  Wind  Components  and  Temperature;  JETNO  4  (Jet  E22 
Pressure  Altitude  34,360  ft.) 
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these  cases.  Each  flight  is  of  rather-  short  duration  so  that 
it  covers  only  a  small  part  of  a  "disturbance "  on  the  synoptic 
scale.  Consequently,  the  synoptic  scale  wind  for  each  compon¬ 
ent  was  approached  by  fit  ting  the  data  for  each  wind  component 
by  a  quadratic  1'unc'l.on  of  rim---.  It;-  residual  wind  c  0.7170  on  ants 
obtained  when  t  he  !  1  -cl  wind  component  values  war-  sub- 
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Fig.  VII-11  Eddy  Velocity  Components,  JETNO- 
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Fig.  VI I -12  Eddy  Velocity  Vectors ,  JETNO-l 


VII -5  8 


These  serial  numbers  and  the  corresponding  time  in  hours, 
minutes,  and  seconds  is  shown  as  the  abscissa  in  Pig.  VII-11. 
Several  sections  are  shown  displaced  vertically,  simply  to 
permit  a  scale  large  enough  to  show  the  eddies  clearly  but  to 
keep  the  whole  figure  to  reasonable  dimensions. 

The  data  values  of  Figs.  VII-11  and  VII-12  are  carried 
somewhat  beyond  the  real  limits  of  accuracy  of  the  tabulated 
data  values.  Wind  component  values  were  tabulated  to  the  near¬ 
est  knot  while  the  values  plotted  in  Pig.  VII-11  are  to  the 
nearest  tenth  of  a  knot  and  those  of  Pig.  VII-12  to  the  nearest 
half  knot.  The  irregularities  of  the  points  in  Fig.  VII-12 
clearly  indicate  this.  It  is  also  evident  that  the  eddy  veloc¬ 
ities  are  of  an  order  of  magnitude  greater  than  these  irregu¬ 
larities  . 

The  size  of  the  eddies  shown  in  Pig.  VII-12  appears  to  be 
reasonably  regular.  By  averaging  the  spacing  between  reversals 
of  the  eddy  velocity  with  respect  to  the  flight  path,  it  amounts 
to  approximately  13  time  intervals  or  about  40  n  mi.  The  magni¬ 
tude  of  the  eddy  components  appears  to  be  3  to  4  knots. 

The  eddy  velocity  correlation  coefficients  for  the  compon¬ 
ents  of  the  velocity  correlation  tensor  were  computed  for  JETNO-1 
through  JETNO-5  and  are  illustrated  in  Fig.  VII-13.  The  data 
itself  In  tabular  form  is  contained  in  Appendix  B-IV  together 
with  the  means  and  standard  deviations  of  the  velocity  compon¬ 
ents.  (These  correlation  coefficients  are  not  the  same  as  those 
appearing  in  A^2  and  A^2  discussed  previously.  The  correlations 
that  appear  :ln  such  expressions  involve  the  total  wind  as  a 
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function  of  distance  while  the  correlations  here  discussed  in¬ 
volve  the  meso  scale  eddy  part  of  the  total  wind.) 

There  are  several  features  in  common  in  the  five  parts  of 
Fig.  VII-13.  Some  of  these  are: 

i)  Both  r^  and  rfct  decay  from  the  value  of  unity  at 
zero  separation  in  an  exponential  manner. 

li)  The  value  of  r  ^  (as  well  as  that  of  r^fc)  decreases 
to  negative  values. 

iii)  Near  zero  separation,  either  r^  or  r^  may  be  the 

larger. 


iv)  The  values  of  and  r^  seem  to  oscillate,  or 
tend  toward  such  behavior  for  larger  separations. 

v)  The  values  of  v ^  and  rfc  near  zero  distance  may 
be  appreciably  different  from  zero. 

vi)  The  values  of  r^  and  rfc  also  tend  to  oscillate 
at  larger  distances. 

Since  the  values  of  r^  and  r^.  are  significantly  different 
from  zero,  especially  at  zero  distance,  and  since  the  values  of 
the  standard  deviations  of  the  longitudinal  and  transverse  com¬ 
ponents  were  also  appreciably  different,  the  one-point  frequency 
distribution  of  the  eddy  wind  must  be  far  from  circularly  normal. 
To  test  whether  there  was  a  preferential  co-ordinate  system 
oriented  in  the  direction  of  the  maximum  eddy  wind  variability, 
the  reference  axes  were  rotated  into  this  direction  and  the 
corresponding  correlation  coefficients  computed.  (The  rotation 
was  carried  out  in  terms  of  the  covariance  forms  so  that  it  was 
not  necessary  to  re-compute  from  the  original  data.)  The  results 
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are  shown  in  Fig,  VII --14  and  the  values  are  tabulated  in  .Appendix 
B-I7  together  with  the  accessory  quantities. 

Since  we  have  used  terms  and  the  subscripts  corresponding 
to  longitudinal  and  transverse  to  refer  a  co-ordinate  system 
oriented  along  and  at  right  angles  to  the  two  points  concerned, 
the  terminology  associated  with  the  components  in  Fig.  VII-14 
is  changed.  The  v-rins  caryll-I  and  normal  to  the  direction 
concerned  will  be  us-.d  an i  indicated  by  corresponding  first 
letter  subscript.'..  Tm  nix  it  ms  characterizing  the  correlation 
coefficients  c.r  Fig.  VII-.. v  ar  :  still  valid  for  Fig.  VII-14  with 
the  exception  of  Item  v) .  -.a  pLaoe  of  Item  v)  above,  the  dis- 

-ussion  of  Fig.  VII-l’i  would  not-;  that  in  those-  cases  in  Fig. 
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be  discussed  separately,  vn  ■  data  in  tills  case  comes  from  a 
particularly  long  flight ,  1  !.  •  orrc-lations  are  based  on  over 
'■'.00  data  pairs,  and  ■  .  •  p-nl  ly  the  correlations  repre-s-r.t  a 

sample  from  a  wide  varie*  y  of  conditions  (summer).  The  data 
entering  the  '  "’me la’  inns  of  .''EFMJ-1  through  .3ETN0--:-i  are-  from 
short  flight  legs  and  r-pr-T.sent  small  samples  of  only  a  f>w 
eddies  each  (il  identifiable  eddies  in  the  case  of  VETMO-l, 
about  half  this  number  if  a  full  "wave  length"  is  considered 
to  constitute  an  eddy). 

The  characteristics  noted  regarding  the  behavior  of  th-. 
correlation  coefficients  corresponding  to  the  components  of  the 
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correlation  tensor  for  eddy  velocities  Imply  several  significant 
items  regarding  the  kinematic  structure  of  the  eddy  motions. 

First,  the  fact  that  r  ^  and  r^  were  distinctly  non-zero 
at  zero  distance  and  since  and  are  quite  different,  one 
would  conclude  that  the  eddy  structures  are  strongly  asymetricalo 
This  is  especially  true  since  there  was  a  greater  difference  be¬ 
tween  c?  and  cn  . 

Second,  since  r^  and  r^^.  (rpp  and  rn  n  )  (and  also  the  cross 
;•  orrelations )  tend  to  oscillate,  it  would  seem  to  indicate  that 
eddies  of  alternate  sense  of  rotation  form  a  reasonably  regular- 
pat  tern  (over  relatively  short  distances). 

Tnird,  the  more  rapid  decay  of  the  value  of  r,,  as  compared 

tt 

with  in  discussing  the  synoptic  scale  eddies  in  the  atmosphere 
Is  strictly  the  case  for  eddies  with  a  stream  function.  If  the 
eddy  velocity  v/ero  determined  by  a  potential  function  of  the  same 
geometrical  pattern,  the  roles  of  the  velocity  components  would 
be  interchanged.  In  such  a  case,  the  longitudinal  component  would 
decay  more  rapidly  than  the  transverse.  The  correlations  observed 
are  some  of  one  and  some  of  the  other  type.  From  this  it  would 
seem  reasonable  to  assume  that  the  meso  scale  eddies  are  asso¬ 


ciated  with  local  two  dimensional  divergence  fields  that  ara 
relatively  much  stronger  than  those  associated  with  synoptic 
scale  disturbances. 


Fourth,  if  one  uses  the  distance  at  which  r^  or  (or 
rpp  or  rnn)  first  reaches  the  zero  value  as  a  scale  parameter, 


this  scale  appears  to  lie  between  20  and  40  n  mi. 
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The  case  of  JETNO-5  Is  somexvhat  different  from  the  first 
four,  particularly  in  that  the  rate  of  decrease  of  correlation 
with  distance  is  smaller.  This  would  indicate  a  larger  eddy 
size  in  general.  This  flight  was  made  in  summer  and  extended 
across  the  continent  so  that  many  more  eddies  were  "sampled" 
under  conditions  quite  different  from  the  flights  discussed 
above . 

The  effect  of  the  time  constant  of  the  doppler  navigator 
system,  without  doubt,  affects  the  correlations  at  the  smaller 
distances.  In  all  cases  (except  JETNO-3  in  parallel  and  normal 
components)  the  correlations  of  like  components  "pinch  together" 
as  the  correlations  approach  1.  It  is  understood  that  the  time 
constant  of  the  doppler  system  is  normally  about  30  to  45  seconds 
but  that  for  the  flights  concerned  the  value  may  have  been  espe¬ 
cially  reduced  to  near  15  seconds.  Since  the  data  readings  were 
made  about  30  seconds  apart,  the  effect  of  the  time  constant  of 
the  system  is  without  doubt  present  if  the  former  is  the  case 
and  is  probably  present  If  the  latter  were  true. 
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CHAPTER  VIII 


APPLICATIONS 

This  chapter  is  devoted  to  a  discussion  of  some  applications 
of  the  ideas  developed  in  the  preceding  chapters.  It  is  not  in¬ 
tended  to  be  complete  or  exhaustive.  The  reader  will,  without 
doubt,  have  problems  at  hand  to  which  application  of  the  basic 
philosophy  may  be  made,  if  not  the  specific  results. 

There  are  two  basic  ideas  behind  the  applications  which 
should  be  made  explicit.  These  are: 

1.  In  spite  of  the  discussions  of  variation  of  parameters 
with  season,  level,  etc.,  the  correlation  coefficients  that 
appear  are  remarkably  conservative.  They  have  very  nearly  the 
same  parameter  values  over  the  entire  hemisphere  (outside  the 
tropics).  Consequently,  instead  of  speaking  of  correlations 
between  properties  at,  say,  Washington  and  Buffalo,  it  is  pos¬ 
sible  to  speak  of  correlations  of  these  properties  between  any 
two  points  separated  by  a  distance  of  2 56  n  mi  and  with  an  azi¬ 
muth  of  the  line  joining  them  of  343°.  Further,  if  reference 
directions  are  suitably  chosen  it  may  be  possible  to  ignore  the 
azimuth. 

2.  The  correlation  coefficients  concerned  deal  with  the 
properties  of  a  continuum.  Consequently,  these  correlation  co¬ 
efficients  must  have  certain  definite  properties  as  functions  of 
space  and  time  separations.  They  are,  for  example,  non -negative 
Hermitian  functions;  they  are  continuous  except  for  a  possible 
jump  discontinuity  at  zero;  etc. 
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This,  of  course,  comes  from  the  fact  that  the  structure  of 
the  correlation  coefficients  reflect  the  basic  dynamical  char¬ 
acteristics  of  the  atmosphere  and  insofar  as  these  characteristics 
are  unchanged  in  space  and  time,  so  is  the  structure  of  the  cor¬ 
relation  coefficients.  In  fact,  to  carry  the  case  further,  the 
correlation  coefficients  in  normalized  form  do  not  reflect  the 
background  parameter  variability  which  appears  in  the  standard 
deviations  (one-point  variances)  so  that  the  space  and  seasonal 
variations  of  these  basic  variabilities  may  be  handled  separately. 
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A.  EXTRAPOLATION  FROM  LIMITED  DATA 

The  term  extrapolation  is  used  here  in  the  sense  of  both 
extrapolation  in  space  and  in  time.  The  term  "forecasting" 
might  well  be  used  when  extrapolation  in  time  (or  time  and 
space)  is  involved.  The  forecasting  techniques  are,  of  course, 
the  most  elementary  of  the  statistical  type.  For  this  reason, 
the  term  "forecasting"  was  not  used. 

1 .  Single  Station 

(a)  The  Basic  Relations 

In  extrapolating  information  from  a  single  station, 
an  extension  of  the  ideas  presented  in  the  preceding  chapters  is 
needed.  This  is  the  form  of  the  correlation  between  wind  and 
height  on  an  isobaric  surface.  The  empirical  form  of  this  cor¬ 
relation  is 

“  0 

r  -  _  %-rV2L* 
nt  "  Tr 

where  r^  =  the  correlation  coefficient  relating  the  height  at 
P  to  the  longitudinal  wind  component  at  P'  and  r^t  =  the  cor¬ 
relation  coefficient  relating  height  as  P  to  the  transverse  wind 
component  at  P' . 

That  this  should  be  the  form  of  these  expressions,  consider 
the  geostrophic  wind  equations  in  the  form 

K  -  -  & 

K  -  sir 


where  the  x-direction  is  along  PP'  and  the  y-direction  is  normal 
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PP' .  Then 


On  the  right  sides,  we  have  (following  Appendex  A-l) 


where  the  correlations  arc  between  the  longitudinal  and  trans¬ 
verse  wind  components  at  P  and  the  height  at  P'. 

The  preceding  relations  may  be  used  to  extrapolate  either 
wind  components  or  height  from  a  single  observation.  Thus,  if 
one  uses  a  simple  regression  form  (standardized  variables) 

h'  =  at  ut  +  aaut  +  a3h, 
u'  =  bj^Uj^  +  baut  +  k3h, 
v'  ~  ciui  +  caut  +  c3h, 
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then  it  follows  at  once*  that 

h '  -  £(r/L)exp(-r2/2L2)^](ut )  +  j^exp  ( -ra/2L2)^j(h) , 

u'  =  [exp(-r2/2L2)](uA),  ^ 

u l  =  [(1  -  ra/La)exp(-ra/2La)](ut) 

+  [(-r/L)exp(-ra/2La)](h) . 

(t>)  Discussion 

i)  Extension  of  Geostrophic  Relation 

The  preceding  extrapolation  relations  are 
interesting  in  two  ways.  First,  the  correlation  coefficients 
arid  rth  may  be  thought  of  as  a  statistical  extension  of 
the  geostrophic  wind  equation.  It  is  well  known  that  the  nu¬ 
merical  value  of  the  height  of  the  isobaric  surface  does  not 
provide  information  on  wind  and  conversely.  The  geostrophic 
wind  equation  requires  height  gradient.  On  the  other  hand, 
every  meteorologist  knows  that  if  the  height  of  the  isobaric 
surface  is  known,  especially  if  it  departs  appreciably  from 
normal,  this  information  alone  makes  it  possible  to  construct 
a  rough  picture  of  conditions  throughout  the  surrounding  area. 
This  is  possible  through  the  correlations  rtS  and  r* t  which 
are  zero  at  zero  separation  but  which  attain  a  maximum  of  0.607 
at  r  =  L.  Thus,  in  estimating  height,  knowledge  of  the  trans¬ 
verse  wind  contributes  appreciably  to  the  estimate  and  in 


*Use  is  made  of  the  correlations  rhh,  r  ,  and  rtt  which  are 
obtained  in  the  previous  chapters.  11 
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estimating  wind,  knowledge  of  the  height  contributes  much  (in 
this  case  the  most,  just  where  the  information  on  wind,  r^t, 
becomes  zero) . 

ii )  Reasonable  and  Silly  Estimates 

In  the  second  place,  it  is  a  common  miscon¬ 
ception  that  a  statistical  forecast  procedure  only  gives  esti¬ 
mates  that  tend  toward  the  mean  and  cannot  give  estimates 
indicating  further  departure  from  the  mean.  The  first  and 
last  of  the  estimates  above  are  examples  to  the  contrary. 

Just  how  this  comes  about  is  readily  seen  in  the  case  of  the 
estimate  for  height,  H'.  This  obviously  arises  from  the  first 
terms  of  a  Taylor  expansion  of  the  height  of  the  isobaric 
surface : 


H'  =  II  -i-  |§Ax  +  — 

OA 

*sTJ  \ 

where  .  If  we  divorce  ourselves  of  physical  ideas  and 

put  the  emphasis  on  the  observed  quantities,  one  obtains  the 
linear  form 


H'  -  aaH  +  aaUt 

where  the  coefficients  a:  and  a8  are  to  be  determined  from  the 
usual  statistical  methods.  The  physical  method  amounts  to 
extrapolation  along  the  tangent  to  the  isobaric  surface  in  the 
direction  concerned.  This  type  of  extrapolation  is  strongly 
limited  in  the  range  of  distance  used  (Ax  =  r  very  small) .  If 
the  distance  is  excessive,  the  resulting  estimates  are  ridicu¬ 
lous.  On  the  other  hand,  the  statistical  estimate  based  on  the 
same  relations  also  starts,  for  small  Ax  =  r,  estimating  along 


•  •  •:  4  - ilr'-Itt'SHf  ait  M-iBI-  ■■ 
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the  tangent  but  departs  therefrom  and  eventually  estimates  the 
mean  value  (zero  in  the  case  of  standard  variables).  Thus,  it 
may  be  said  that  the  statistical  estimate  behaves  like  the 
physical  estimate  if  not  pushed  too  far  and  that  in  any  case 
is  prevented  from  giving  silly  answers  (see  Pig.  VIII -1 ) . 

The  estimate  of  height,  H',  also  provides  an  estimate  of 
the  location  of  the  nearest  High  or  Low.  This  is  obtained  by 
simply  differentiating  the  expression  for  H'  with  respect  to 
distance,  setting  the  result  equal  to  zero,  and  solving  for  r. 
This  involves  only  a  simple  quadratic  in  distance  and  may  be 
solved  graphically. 

iii)  Where  the  Variables  are  Located 

The  origin  of  the  estimate  for  the  transverse 
wind  component  is  not  Immediately  obvious.  Were  it  to  have 
come  from  a  Taylor  expansion,  the  first  term  in  Ut  would  be 
present  but  the  next  term  would  have  been  3U^/5x,  a  quantity 
not  readily  observed  at  a  point.  In  place  of  9Ut/ax  there  is 
the  totally  foreign  term  in  height,  H.  To  trace  the  origin  of 
this  term  it  is  necessary  to  go  back  to  the  Taylor  expansion 
of  height  in  the  "statistical"  form 

H'  =  axH  +  a8Ut  . 

Note  here  that  ax  and  aa  are  functions  of  the  distance  between 
P  (where  H  and  are  measured)  and  P'  (where  H'  is  estimated). 
Even  in  the  standard  Taylor  expansion,  where 

aj,  s  1, 

aa  =  (Vg)(Ax)  =  (Vs)r, 
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Fig,  VIII-1  Physical  vs.  Statistical  Extrapolation 
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these  coefficients  are  the  functions  of  the  variables .  When 
the  expression  is  differentiated  with  respect  to  distance  to 
obtain  for  the  transverse  wind  component. 


\j'  MiUJ  4  ,$.a2.TJ 

g°t  ~  aT1  ax^t 

and  so,  if  new  coefficients  are  introduced  to  simplify  the 
relation,  then 

=5  bxH  +  b8Ut 


and  the  same  "observables"  are  retained.  In  the  physical  case, 
Ba^Sx  =  0,  which  is  not  the  case  for  a  statistical  extrapola¬ 
tion.  If  the  extrapolation  for  height  is  differentiated  with 
respect  to  distance,  the  extrapolation  for  transverse  wind  is 
obtained  at  once,  except  for  a  factor  l/L  ’which  would  be  elimi¬ 
nated  in  going  to  the  standardized  variables  since  goh/Aau  a  L. 
(c )  Refinements  and  Extensions 

The  simplified  estimates  of  equation  (l)  using 
the  empirical  form  of  the  correlation  coefficients  may  be 
sharpened  somewhat  by  including  the  terms  in  the  gradient  of 
the  standard  deviation  of  height.  These  terms  were  included 
in  the  derivation  of  the  covariances  from  which  the  correlation 
coefficients  came,  but  were  dropped  as  being  negligible.  When 
this  is  done,  the  terms  are  no  longer  dependent  on  distance 
only,  but  also  depend  on  the  orientation  of  the  line  joining 
the  point  of  observation,  P,  and  the  point  of  estimate,  P' . 

The  time  factor  may  also  be  included,  but  in  this  case  it 
is  necessary  to  go  to  a  standard  geographical  co-ordinate 
reference  system.  Thus 

h'  =  aau  +  aav  +  a3h,  etc. 
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The  idealized  correlations  have  been  employed  in  the  sense 
that  when  like  properties  are  correlated,  the  coefficients 
approach  unity  as  separation  approaches  zero.  It  would  be  more 
realistic  to  degrade  these  coefficients  by  multiplying  by  a 
factor  slightly  less  than  one  to  account  for  observation  errors 
and  small  eddy  effects. 

2 .  Some  Geometrical  Problems 

(a)  A  Weather  Reconnaissance  Problem 

Suppose  that  there  are  observations  of  H  at  Px , 

- ,  Pn  ,  equally  spaced,  all  lying  on  a  straight  line.  These 

data  are  to  be  used  to  estimate  H'  at  P'  by  a  linear  regression 
(standardized  variables) 

h'  -  ax  h,  i - +  an  iin  . 

The  normal  equations  are 


ri  y 

-  ax 

aa  ri  s 

+  —  -i-  aln 

ra  y 

=  airai 

a3 

a3n 

r3n 

y 

=  alr»l 

1  a3 

rn  a  +  —  + 

an 

where  ryl  Is  the  correlation  coefficient  between  h '  and  ht  and 
rt i  is  that  between  ht  and  hj  (rtJ  =  r }  s ) .  The  coefficients 
are  given  by 
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where  M  is  the  determinate  of  coefficients 

M  = 

Now  let  rtJ  =  exp^ -( j -i )  2d2/2L2^J  where  d  is  the  distance  between 
successive  points  along  the  line  and  let  P'  be  a  distance  D  from 
the  line  and  opposite  point  ?3  (see  Fig.  VIII-2 ) ,  then 

r3y  =  exp{  -  [d2  +  (j-i)2d2]/2L2},  rt  3  =  exp[  -  D2/2L2]. 

Thus  if  j  =  i,  then  the  column  of  elements  r3  y  is  the  same  as 
rt j  in  M  except  for  the  factor  exp  ( -De/2  L2 )  and  a3  -  exp( -D2/2L2) . 
If  j  /  i,  then  the  column  of  elements  r3  y  is  a  multiple 
exp(-D2/2L2)  of  another  column  and  so  a3  =  0. 

Thus,  the  whole  estimate  of  H'  at  P  is  due  to  the  observa¬ 
tion  at  P3  and  all  other  observations  at  the  other  points  con¬ 
tribute  nothing  to  the  estimate  of  H'.  It  would  appear  that,  as 
far  as  linear  extrapolation  of  height  of  an  isobaric  surface  is 
concerned,  the  usual  system  of  weather  reconnaissance  may  leave 
something  to  be  desired. 

(b )  An  Example  from  Time  Series 

The  example  above  has  its  intrinsic  properties 
based  on  the  functional  form  of  the  correlation  function. 

Another  well-known  example  from  time  series  analysis  is  the 
following.  Consider  the  extrapolation  of  the  value  of  H  from 
a  sequence  of  values 


^  3  >  >  rln 

r3 1  ,  1  »  >  ^an 

rni  vm  >  »  1 


h '  =  +  aah8 


<|f!i 
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(standard  variables).  Then  it  is  well  known  that  the  mean 
square  error  of  estimate  is  given  by 


=  1 


ai  ryl 


-  3-a  rya 


1  -  r 


xa. 


-  r 


ix. 


-  r 


2  y 


+  2r 


Xy.Aaji.1 2. 


■  1  a 


Suppose,  now,  that 


■  1  a 


=  e-*T, 


ri„  =  e 


a  y 


=  e 


"0t 

} 

-J3(t+T) 


where  t  =  time  of  estimate  following  of  the  observation  hx  and 
T  =  time  elapsed  between  the  observation  ha  and  hx .  Then 

1  _  e-2^  -  e~2^  -  e~2^t+'r)  +  2e"^t+t+T+T) 


s2  = 


.  P~2^ 


-  1  -  e“2^ 

so  that  the  whole  estimate  of  h'  at  t  units  after  the  observa¬ 
tion  hx  is  dependent  on  hx  and  the  observation  h3  at  T  units 
prior  to  hx  contributes  no  information  to  the  estimate. 

The  above  example  has  applications  of  great  Importance  in 
the  time  sequence  extrapolation  of  meteorological  parameters. 

It  has  been  proposed  that  the  form  of  the  time  lag  correlation 
of  wind  follows  this  functional  form'  This  form  Implies 

that  for  extrapolation  purposes  the  latest  observation  contrib¬ 
utes  all  of  the  available  information.  There  is  no  doubt  that 
the  latest  observation  contributes  by  far  the  most  information. 
From  simply  the  philosophical  point  of  view,  i.e.  to  not  close 
the  door  too  soon,  it  would  seem  wise  to  require  the  greatest 
amount  of  proof  before  accepting  such  a  correlation  function. 
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The  choice  of  the  empirical  function  to  represent  the  cor¬ 
relations  is  not  a  simple  one.  In  Chapter  V  the  problem  was 
discussed  in  connection  with  the  longitudinal  and  transverse 
correlation  functions.  The  choice  there  was  made  on  the  basis 
of  three  properties;  it  must  be  reasonably  simple,  it  must  ex¬ 
plain  the  correlations,  and  the  functions  for  the  longitudinal 
and  transverse  correlations  must  be  related  in  a  specific  way. 
Several  examples  were  given  which  adequately  fulfilled  the  first 
and  last  conditions,  but  where  they  (all  but  one)  failed  was  to 
estimate  the  minimum  value  of  the  transverse  correlation  co¬ 
efficient  which  occurs  at  a  considerable  distance.  The  appear¬ 
ance  of  the  correlation  values  for  the  longitudinal  correlation 
coefficient  when  plotted  as  a  function  of  distance  strongly 
suggest  the  simple  exponential  decay  law.  In  fact  it  is  doubt¬ 
ful  whether  ordinary  test  for  goodness  of  fit  could  distinguish 
between  such  a  decay  and  the  squared  exponential  decay  selected. 
For  this  reason,  others  have  used  quite  inconcistent  empirical 
expressions  in  their  attempts  to  fit  the  data  to  empirical 
formulas  (Chapter  V,  references  2-4 J. 

( c )  Silent  Area  Problem 

The  two  examples  discussed  above  both  lead  to  the 
same  conclusion,  that  the  "nearest"  data  point  contributed  all 
of  the  information  to  the  estimate.  The  two  problems,  though 
they  have  the  same  answers,  are  basically  quite  different. 

What  is  going  on  will  be  explored  in  two  dimensions  in  what 
might  be  called  the  Silent  Area  Problem. 
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Suppose  that  it  is  required  to  estimate  H'  at  P',  that  an 
observation  of  Hx  at  Pj_  is  available,  that  no  observations  may 
be  made  between  Px  and  P'  (to  eliminate  the  trivial  case),  and 
that  one  other  observation  of  Ha  at  Pa  may  be  made.  Where  is 
the  best  location  for  P8  in  the  sense  that  the  estimate  of  H' 
by  a  linear  regression  on  Hx  and  Ha  has  the  least  error?  (See 
Pig.  VTII-3.)  In  standard  variables,  the  regression  is 

h '  =  axhx  +  aaha 

and  the  error  of  estimate  (standardized)  is 


Let  PXP'  =  D,  PxPa  =  d,  PaP'  =  L,  angle  P'PxPa  =  e. 


< 

1 
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Fig.  VIII-3 


Geometrical  Relations  in  the  Silent  Area  Problem 
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(note  that  the  L  here  is  not  the  size  parameter  (l/a)  for  which 
this  letter  was  formerly  used).  Then 


a  -a£D£  e-a2(De+13a)  (x  .  e-»8M' 


cos8  2^ 


(3) 


1  -  e 


a,  =  e 


^_ffiaD2/2  !  _  g~a2(d2  -  Dd  »  cos6 ) 

-a2d2 


1  -  e 


-o£  (D£-hJ£ ) /2  -  1 

'  '  j  _  e-a6ds 

Along  the  line  9  =  +  tt/2,  a8  =  0.  This  is  the  line  through  P4 
perpendicular  to  Pj^P7.  Along  this  line  the  information  at  P3 
contributes  nothing  to  the  estimate  at  P'  (see  section  (a) 
above).  On  the  circle  D*cos9  =  d,  the  circle  with  PXP'  as 
diameter,  Pa  is  the  closer  to  P'  and  consequently  all  of  the 
estimate  comes  from  Pa .  This,  however,  is  the  trivial  solution. 

Consider  Pa  lying  on  a  circle  of  radius  d  about  the  point 
Px  .  The  factor  (l  -  ea  Dd’c0°9)  in  the  expression  for  s2  takes 
on  its  largest  value  in  the  range  +  tt/2  <  0  <  3ir/2  (to  exclude 
the  trivial  solution)  when  9  -  ir.  Thus  ss  is  least  and  the 
smallest  error  of  estimate  occurs  when  P3  lies  on  the  line  P'Pj 
extended  (beyond  Px ) .  (Of  course,  the  error  of  estimate  is  less 
if  Ps  lies  between  P:  and  P',  but  this  is  the  trivial  case.) 
ii)  Simple  Exponential  Decay 

Consider  now  the  case  in  which 


r 


i  a 
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a  simple  exponential  decay  as  a  function  of  distance.  Then 

:,2aD  (e-^-e-^))° 

1  -  e-^d 

It  is  at  once  apparent  that  the  last  term  is  zero  when  L  =  D  +  d 
i.e.  when  P3  lies  on  PjP'  extended.  This  is  the  line  along 
which  the  observation  at  Pa  contributes  nothing  to  the  estimate 
at  P'  and  the  case  corresponds  to  that  in  section  (b).  Let  P3 
lie  on  a  circle  of  radius  d  centered  at  P:  .  Then  from 

L2  =  D2  +  d2  -  2Dd‘COs9  -  (D+d)2  =  2Dd(l  +  cos9) 
it  is  seen  that  as  0  decreases  from  m  to  0  the  value  of  L  de¬ 
creases  and  consequently  the  contribution  of  the  observation  at 
Pa  is  more  and  more  effective  in  reducing  the  error  of  estimate 
at  Pj  .  The  decrease  of  9  is  stopped  at  9  =  m/2  to  avoid  the 
trivial  case.  At  9  =  tt/2,  the  observation  at  P3  makes  its 
largest  contribution  to  the  estimate  at  P'  (within  the  bounds 
Imposed  on  the  problem).  Consequently,  where  the  "normal  law" 
decay  contributes  nothing  the  "simple  exponential"  decay  con¬ 
tributes  the  most. 

To  summarize  the  situation,  the  answers  to  the  best  loca¬ 
tion  of  the  second  observation  are  exactly  opposite  to  each 
other  depending  on  whether  the  variation  of  the  correlation 
coefficient  is  a  "normal,  law"  or  "simple  exponential"  decay. 

Idealized  contours  of  the  root  mean  square  error  of  esti¬ 
mate  are  shown  in  Pig.  VII I -4  . 


vm-i9 


Fig,  VIII-4  Contours  of  fz3  (Schematic)  in  the  Two  Examples,  k  =  1 
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ill )  Some  Limit  Relations 

Equation  (2)  may  be  expressed  in  the  form 


Now  consider  P2 


Pa  and  let 


ray  =  riy  +  |Vrx  y  |  *d -cos  (8  -co)  + - 

where  |Vrx y |  is  the  magnitude  of  the  gradient  of  the  correla¬ 
tion  surface  and  cp  is  the  direction  of  the  gradient  at  Px  (this 
correlation  is  thought  of  as  centered  at  P')  while  9  is  the 
azimuth  of  P2  from  Px  .  Also  let 


1-13  -  1  “  C|d|7  +  — 

(this  correlation  surface  is  thought  of  as  centered  at  Px  ) . 
Then 


It  then  follows  that  if  the  field  of  property  (wind  component, 
height  of  isobaric  surface,  etc.)  is  a  differentiable  function 
(y=2 ) ,  then  the  value  of  s2  has  a  limit  and  this  limit  is  less 
than  the  value  for  the  case  of  only  a  single  observation  at  Px 
If  the  field  of  property  is  not  differentiable  (0  <  y  <  2), 
then  the  information  provided  by  the  additional  observation  at 
Pa  approaches  zero  as  Pa  approaches  Pj  . 
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The  situation  described  in  the  above  is  quite  ideal.  Even 
in  the  differentiable  case  (in  the  limit  P8  -»  Px )  all  that  one 
has  is  two  observations  to  average  and  thus,  perhaps,  reduce 
the  error  of  observation  a  little.  The  mathematical  formulation 
above  is  valid  only  if  it  is  considered  possible  to  measure 
gradients  with  infinitely  close  observations. 

Note  that  the  limit  for  y  =  2  is  a  function  of  (9-co),  the 
orientation  of  the  line  ?!  P8  v/lth  respect  to  the  line  PlP/, 

The  function  ss  is  then  discontinuous  at  the  point  Px  and  the 
function  value  there  may  range  from  1  -  rn2  to  1  -  rn  2  -  |Vrx  ,  | 2 
depending  on  the  direction  in  which  P8  approaches  Px .  If  this 
is  along  the  gradient  of  the  correlation  surface,  the  smaller 
value  is  obtained j  if  perpendicular  to  the  gradient,  the  larger 
value  is  the  limit.  This  corresponds  to  the  directions  of  maxi¬ 
mum  and  no  additional  Information  in  part  (l)  above  and  general¬ 
izes  the  analysis  there  into  the  case  v/here  the  contours  of  the 
correlation  surface  need  not  be  circular. 

Consider  the  coefficients  ax ,  a8  as  P8  -*  Pt  .  Then 

a  =  u1vCldly  -  IVr^  > 1 ‘d • cos  ( 9  ■— cp)  +  - 

1  2C|d|7  +  --- 

r.yC|dly  +  |7ri  y  1  .d  .cos  (8 -co)  + - 

a®  2C  jd  |  y  + - 
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so  that 


a,  =  -  °°, 


aa  =  +  ", 


0  7^  cp  +  yp 


6  =  cp  +  p-,  l<ys2, 


cos  (8' 
2C 


y  =  i  , 


IJLX.  + 
2 


-  0  , 


0  -  CD  » 


8  /  cp  ,  0  <  y  <  1 , 


o  -  1  J- 

da  ?)) 


=  0  , 


The  behavior  of  the  coefficients  Is  highly  dependent  on 
the  direction  of  approach,  9,  and  on  the  value  of  y.  The  most 
irregular  behavior  occurs  In  exactly  the  case  of  most  Interest, 
y  •=  2 ,  8  =  cc,  cD  +  fr. 

The  reason  for  this  behavior  is  seen  by  rearranging  the 


terms  of  the  regression  into 


h'  =  (ax+as  )hj_  +  (aad ) 


h=  -hi 


and  nothing  that 


at  “t  a8  — 


ri>  11 1/^L.-— 

\d\y  +  — 


a-d  = 


VrT  y  ]d5cos  18 -cp)  +  riyCldK-1'  + - 

2C|d|r  +  --- 


Then 


+a  )  =  v  o  ^  v  ^  ° 
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while 


lim 

d-0 


(aa *<0 


=  i^i-Uos(9-co), 

-  0, 


y  =  2  , 

1  <  y  <  2 . 


In  the  limit,  the  second  term  of  the  regression  pertains 
to  the  directional  derivative  of  the  field  of  property  in  the 
direction  Px  Pa  while  the  first  term  still  pertains  to  the 
property  itself  at  Px 


h' 


bjhj  +  ba 


where 


2  (ai+as)>  ba  =  (a8-d) 


Then  b1  =  r1 y ,  ba  =0  for  0  <  y  <  2  so  that  the  situation  is 
unchanged  from  that  considered  originally,  i.e.,  for  a  non- 
dlfferentiable  field  of  property,  the  measurement  of  a  property 
gradient  (3h1/ad)  contributes  nothing  to  the  estimate  at  P' . 

For  y  =  2,  then  bt  =  rl  y  and  b8  =  (  |^rly  |/2C)cos  (9-co)  from 
which  it  follows  that  for  a  field  of  differentiable  property 
the  estimate  at  P'  is  Improved  by  the  measurement  of  the  direc¬ 
tional  derivative  of  property  and  that  the  amount  of  improve¬ 
ment  of  the  estimate  at  P'  is  dependent  on  the  direction  in 
which  the  derivative  is  measured,  a  maximum  if  on  line  PlP/, 
and  zero  in  the  direction  perpendicular  to  PjP'.  In  the  process 
of  taking  the  preceding  limits,  explicit  expressions  for  the 
regression  coefficients  have  been  derived  for  the  case  of  a 
differentiable  field  of  property  in  terms  of  the  characteristics 
of  the  correlation  function  of  this  property. 
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iv)  Observations  with  Errors 

The  material  in  the  preceding  section  is  of 
theoretical  interest  but  the  conclusions  are  not  particularly 
helpful.  The  relations  do  tell  that,  for  example,  if  we  are 
to  estimate  height  at  P'  from  height  at  Px ,  it  would  be  better 
to  measure  the  wind  component  perpendicular  to  PXP'  at  than 
make  another  height  measurement  anywhere.  On  the  other  hand 
if  the  field  of  property  is  wind  component,  we  cannot  measure 
wind  component  gradient  and  the  conclusions  are  silly  because 
as  P3  approaches  Pt  the  difference  in  wind  component  becomes 
lost  in  the  error  of  observation  or  the  effect  of  small  scale 
eddies.  To  considex-  these  effects,  the  correlations  are  de¬ 
graded  by  a  factor  1c,  0  <  k  <  1  which  measures  the  effect  of 
such  "errors".  Then  (2)  becomes 


v/here  the  correlation  coefficients  are  still  the  "theoretical" 
ones  (i.e.  without  "errors"  of  observation  considered ) .  In  the 
limit  as  Ps  -  Pa  there  are  no  indeterminate  expressions  and  the 
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relations  are 


„2  •)  2  k  2 

s  1  "  l+kri*  1 


ai  ”*  l+k1"1  s' * 


a3  “*  l+k^ 


so  the  limiting  form  of  the  regression  is 


h'  =  2k 


1+k 


r{ 


where  hx  and  h3  are  two  Independent  measures  at  Px . 

The  best  location  of  Pa  would  minimize  s2  or  maximize  the 


term 

p  =  ~  kr-lgriy)g 

1  -  k2rla2 

in  v/hich  the  co-ordinates  of  P3  appear.  To  determine  the  loca¬ 
tion  of  Pa  requires  a  detailed  knowledge  of  the  correlation 
functions  concerned.  Consequently,  the  specific  examples  will 
be  used. 

For  the  "normal  law"  correlation  surface 

P  -  P.  -  [e^/2  -  ke-"a( ■“•■«*)/»]  V(x  -  k’e-^2) 

=  e-as(D=4d2)/2  [^Dd-cose  .  k]s/(1  .  k£e-a2'J2) 
There  v/311  be  a  line  along  which  Fn  =  0  given  by 


d  -cosQ 


log  k 
asD 


This  line  will  be  perpendicular  to  PXP'  and  will  cut  it  a  dis¬ 
tance  | log  k|/a2D  from  Px  on  the  side  opposite  P'.  This  corres¬ 
ponds  to  the  line  through  Px  perpendicular  to  Px  P'  in  the  ca  se 
of  k  =  1.  It  is  the  line  along  which  an  observation  at  P8 
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contributes  no  information  to  the  estimate  at  P'  in  addition  to 
that  given  by  the  observation  at  Pa  .  The  point  at  0  =  ir  where 
this  line  intersects  PjP'  may  be  thought  of  as  a  point  where 
the  errors  of  observation  bring  about  so  much  error  in  estimating 
the  derivative  of  property  that  this  part  of  the  estimate  cancels 
the  advantage  of  using  two  independent  observations  instead  of 
one. 

The  quadratic  term  in  the  expression  for  Fn  (for  fixed  d) 
will  be  a  maximum  when  9  =  0  or  t r.  We  consider  only  9  -  tt  since 
for  0  =  0,  P8  lies  between  and  P' .  Then 


p.(a,ir)  =  e-a2(D2M  =  )/2[e-t<2M  .  k72/(! 


and 


3d 


2^-a2(Ds+d2)/2£g-a2Dd 


a  e 


ill 


fl  -  k*e-“*da)2 


k2e 


x  [cl(e-°,SDd  -  k)(l  +  'A*0*'3’)  +  2Dc-“eDd(l  -  k=e-“2<li!)]  . 
The  best  location  for  P3  is  given  by  the  zero  of  the  factor  in 
square  brackets  above.  This  root  (located  on  P^'  extended) 
is  a  true  minimum  of  s2,  a  best  location  for  Pa  (see  Pig.  VIII-5). 

In  the  simple  negative  exponential  form  of  the  correlation 
surface 


F  =  F, 


[e‘dL  -  ke-0(D4d’]2 

1  - 


where  L2  ^  D2  +  d2  -  2Dd-cos0,  It  is  seen  at  once  tiiat  the 
numerator  of  Fe  can  now  never  be  zero  since  L  s  D-fd  and  the 
first  term  will  always  be  greater  than  the  second  term.  Con¬ 
sequently,  a  second  observation  point  even  on  PjP'  extended  will 


Pig.  VIII-5  Contours  of  /s®  (Schematic)  in  the  Two  Examples,  k  <  1 
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provide  some  additional  information  -go  Improve  the  estimate  at 
P'.  However,  the  line  PXP'  remains  the  worst  location  for  a 
second  observation  point  only  now,  as  P8  recedes  from  Px ,  the 
value  of  ss  increases  from  1  -  2karxy/(l+k)  at  Px  to  1  -  karx, 
at  large  distances. 

The  situation  is  shown  graphically  in  Pig.  VIII-5  where 
figures  corresponding  to  those  of  Fig.  VIII-4  are  shown  but  in 
which  k  <  1.  In  Pig.  VIII-5a  it  is  seen  that  the  knot  point 
at  Px  Is  no  longer  present  and  that  the  loops  that  extended  to 
the  left  of  Px  are  now  replaced  by  a  center  where  the  estimate 
at  P'  is  best.  In  Pig.  VIII-[3b  the  geometrical  situation  is 
more  or  less  unchanged  as  compared  with  Fig.  VIII-4b.  Instead 
of  a  ridge  of  constant  value  for  the  error  of  estimate  extend¬ 
ing  to  the  left  of  P'  the  "ridge"  has  become  somewhat  "sloping". 

In  the  case  of  Pigs.  VIII~4b  and  VIII-5b  there  the  line 
PxPa  perpendicular  tc  PXP'  contains  two  points  (one  on  each  side 
of  Px )  where  the  location  of  Pa  is  best.  This  location  has  not 


been  discussed. 


...lUUlli.  ill  silltiitllf-  i- 
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B .  WIND  SHEAR  VECTOR  CORRELATIONS 

The  part  of  the  vertical  wind  shear  that  is  due  to  the 
synoptic  scale  variability  of  wind  may  be  obtained  from  the 
geostrophic  wind  shear  equations  in  the  form 


-Vo  -\aU 

As*  =  Aa7 


g  ST 

Te  ‘  ay  ’ 


t,  o  _  _  _°L  5T 

ASy  -  Agz  ~  3x- 

Where  T0  is  the  temperature  of  the  standard  atmosphere  at  the 
level  where  the  height  increment  is  measured  (assumed  constant 
in  the  horizontal). 

These  relations  are  identical  to  the  geostrophic  wind 
equations  with  the  exception  that  g/T„  now  replaces  g  and 
temperature  now  replaces  height  of  the  isobaric  surface.  Note 
that  x  and  y  are  to  be  measured  on  the  isobaric  surface  so  that 
the  gradient  of  temperature  is  not  in  the  true  horizontal. 

These  relations  may  be  used  to  construct  the  components 
of  the  wind  shear  correlation  tensor  in  terms  of  temperature 
correlations  in  exactly  the  same  way  that  components  of  the 
wind  correlation  tensor  were  expressed  in  terms  of  the  height 
correlations . 
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Trie  standard  deviations  of  the  wind  shear  correlations  may 
be  expressed  as^^ 

?\2(s7a')  =  (^-)  [  |VaTj2sin2i  +  2(ba  +  ba sina60  )°TS]> 
A2(“y’s‘)  =  [|7oTl2cos2<j<  +  2(bj  +  bscos290  )aT2], 

A2(st  sy )  =  -  [  |7oT|2sin^cos>l!  +  2basinQ0eos90  , 


where  it  has  been  assumed  that  the  temperature  correlation 
function  on  an  isobaric  surface  may  be  written  in  the  form 

RTT  =  1  "  [bi  +  b3cos2(9-e0  )]r2  +  — 

where  9  is  the  orientation  of  P'  with  respect  to  P  and  r  is  the 
distance  PP'.  The  parameter,  L^,  is  a  size  parameter  which 
measures  the  average  dimensions  of  temperature  systems  on  an 
isobaric  surface  (L^"2  -  2b*  ) . 

The  above  relations  have  not  been  used  to  estimate  the 
standard  deviation  of  the  wind  shear  component  due  to  the  fact 
that  data  on  the  temperature  correlations  are  not  generally 
available.  If  the  structure  of  the  temperature  correlations 
in  the  horizontal  were  known,  the  Information  provided  would 
be  useful. 

The  relations  are  not  as  complex  as  they  appear  since  they 
contain  several  small  terms.  The  value  of  |7o^,|  is  reasonably 
small  and  the  value  or  ba  may  also  be  small.  In  the  later  case, 
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this  amounts  to  the  assumption  that  two-point  temperature  cor¬ 
relation  contours  are  essentially  circular. 

The  analysis  of  the  components  of  the  wind  shear  correla¬ 
tion  tensor  proceeds  exactly  as  that  of  the  wind  component  cor¬ 
relation  tensor.  The  tensor  terms  should  have  almost  the  same 
"structure",  i.e.  the  contours  of  constant  shear  correlation 
tensor  component  should  be  nearly  coincident  with  the  corres¬ 
ponding  wind  correlation  tensor  component  if  the  assumption  that 
the  temperature  correlations  would  be  like  the  height  correla¬ 
tions  is  made.  This  is  not  necessarily  the  case. 

In  the  case  of  the  height  correlations,  the  assumption 
that  height  was  a  differentiable  function  was  made.  This 
assumption  is  somewhat  less  plausable  in  the  case  of  tempera¬ 
tures.  The  nature  of  the  correlation  structure  is  intimately 
dependent  on  the  character  of  the  field  of  property  (differen¬ 
tiable  or  not  and  if  not  the  nature  and  "density"  of  its 
Irregularities ) . 

It  was  seen  in  Chapter  VII  that  the  standard  deviations  of 
wind  shear  may  be  expressed  as 

s*  *  =  (IsH  +  2au2as  , 

Sy2  =  (fz2-)  +  2avat>8  > 

where  aa  and  ba  are  the  coefficients  of  the  second  order  terms 
of  the  correlation  of  wind  component  as  a  function  of  height: 

ruu  =  1  -  asz2  + - , 

rvv  =  1  -  b8z2  +  - . 
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c.  "universal"  extrapolation  and  forecast  methods 

The  usefulness  of  correlation  structure  In  the  extrapolation 
of  weather  parameters  and  In  forecasting  those  parameter  values 
may  be  illustrated  by  an  analysis  of  an  elementary  least  squares 
extrapolation  or  interpolation  problem.  Consider  observations 

at  P1  , - ,  Fn  where  x1  , - ,  xn  are  observed  and  let  P'  be  the 

point  at  which  an  extrapolated  value,  y,  is  required  and  assume 
that  the  extrapolation  or  forecast  is  to  be  a  linear  combination 
of  the  Xj  's  in  the  usual  way.  Then 

y  =  b1x1  + - +  bnxn  (4) 

(standard  variables  for  convenience) .  The  normal  equations  are 


+ 

baria 

+ - + 

b»rln 

>* 

Sn 

II 

a  i 

+ 

ba 

+ - + 

bn  ra  n 

-  ray 

'n  1 

+ 

barns 

+ - + 

b»  r» 

II 

a 

X 

where  ru  =  rj t  are  the  correlation  coefficients.  If  these  are 
rewritten  in  general  notation, 

-  gi 

aij  =  riJ>  (rSl  “  Ci  =  ri  y 

then  the  solution  for  the  coefficients  is 

bj  =  gta1J 

where  a1J  represents  the  element  in  the  i'th  row  and  j 1 th  column 
of  the  inverse  of  the  matrix  of  coefficients.  Now  since 

y  =  bjXj 

then 

y  =  Xjgtal)  =  gj  (xja1J  ) .  (5) 
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Consider  now  what  is  involved  in  the  various  quantities 
that  have  been  used.  First,  the  values  of  xt  are  the  observa¬ 
tions  made  at  a  number  of  fixed  points,  P: , ,  Pn  .  Second, 

the  a1J  are  the  elements  of  the  inverse  of  the  matrix  at i  which 
are  in  turn  the  n2  correlation  coefficients  that  relate  the 
observations  at  the  same  set  of  points.  These  quantities  are 
known  a  priori  from  the  past  record,  or  from  the  structure  of 
two-point  correlations  of  the  quantities  concerned.  As  soon 
as  the  quantities,  xt ,  are  kn own,  then  the  terms  (xja1J)  may  be 
computed.  Third,  the  correlations,  gt  =  rly,  are  the  two-point 
correlations  relating  xt  at  P}  with  y  at  P' ,  In  view  of  the 
fact  that  one  is  dealing  with  a  field  property,  these  correla¬ 
tion  coefficients  may  be  considered  as  known  functions  of  the 
position  of  P'  with  respect  to  Ps .  As  a  consequence,  the  loca¬ 
tion  of  P'  may  be  considered  variable  and  equation  (5)  may  be 
considered  as  a  least  squares  interpolation  formula. 

In  the  elementary  form  considered  above,  the  interpolation 
aspects  of  the  relationship  are  emphasized.  The  same  relation 
may  be  used  as  a  forecast  formula  if  the  time  lag  relations  are 
included  in  the  structure  of  the  factors  gt  -  riy. 

In  the  example  above  only  observations  xt  at  point  Pj  were 
considered.  This  was  to  simplify  the  mode  of  expression .  The 
technique  is  quite  applicable  to  the  case  in  which  several 
quantities  are  observed  at  each  Pt .  For  example,  height  and 
two  wind  components  may  be  used  and  from  these  estimates  at  P' 
of  height  or  either  (or  both)  wind  component  may  be  made.  In 
such  an  instance  there  would  be  n  observation  points  but  3xn 


predictors  and  three  equations  leading  to  the  prediction  of 
three  quantities.  In  this  particular  instance  the  three  sets 
of  equations  would  have  a  common  matrix  of  coefficients  and  would 
differ  only  in  the  terms  gj  =  rt r .  Thus  once  the  common  inverse 
elements  of  (2)  are  obtained,  there  would  be  three  sets  of 
multipliers,  gt ,  one  for  each  of  the  predicted  quantities. 

Several  applications  may  be  made  to  problems  of  practical 
importance.  Some  of  these  are  listed  below: 

1 .  Checking  for  Erroneous  Observations 

In  this  case  the  value  at  P'  is  observed  and  a  value 
interpolated  from  surrounding  observations.  If  the  value  ob¬ 
served  differs  from  the  interpolated  value  by  more  than  a  pre- 
asslgned  amount,  the  observation  may  be  rejected  as  erroneous. 
Somewhat  similar  methods  are  used  in  actual  practice,  but  are 
usually  not  based  on  the  statistical  structure  of  the  field  of 
property.  For  example,  a  relation  of  the  form  j^l  +  (d/L0)sJ-1 
is  used  where  the  value  of  L0  is  allowed  to  Increase  as  further 
refinement  of  the  test  is  required.  This  formula  is  acceptable 
as  far  as  functional  form  is  concerned,  but  to  allow  L„  to  ex¬ 
ceed  the  value  imposed  by  the  correlation  coefficient  structure 
amounts  to  smoothing  the  data  beyond  that  which  could  be  natu¬ 
rally  present.  Other  smoothing  techniques  use  a  form  in  which 
the  exponent  of  d/L0  is  4  instead  of  2.  Such  formulas  force  a 
smoothing  (or  rejection)  which  is  quite  impossible  since  it  can 
be  easily  shown  that  the  exponent  cannot  exceed  2. 
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2  .  A  Sage  Problem 

In  transforming  data  from  synoptic  observations  Into 
values  on  a  three-dimensional  grid  for  the  purpose  of  interceptor 
missile  control,  the  preceding  type  of  an  extrapolation-interpo¬ 
lation-forecast  method  is  unusually  useful.  It  lends  itself 
readily  to  computer  applications  in  that  the  only  stored  elements 
need  be  the  observation  point  and  grid  point  co-ordinates  to¬ 
gether  with  the  simple  two-point  correlation  coefficient  formula. 
Observations  fed  into  the  system  can  be  transferred  to  the  grid 
by  direct  interpolation -extrapolation  and  the  time  element  in 
the  gt  *s  may  oe  used  to  up-date  the  information  hour  by  hour 
(or  minute  by  minute)  until  new  observation  information  is  in¬ 
troduced  . 

3.  Artillery  Problems 

For  Army  operations  in  which  data  are  available  over 
a  corps  area,  the  above  type  or  extrapolation  provides  an  easily 
computerized  procedure  for  arriving  at  expected  parameter  values 
in  time  and  space,  not  only  in  the  area  of  observations,  but 
also  extending  into  the  surrounding  areas  .  The  method  is  in¬ 
dependent  of  the  location  of  the  observation  points  so  that  it 
is  applicable  to  a  "fluid"  situation  in  which  the  observation 
points  may  be  changing  from  day  to  day. 

4 .  Fleet  Operations 

Extrapolation  from  observation  points  that  are  in  a 
state  of  continual  change  may  be  handled  with  nearly  the  same 
ease  as  the  case  of  fixed  observation  points.  The  major  require¬ 
ment  in  this  situation  is  merely  that  the  observation  must  be 


I  cl, i  :  .4‘i  lift.  ftb.i.Aif'im  i  ttM'W  • 
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accompained  by  the  co-ordinates  of  the  point  of  observation 
when  introduced  into  the  computing  procedure.  In  the  other 
case  (fixed  observation  points)  these  co-ordinates  were  con¬ 
sidered  as  already  tabulated  in  the  procedure. 
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APPENDIX  A -I 

WIND  COVARIANCE  TENSOR  COMPONENTS  IN  TERMS  OF 
HEIGHT  CORRELATIONS  AND  STANDARD  DEVIATIONS 

The  equations  (7)  of  Chapter  III  are  derived  here.  The 
notation  used  is  the  same  as  that  of  the  preceding  chapters 
with  additions  where  necessary.  Let  the  geostrophic  wind  equa¬ 
tions  be  written  as 

ela„\,u,  =  -  g.|a-,  1  =  1,  2,  (1) 

where  A3  =  >,  =  Coriolis  parameter,  uk  =  wind  component,  g  = 
acceleration  of  gravity,  h  =  height  of  the  isobaric  surface, 
xt  =  horizontal  co-ordinate  on  an  isobaric  surface,  el3lt  =  +  1 
if  13k  are  in  natural  order,  -  1  if  in  reverse  order,  0  if  other¬ 
wise.  The  summation  convention  on  repeated  subscripts  is  being 
used.  Here,  and  in  what  follows,  the  subscript  conventions  of 
ordinary  cartesian  tensors  are  employed.  The  reason  for  this 
becomes  immediately  evident  since,  if  the  calculations  are 
carried  out  in  component  form,  the  bulk  of  the  expressions  is 
extremely  cumbersome. 

Since  the  equations  (1)  are  linear  in  uk  and  h,  the  symbols 
uk  and  h  may  be  understood  to  represent  mean  values,  observed 
values,  or  departures  from  the  mean.  In  the  following,  the  sym¬ 
bols  will  be  accepted  as  departures  from  the  mean  to  simplify 
the  ideas,  means  will  be  considered  as  taken  over  a  sample  of 
years  or  seasons  or  days  at  a  fixed  point  (geographical  location 
and  isobaric  level)  though  this  restriction  is  somewhat  stronger 
than  is  actually  required. 


"  Y!  7!”  T1' Jn  T’1.'  p7?TiT 
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The  equations  (l)  are  equally  valid  at  all  points  so  that 
the  equations  are  understood  to  refer  to  a  point  P  as  written 
and  to  a  point  P7  if  the  symbols  are  suitably  modified  as  in 

ej  3  =  -  S-fpv  J  =  1,  2.  (!') 

On  multiplying  equations  (l)  and  (l7)  together  and  taking 
mean  values,  indicated  by  the  bar  over  the  symbols. 


eiakeja^s^fukul)  -  ga  ~>  o  -  1>  2,  (2 ) 

where,  on  the  right  of  (2),  use  has  been  made  of  the  fact  that 
h  is  a  function  of  and  h7  is  a  function  of  xf  only.  The 
correlation  function,*  Rm,  is  defined  by  the  relation 

<?a'Rhh  =  (hPJ  (3) 

where  a  is  the  standard  deviation  of  height  at  P  and  o'  is  the 
standard  deviation  of  height  at  P7.  As  pointed  out  in  Chapter 
III,  the  correlation  function,  R61,  is  a  function  of  the  loca¬ 
tion  of  P  and  P7.  On  the  other  hand,  a  is  a  function  of  the 
location  of  P  only  and  a'  of  P7  only.  The  differentiation  of 
the  right-hand  side  of  (2)  takes  this  situation  into  account: 


3b  -  3g/n 

wr wf* 


The  height  correlation  function  may  be  thought  of  as  pri¬ 
marily  a  function  of  the  location  of  P7  with  respect  to  P  which 
involves  the  location  of  P  itself  in  a  secondary  sense.  For 


*The  symbol  Rhh  is  used  instead  of  the  usual  symbol  for  such  a 
correlation  coefficient,  rhjl»  since  the  lower  case  r  will  be 
used  later  as  the  radial  distance  from  P  to  P7. 


iNryi-  !H|ipM 
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example,  if  (r,8)  are  the  polar  co-ordinates  of  P'  with  respect 
to  P  and  if  one  expresses  the  height  correlation  fiinction  in 
terms  of  a  series  of  powers  of  r  with  coefficients  as  trigono¬ 
metric  polynomials  in  0  such  as 


Rhh  =  1  -  [ai  +  ascos2(e-90 )]r2  +  —  (5) 

the  parameters  a1 ,  a3 ,  60,  would  be  functions  of  the  location 
of  P  itself.  The  values  of  these  parameters  might  be  thought 
of  as  being  slowly  varying  functions  of  the  location  of  P.  Let 
the  co-ordinates  of  P'  with  respect  to  P,  in  cartesian  form  be 
given  by  where 

Xj'  =  x,0  +  h,  (6) 

Xt  -  X*0. 

Then 


_b _ 3_  _J_  _5_  _  a 

~  9xt0  “  35 1'  3x, 

so  that  (using  (3),  (4),  and  (7))  (2)  becomes  (dropping  the 
superscript  0  as  being  superfluous) 


(7) 


ei3kej3iA,^(u1{ui)  =  g2{|2_.|£^6ll  +  o'&L-  ( 

aCT  a5t3§,  Uxj  3Xi  CT  3?jax1J;* 

The  final  term  in  brackets  of  equation  (8)  involves  the 
fact  that  the  correlation  function  Rbh  is  a  function  of  the 
location  of  the  point  P.  The  relative  location  of  P'  with  re¬ 
spect  to  P  is  not  shown  by  the  notation  in  the  derivatives  of 


o  and  o',  the  variables  locating  P  and  P'  separately  are  re¬ 
tained  to  emphasize  the  functional  dependence  of  these  q\jantities. 
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As  noted  by  (5)>  the  correlation  function  has  a  tendency 
to  have  an  axis  of  symmetry.  If  such  a  direction  is  taken  as 
that  of  60,  then  it  is  convenient  to  express  Rhh  in  terms  of 
co-ordinates  (r,s),  s  =  r*cos(8-0o)  (see  Fig.  III-2).  Let  Xj 
denote  a  unit  vector  in  the  direction  of  90  and  let  be  the 
unit  vector  in  the  direction  PP'.  Then  s  =  SjXj,  r2  =  = 

r2^2,  and  then 


9?!  Sr'”’1 


1 


and 


Var2 
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_  t 
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ar 


^vivi  +  F* 


bH 


Jub. 


ar 


«1J 


gfe*''.*. 


+  VjX,  ) 


(9) 


(10) 


where 


fi,,  =  1  if  1  =  j  and  6ij  =  0  if  i  /  J. 


First,  the  expressions  (9)  and  (10)  are  substituted  into 
(8),  then  the  identity 

ei3iteJ3i(ufcuJ[)'  -  6 1  j  ( uk  u^  )  -  (u4u/) 

is  used  to  obtain  the  wind  component  covariances  on  the  left- 
hand  side  without  the  awkward  triple -subscript  expressions 
present,  and  finally  the  derivatives  of  the  standard  deviations 
of  height  are  replaced  by 


fr=  l7a||j'1'  !xj" =  iVc?/i^ 

where  | Vcr  | ,  | Vcr 7 1  are  the  magnitudes  of  the  gradient  of  the 
standard  deviation  of  height  and  m  ,  [i.j  are  the  unit  vectors 


so' 
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in  the  direction  of  this  gradient  at  P  and  P'  respectively.  The 
above  three  steps  are  performed  with  no  conceptual  difficulties, 
but  the  algebra  is  admittedly  tedious.  The  final  result  is 


(id 


+  •^ttal7a/|vJ^/  -  o'  IVcIvjHj] 

+  ^[gl7g'^lx»  ~  cr^  lVal^iPLj  ]  -  |Va|  17a7  |RhhM-:  M-i 
u {- 


+  6 

+  IVol  iVa7  lRahHi.  M*'  -  ^^glVcr'Kh.'  -  1 W  |  vfc  ^  ] 

-  ^4gl7g/  1A*P*  -  g'lVgl\^t]  • 

The  terms  of  (ll)  are  underscored  in  the  same  way  as  those  of 
Chapter  III,  equations  (7).  The  significance  of  the  underscoring 
is,  briefly, 

Terms  remaining  if  R„  n  =  Rhh(r)  and  |Va|  ■  0 
ru, ,  r jru-L.  Additional  terms  if  Rhh  =  Rhh(r,s),  |Va  |  s  0 
— —  Additional  terms  if  Rhi,  =  Rhh  (r),  |7q|  yo 
— -,r— rr  Interaction  terms  when  both  Rhh  =  Rhh(r,s),  |Vcr|  /  0 
When  Rj,!,  =  Rbh(r),  |7aj  e  0,  the  wind  covariance  tensor  is  given 
by 


[(ay .  .  6,,^* 


(12) 
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This  tensor  form  is  that  of  a  two-dimensional  isotropic  turbu¬ 
lence.  This  situation  is  discussed  in  some  detail  in  Chapter 
III. 

When  Rhh  =  Rb)l(r),  |Vo|  ^  0,  the  wind  covariance  tensor 
becomes 


pi(Ulu;)  =  -  r^-)v,v,  +  -  CT'l7a| V,|X4 

-  |Vff||ttr'|RnH,n»'  +  +  |Va |  |Va'  |Rhh u* pk' 

L  Br* 

-  -^[crlVa7  Km*'  -  a'l^KM,]}  (13) 

where  the  effect  of  curvature  or  unequal  spacing  of  the  contours 
of  o  =  constant  is  retained  by  the  distinction  betv;een  a  and  a' , 

M-i  and  m-j,  etc.  The  last  term  in  braces  may  be  considered  as 
quite  small.  The  expression  ( 13 )  resembles  the  covariance  tensor 
for  a  two-dimensional  cylindrical  or  axlsymmetrlc  turbulence, 
except  that  the  factor  in  brackets  in  the  second  term  on  the 
right  contains  a  negative  sign  instead  of  a  plus  sign. 

If  one  considers  Rtoh  -  Rhh(r,s),  |Vcr|  a  o,  then  the  wind 
covariance  tensor  is  given  by 


(1^) 


which  is  precisely  the  form  of  a  two-dimensional  homogeneous, 
axlsymmetrlcal  turbulence. 


A -I -7 


THE  COVARIANCE  TENSOR  IN  LONGITUDINAL  AND  TRANSVERSE  COMPONENTS 
To  obtain  the  covariance  tensor  in  terms  of  longitudinal 
and  transverse  co-ordinates,  the  rotation 

u^  =  us  VI 
Ut  = 

is  used.  Let  the  direction  of  Vj  be  0,  of  be  0  0 ,  and  of  (ij 
and  u..'  hp  ijr  and  $ 7 .  Then  multiplying  (ll)  by  vlvi  and  carrying 
out  the  consequent  subscript  summations 

^V?  -  -  +  (15) 

+  |Vcr  1  |Va7  |Rhh^cos  (*-!>'  )  -  cos  (9  )  cos  (9 -\|; ' )  J 

-  ^||*-{cf|Vo7  |[cos(0o -i'7 )  -  cos(9-9e  )cos(9-^,)J 

-a' |7a|[cos(90 -■!/)  -  cos  (9  -0O  )cos  (e -\Jf )]} . 

To  obtain  (u^u^),  multiply  (ll)  by  e13vvJej3l!v1  and  carry  out 
the  summations  to  obtain 

17^  +  2|?i^s(9-e0)  +  ~tA^os2(e-0o 

(16) 

+  |Vaj  |Va7  |R„  „[cos  (^-ii; 7 )  -  sin(9 -i|i )sin(0 -iji ' )] 

-  ^(a|7a,|cos(e-i'/)  -  a7  |7a |cos (e ) ] 

"  |^o,/  j[cos(0e  ->i7)  -  sin(0-0e  )sin(9-'^7^] 

-  a'  jVc  jjjjos  (0„  -i‘)  ~  sin  (9 -0 0  )sin (8 )  j|. 

The  covariance  (u^u£)  is  obtained  by  multiplying  (l)  by  v1eJ3Itvk 


-  00 


■ 
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and  is  given  by 


77' 


=  ^^-^^oo3(6-eo)sin(e-0o)  +  ^|tcin(B-90)j  (17) 

-  \Va\ |Vcr'  |Rhhcos(e-*,)sin(6-iif)  -  o#- |Vff|^^fiin(0-*) 

cos  ( 9 '  )sin(e  -90 )  -  a' |W|cos(9-e0  )sin(e-il()^, 


The  covariance  (u^u'),  from  multiplying  (ll)  by  elRicv^\>s  is 


given  by 


77/ 

„2 


(Vl)  =  aa'  B-a-R“cln(e-9a)cos(9-88)  +  |^*-cin(9-9, )]  (18) 

-  |Va|  |Va'  |RSbcos(9-'i)sin(9->j/)  +  ajVa'  l-^cin  (9 ' ) 

+  -^r^a |Va' 1 coo (9 -90  )sin (e -■') ' )  -  a'  |Va|sin(9-90  )cos  (e->jr)]. 


If  now  we  use  the  approximations 


*  7<^  -  X’aJ’  Lt 


£S_  *.  p,s 
7a.  “  tToT 
t  t 


where  L  is  a  size  parameter,  and  if  we  assume  that  $  a:  i[i 7  and 
i^H  -  |^pH  ,  then  (15)  through  (18)  may  be  expressed  in  terms 
of  correlation  coefficients.  Thus 


~i~  ~  r  •  +  ‘if^3 ln 5 (9  ~9»)  +  \^-\\  1  sin2 (6-*),  (15a) 
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-  |~|  Rbhsin(e-t)cos(e-^) 

+  1^1  [^ABin(e-'^)  +  Gin  (8  9  -fr)]  . 


(l6a) 


(17a) 


(18a) 


In  the  case  of  (17a)  and  (l8a),  it  is  seen  that  the  terms  are 
grouped  in  pairs  which  may  be  separated  in  the  expressions 


D 

¥7 


i£i  [J 


;in(e-t) 


(20) 


The  underscoring  of  the  terms  in  equations  (lRa),  (l6a), 
(19)  and  (20)  correspond  to  that  previously  used  on  the  terms 
of  equation  ( 1 1 ) .  The  principle  terms  appear  only  in  (15a)  and 
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(l6a),  but  each  of  these  includes  one  group  of  terms  due  to 
the  fact  that  Rbli  may  depend  on  -s  and  a  term  that  depends  on 
the  fact  that  |Va|  /  0.  The  expression  for  S  contains  terms 
for  each  of  the  latter  types.  The  expression  for  D  contains 
one  term  due  to  the  fact  that  |Va|  /  0  and  another  term  in  the 
interaction  between  l*o  i  and  Rhh(r,s).  Thus  in  (20)  the  expres¬ 
sion  for  D  vanishes  if  |Va j  a  o.  Note  further  that  (20)  is  the 
only  one  of  the  four  expressions  in  which  |^-|  appears  to  the 
first  power.  If  L*  |^-|  is  reasonably  small,  one  might  say  that 
(20)  contains  the  primary  effect  of  |Ve|  /  0  and  that  this 
effect  appears  to  "higher  order"  in  the  other  expressions.  If 


we  let 


where 


Rji a  =  expj^  -  i(asr2  +  02s2)J 


1  nS  ...  1 

V'  *  h* 


Rh*  =  i 


-  ft®  *  ft  -  ^X®]  *  - 

"  [v  +  (v  "  lt)cos2(c"9o)]%*  + 


which  introduces  elliptical  contours  of  Rhh  =  constant.  Then 


=  '  aS*exp[  “  $(aar2  +  02s2)] 

^  -  «2r2)exp[  -  *(a2r2  +  02s2)] 

=  a202rs-exp[  -  £(a2r2  +  02s2)] 

5-R|*  =  -  02(l  -  02s2)exp[  -  |(a2r2  +  02s2)] 

^  -  -  0ss.exp[  -  j(,V  +  £*»*)] 


-  a" -exp 
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so  that 

vu  -  V01*^1  +  fi--sinsCe-s. )  +  ^|2i|%in*(e-t)] 

-  r^|  Glng(e-e0)cosg(9-80)j-  -exp^-  £(aar2  +  £2s2)J, 
rtt  =  Lt2a2{[l  +  ^cos2(e-0e)  +  i|^|2cos2(0-t)] 

-  a2r2[\  +  ^cos2(e-0o)]  j-exp£-  £(a2r2  +  02s2)j 
S  =  -  LzLta^[\^-\  sin2(9-t)  +  02sin2  (e -0 0  )j 


*/32r2  £l  +  ^cos2(0-0o  )j  sin2(e-0„  )j-*exp[-  £(a2r2  +  j32s2)], 

D  =  -  2LALta2r  l^-l^sin(G-'Jf)  -  |scos(e-90  )sin(9o x 

?(a2r2  +  02s2)]. 


x  exp 


Since  '^ni  r  =1  r  -  1 

bince  r_<0  rJil  -  i,  r-0  rtfc  -  i, 


[Lj'oi2  l]  “  1  +  In5 (e -9a)  +  JfIitI  sina ( 9 , 

[htsus  "']  -  1  +  ^cos={8-9.)  +  Jy|S2.|%os*(e-t), 

which  indicate  a  slight  dependence  of  and  L^.  on  the  orienta¬ 
tion  angle  6,  the  azimuth  of  P'  relative  to  P.  Then 
f  L  2S4r2  1 

?SLl  =  [l - ^lp - sin22  (9  -90  )J  exp[-  i(a2ra  +  ,S2s2)],  (21) 

rtt  =  [l  -  Lt2a4r2(l  +  ^cos2(e-90 )]  exp^-  £(a£r2  +  j32s2)].  (22) 

From  (21)  and  (22),  the  effect  of  the  elllpticity  of  the 
height  correlation  contours  introduces  some  modifications. 
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These  would  have  complex  effects  on  the  shape  of  the  contours 
of  equal  correlation  coefficient.  The  exponential  term,  of 
course,  preserves  the  over-all  ellipticity  of  the  height  corre¬ 
lation  contours,  but  the  trigonometric  terms  in  the  first  factor 
would  induce  some  additional  irregularities  at  the  larger  dis¬ 
tances.  These  would  appear  at  those  distances  where  the  cor¬ 
relations  are  small  and  where  the  numerical  values  of  the 
observed  correlations  are  uncertain.  The  effects  appear  to  be 
due  primarily  to  the  ellipticity  of  the  height  correlation 
contours,  #  /  0.  The  effect  of  the  gradient  of  the  standard 
deviation  of  height  are  hidden  in  the  factors  and  L^. 

The  effects  of  non-homogeneity  (|7a|  /  0)  and  lack  of  isot¬ 
ropy  (j3  4  0)  appear  clearly  as  the  reason  for  3^0  and  D  jL  0. 

When  only  the  dominant  terms  are  retained,  e:  L2, 

S  =  -  [La|¥2-|asin2(e-*)  +  La02sln2(9-9O)] -exp[-  &  («2r2  +  J82s2)],  (23) 

D  =  -  2rL|^-|  *sin(9 -'Ji )  *exp^-  £  (a2r2  4  /32s2)^J,  (24) 

Prom  (23)  the  effect  of  |7a{  and  j3  are  combined  to  form  an 
expression 

3  =  -  A*sin2  (e-v)  *exp[  £(a2r2  4  j32s2)J, 

where 

A2  =  L4{[02cos29o  +  |^|2oos2*]2  4  [/32sin20e  +  |^-|  sin2*]2}, 

j32sin20 0  4  |^2.|  sin2i|f 
tancp  =  •  1  ■  —  1  •  a-"  . 

j32cos29  0  +  ]— |  cos2\Ji 
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The  function  S  is  discontinuous  at  r  =  0  and  that  the  contours 
of  constant  S  form  a  four-leafed  clover  pattern  with  S  =  0  on 
the  directions  0  =  cp  +  k^,  k  =  0,  1,  2,  3j  and  with  values  of  S 
alternating  sign  in  the  quadrants  between  them.  The  extreme 
values  of  S  are  at  the  discontinuity  point.  The  magnitude  of 
the  discontinuity  at  r  =  0  is  due  to  both  the  gradients  of  the 
standard  deviation  of  height  and  the  "ellipticity"  of  the  con¬ 
tours  of  the  height  correlation. 

The  expression  fur  D  in  (24)  is  dominated  by  the  factor  in 
gradient  of  the  standu;!  deviation  of  height.  If  the  height 
field  were  homogeneous,  unis  expression  would  be  zero.  The 
expected  contours  of  constant  D  form  a  pattern  of  two  centers, 
one  of  positive  value  and  the  other  of'  negative  value  separated 
by  a  line  along  which  D  =  0  that  lies  in  the  direction  of  (i, 
the  direction  of  the  gradient  of  tho  standard  deviation  of  height. 
The  extreme  values  arc  approximately  ±  2L;|^j-e‘'n  and  are  lo¬ 
cated  at  a  distance  of  r  =  L. 
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SPECIAL  CONSIDERATIONS  REGARDING  THE  DETERMINATION  OP 
L  AND  k  FROM  THE  UNIVERSITY  OP  WISCONSIN  DATA 

The  University  of  Wisconsin  wind  component  correlation 
tensor  data  requires  special  consideration.  The  majority  of 
the  data  on  wind  correlations  has  been  calculated  from  the 
actually-observed  winds,  or  at  least  from  data  that  purports 
to  closely  approximate  such  winds.  The  University  of  Wisconsin 
correlations  were  calculated  from  basic  height  data  using  height 
differences  between  points  10°  latitude  or  longitude  (generally) 
apart.  These  height  differences  were  multiplied  by  the  appro¬ 
priate  constants  to  convert  them  into  geostrophic  wind  compon¬ 
ents.  Consequently,  the  correlations  are  not  between  geostrophic 
wind  components  at  two  points,  but  between  average  geostrophic 
wind  components  over  two  (sizable)  intervals.  The  correlations 
thus  calculated  tend  to  be  larger  than  the  correlations  that 
would  be  obtained  were  it  possible  to  use  geostrophic  wind  com¬ 
ponents  at  a  point.  The  amount  of  the  difference  between  such 
averaged  components  and  point  components  is  an  Increasing  func¬ 
tion  of  the  length  of  the  interval  over  which  the  averaging  is 
done . 

The  correction  for  the  fact  that  average  instead  of  point 
winds  are  used  may  be  derived  in  the  following  manner.  Let 
u4  ,  u'  be  the  components  of  the  averaged  geostrophic  wind,  let 
u,  u7  be  the  components  of  the  true  geostrophic  wind,  let  the 
averaging  Interval  be  of  lengths  p  and  ?  centered  at  P  and  P7 
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respectively.  Then 


u. 


S.. 

A 


n 

2 


ri/2  ri/2 

H/w ay  -  */u  ay- 

-n/2  -n/2 


v.  = 


Ehl  -h -I 

A'  ? 


5/2 

dx 

A  ?  1  3x  dX 

-?/2 


?/2 


if 

-5/2 


v  dx, 


Then 


n/2  V/2 


u.  u; 


l 

rjrp 


(uu' )  dy  dy' . 


-r,/2  -riV2 


As  long  as  P  and  P'  are  separate,  the  covariance  of  the  averages 
(left  side)  is  nearly  the  same  as  the  average  of  the  covariances 
(right  side) 

u,  u,'  a-  uu ' 

However,  when  P'  coincides  with  P,  this  approximation  is  no 
longer'  valid.  (Actually  it  becomes  valid  as  P'  approaches  P, 
but  in  using  grid  data,  this  situation  is  seldom  encountered.) 

In  this  case  p  =  V  but  the  double  integral  form  is  maintained 
since  on  the  right,  the  points  at  which  u  and  u'  are  determined 
are  independently  distributed  on  the  interval  ( -ri/2,  +  r/2 ) 
about  P. 

uu''  -  crua'rttu 

-  tfu(y)cu'(y')i',uu(y>y') 

=  a*  (y)du'  (y')r(y'-y) . 


Now  let 


fi:  *•!  1J.  f« 
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Substituting  this  expression  into  the  integral  and  using  mean 
values  for  au  and  a' ,  one  has,  approximately 


+t|/2  ri/2 


(y'-y)  dy  dy' . 


-r/2  -tt/2 


Changing  variables  to  Z  =  y'  -  y 


_ _  r 

Hi 


h/2  -y+h/2 


(z)dz 


-n/2  -y-n/2 


dy. 


Note  that  ruu(Z)  is  the  transverse  wind  component  correlation 
function.  A  rough  approximation  for  this  function  appears  to 
be 

ru,(-  rtt)  —  ( 1  -  aaZ8)c"®  Z  /2. 

On  making  this  substitution  and  integrating 

-^r?/2)-  ^[f <«)]• 

Table  I  gives  values  of  f(ar|)  as  a  function  of  «r|. 


TA3LE  I 

VALUES  OP  f(ori). 


ar\ 

[f<<rO] 

0.0 

1.00 

0.2 

1.012 

3.4 

1.041 

0.6 

1.082 

0.8 

1.128 

1.0 

1.289 
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We  have  seen  that  L  “  has  values  In  the  range  of  500  to 
600  n  mi.  If  the  length  of  the  averaging  interval  is  10°  lati¬ 
tude,  then  ar\  is  in  the  range  0.8  to  1.0  and  the  standard  devia¬ 
tion  of  the  average  is  from  10$  to  30$  higher  than  the  standard 
deviation  of  the  point  average. 

In  the  case  of  the  meridional  wind  component,  the  value  of 
ap  is,  of  course,  modified  by  the  factor  coscp.  Consequently,  at 
60°  north,  the  underestimate  of  the  standard  deviation  by  using 
an  average  geostrophic  wind  is  much  less  for  the  meridional 
rather  than  the  zonal  component.  (From  the  above  table,  if 
L  =  8°  latitude,  the  excess  is  4$  for  the  meridional  component 
as  compared  with  13$  for  the  zonal  component  where  the  grid 
point  spacing  is  the  full  10°  of  latitude.) 

In  the  applications  to  wind  component  correlations  in  trans¬ 
verse  and  longitudinal  form,  this  correction  appears  in  the  stand¬ 
ard  deviations  (variances): 


uf2  =  (u2)cos20  +  2  (uv)slnQcos0  +  (v5')sln2e 

ufc2  =  (u5')sin29  -  2 (uv")sin0cos9  +  (v2)cos20 

as  correcting  the  first  and  last  terms  on  the  right.  If  several 
approximations  are  made,  u2  es  v2  »  uv",  then 


<v>. 

cos29  , 

sins9 

f  (an) 

f  (apcosco) 

sin20  , 

cos29 

f(ap) 

u 

f  (apcosco) 

—  (b.j.  )  a  '  Cr3 


where  9  is  the  orientation  angle  of  P'  with  respect  to  P  (measured 
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counterclockwise  from  east),  (u^a)a  and  (^t2)»  are  the  variances 
computed  from  the  averaged  geostrophic  winds,  and  in  which  we 
have  written  pcoscp  as  the  length  of  the  averaging  interval  for 
the  meridional  component  which  effectively  designates  p  as 
measured  in  degrees  (in  the  one  case  degrees  of  latitude  and  in 
the  other  case  degrees  of  longitude). 

It  is  readily  seen  that  the  correction  factor  as  a  function 
of  0,  oscillates  between  l/f^p)  and  l/f  (apcoscp) .  The  correction 
applied  to  the  correlations  appears  in  the  forms 

ri,  =  ru  -  (rtt).v6^T 

rAt  =  ( c f  t,  -/,yi  W  rtt  =  Vtt).\/GaGa' 
in  which  (ra£')»  indicates  the  value  computed  from  average  com¬ 
ponents  and  the  prime  on  or  03  refers  to  the  point  P'. 

The  implications  of  the  preceding  analysis  to  the  University 
of  Wisconsin  correlation  data  lies  principally  in  the  interpreta¬ 
tion  of  the  value  of  k.  This  parameter  was  previously  described 
as  containing  the  effects  of  smaller  scale  eddies  and  observation 
error.  In  connection  with  the  correlations  derived  from  the 
University  of  Wisconsin  data,  it  must  be  Interpreted  in  terms  of 
the  differences  between  averaged  and  point  geostrophic  wind  com¬ 
ponents,  and,  of  course,  the  fact  that  averaged  geostrophic  com¬ 
ponents  instead  of  real  wind  components  are  used.* 


*0f  course,  winds  are  observed  by  techniques  that  make  use  of 
an  average  of  some  kind,  altitude,  time,  etc.  The  point  becomes 
important  for  practical  purposes  when  the  averaging  interval 
reaches  a  size  comparable  with  the  definition  required  for  the 
study  at  hand. 
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The  values  for  L~a  derived  from  the  location  of  the 
zeros,  minima,  etc.  of  the  correlation  functions  are  not  in¬ 
volved  in  this  correction  for  using  averaged  winds  since  the 
numerical  values  of  the  correlations  are  not  involved  direction. 
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REPRESENTATION  OP  THE  CORRELATION  COEFFICIENT 


A.  CONDITIONS  IMPOSED  FROM  SYMMETRY 

Let  H  be  an  observed  quantity,  a  function  of  a  single  vari¬ 
able,  H(x),  and  let  the  values  considered  by  the  departure  from 
a  suitable  mean  value.  Let  r(h,x)  be  the  correlation  coefficient 
of  H  with  itself  at  two  points,  x  and  x  +  h.  Assume  that  r(h,x) 
has  an  infinite  series  expansion  in  powers  of  h  with  coefficients 
which  are  functions  of  x.  Then 
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On  collecting  terms  in  like  powers  of  h  from  this  second  expres¬ 
sion  and  equating  these  to  the  coefficients  of  the  first  expan¬ 
sion  in  powers  of  h,  one  obtains  the  system  of  equations: 


It  is  seen  at  once  that  the  even  subscripted  terms  are  independent 
but  that  the  odd  subscripted  terms  depend  on  the  preceding  even 
subscripted  terms.  Thus 


ax  -  £*a'. 


Since  as  a  rule  one  is  interested  in  the  case  of  a0  -  1,  then 
r(h,x)  =  1  +  |^-h2  +  £*§^-h3  +  +  ( £  *  ^T‘| 

+  +  (i-sf  -  Trft +  ^tr)h7  +  — • 

It  is  seen  that,  if  the  "shape"  of  the  correlation  function 
changes  from  place  to  place  (i.e.,  the  coefficients  are  indeed 
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non-trivial  functions  of  x,  a„  ^  constant,  n  =  2,  3,  - )  then 

the  correlation  function  expansion  contains  odd  powers  of  h  and 
the  correlation  curve  is  not  symmetrical  with  respect  to  h  =  0. 

The  situation  is  somewhat  more  general  than  indicated  in 
the  above.  If  it  is  assumed  that 


r(h,x)  =  1  -  A(x)|h|a  -  |h|am(h,x),  0  <  a  s  2 

where 

co(h,x)  -  b1(x)h  +  b2(x)|j  +  --- 
then  the  condition 


requires  that 


r(h,x)  =  r(-h,x+h) 


A  =  A, 


bj  =  -  b4  +  A', 


bo  bo 

Of  ”  Of 


b' 


+ 


A" 

2p 


£a.  = 

3T 


3T 


+ 


<  .  £1 
2 T  2T 


+ 


A'" 

3 V 


These  are  the  same  equations  as  obtained  before  with  n0  =  1,  ax  =  0, 
a2/2 1  =  -  A,  a3/3i  =  -  bx ,  a*/4J  =  -  b3,  etc. 

Consequently,  the  correlation  function  need  not  approach  ae 

(or  l)  as  a  quadratic  (l  -  hs/2Aa  +  - )  but  may  approach  its 

limit  value  in  any  of  the  permissible  (0  <  a  s  2)  ways.  When  the 
field  of  variable  correlated  is  not  homogeneous  (i.e.,  when  the 
coefficients  of  the  expansion  are  not  constant,  independent  of  x) 
then  the  correlation  function  will  contain  "odd"  power  terms  and 
will  not  be  symmetrical  with  respect  to  the  value  h  =  0. 
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A  somewhat  different  aspect  of  this  property  is  discussed 
in  the  next  section  in  the  case  of  a  differentiable  field  of 
property  being  correlated. 


A-III-5 


B.  EXPANSION  FOR  AN  ANALYTIC  FUNCTION 

If  H(x)  possess  a  power  series  expansion  (is  an  analytic 
function),  then 

H(x+h)  =  H  +  H'h  +  ^  +  — 

c  •  n  i 

and  the  covariance  of  H  is  given  by 

H(x)H(x+h)  =  (H2)  +  (ffiP^h  +  +  —  +  (HH^  ^  )~y  +  — 

where  it  is  assumed  that  the  means  appearing  as  coefficients 
exist.  In  terms  of  the  correlation  coefficient 

H(x)H(x+hj  =  o(x)a(x+n)r (h,x) 

=  c0(ao  +  aob  +  0ojrr  +  — )(l  +  aih  +  as|r  +  — )• 
Multiplying  the  two  series  and  equating  the  coefficients  with 
those  of  the  other  expansion  above, 

ii*  =  a02, 

(HH' )  =  a0o0a  +  a0a', 

(HH7'  )/2  !  =  a3  o02/2  I  +  a^a,'  +  a0cx"/2i, 

(HH'")/3  •'  =  a3a02/3i  +  a3c0c0'/2i  +  a  ^3qI/2\  +  a0a;"/3J, 

From  the  first  and  second  equations,  it  is  seen  at  once  that 
aa  =  0  (as  before).  The  left-hand  terms  may  be  developed  as 
follows : 

(iF)  =  o02 

HiT  =  a0a'  =  £(c52)' 


(HH")  +  (H'H')  =  |(  c  0  2 ) " ,  (HH")  =  -  ffx2  +  £(o02) 
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where  )  ==  a^s .  Continuing, 

(HIT)  +  (H'H"  )  =  -  K2)'  +  i(cr0a)/'/,  (ipH""] )  = 
so  that 

(HHW)  =  -  |-(a12)'  +  %(o0s)m. 

Further  _ 

(HH*’)  -  °„s  -  2(0'  +  mo1' 

(hF)  =  |(0'  -  |(0" 

where  (H"H" )  =  a32,  etc. 

Substituting  the  values  for  the  covariances  (HH^n  ^ )  and 
solving  for  the  an  ‘a,  one  obtains 


The  expressions  for  the  coefficients  become  quite  complex  when 
the  conditions  for  homogeneity  are  not  satisfied,  i.e.,  o'a  /  0, 
etc. 
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deviation  of  wind  and  the  standard  deviation  of  height  where 
L  plays  the  role  of  the  parameter  that  describes  the  size  of 
height  disturbances. 

If  II  represents  a  wind  component,  say  the  east  pointing 
component,  u,  and  x  Is  In  the  vertical  direction,  then  is 
the  standard  deviation  of  the  east  pointing  component  of  the 
wind  shear.  The  value  of  L2  may  be  obtained  from  the  radius  of 
curvature  of  the  correlation  of  wind  with  height  at  h  =  0  (this 
corresponds  to  a  "depth"  of  wind  systems).  It  Is  seen  that  the 
standard  deviation  of  wind  shear  is  made  up  of  two  parts,  one  is 
the  combination  of  the  standard  deviation  of  wind  divided  by  a 
depth  parameter,  the  other  is  due  to  the  non -homogeneity  of  wind 
with  height.  Note  that  in  the  above  application  one  is  finding 
the  standard  deviation  of  -wind  shear  due  to  the  differentiable 


part  of  the  wind.  It  Is  necessary  to  consider  the  relation 


r*  =  k0(l  -  h2/2L2  +  —) 
in  finding  L  and  to  note  that 

k0  =  1/(1  +  e2/ du  2 ) 
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where  es  is  the  mean  square  value  of  the  observation  error  and 
the  non-differentiable  part  of  the  wind  (i.e.,  the  part  due  to 
small  eddies  of  depth  less  than  a  depth  commensurate  with  the 
averaging  depth  of  wind  observations). 


APPENDIX  A -IV 


GAUSS'  PROCEDURE 

In  the  case  of  the  parameters  A,  B,  C,  D,  E,  and  F  of  the 
equations  on  page  V-27,  a  simple  least  squares  procedure  is 
impossible.  As  a  consequence,  an  iterative  method  of  solution 
is  used.  The  major  steps  of  the  procedure  are  outlined  here. 

The  equations  consist  of  pairs  which  may  be  represented 
schematically  as  follov/s: 

Yi  —  P*!  (r ,  9  J  A^ ,  —  A0  ) , 

Ya  =  Fa  (r,e;Ax,—  ,A6  ) . 

The  estimates  on  the  right  are  considered  equally  important  in 

the  determination  of  the  parameters  A1 ,  - ,  Aa  (corresponding 

to  A,  — ,  F  respectively).  Consequently,  the  sum  of  the  squares 
of  the  residuals  from  the  two  equations  together  are  minimized. 
Thus 

£[(yi"Fi)2  +  (Ya-P’a)2]  =  Minimum. 

The  usual  procedure  of  differentiating  with  respect  to  the 
parameters  and  setting  these  derivatives  to  zero  leads  to  the 
equations 

Z  “Pl  )  aft"  +  "Fa  =  °’  1  =  lj  * 

This  produces  six  equations  in  the  six  parameters,  Ax ,  - ,  A0 . 

The  parameter  values  are  involved  in  complex  expressions  which 
cannot  be  solved  simply;  i.e.,  each  of  the  six  equations  is  so 
complex  that  it  cannot  be  solved  for  the  parameter  values  in 
the  usual  explicit  form. 
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Consequently,  an  iterative  method  of  solution  is  used . 


This  is  obtained  by  using^  the  Taylor  expansion  of  the  functions 
Fj  and  Fa  . 


Pi  « 


-j-  ^Pl  Al\  -I- _ ^Fj  A  A  +  _  _  _ 

+  3A^Al  +  + 


P  _  p  (°  )  O.  ^Pg  Aft  J. 4.  SFb  a  ft  I 

*a  “  +  SA/^3  Me^A®  * 

where  Fx  and  Fa  ^  Indicate  that  the  function  is  evaluated 

at  the  initial  parameter  values  A1^°'IJ - ,  Ae^.  It  is  assumed 

that  the  derivatives  are  evaluated  at  the  same  point,  but  this 
is  not  shown  specifically  to  keep  notation  complications  within 
reason.  When  this  substitution  is  made  and  the  terms  rearranged, 
the  preceding  equations  become: 


The  bars  inside  the  square  brackets  are  to  call  attention  to 
the  fact  that  these  expressions  arc  summed  over  the  data  values. 

The  higher  order  terms  on  the  left,  indicated  by  - ,  are  neglected. 

All  of  the  coefficients  are  evaluated  at  the  initial  parameter 
estimate,  Al^0^ . - ,  Aa^°^.  The  set  of  equations  is  then  solved 
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for  AA^0), - ,  AAa^0^.  These  determine  the  new  estimates, 

Ax  ^  =-•  Ax  (°  ^  +  Mx  (°  ^  — ,  A6  ^  )  =  Aa  (°  ^  +  AAa  (°  ^ .  The  co¬ 
efficient  values  are  then  recomputed  from  the  new  parameter 
values,  Aj^1).  ,  A0^  and  further  increments  are  obtained 


+  AA. 


(O  ...  A  (■}  = 


The  process  is  continued  until  the  additional  increments  are 
considered  sufficiently  small  or  until  some  other  convergence 
criterion  is  satisfied. 

In  the  particular  case  at  hand,  the  iterative  sequence 
was  stopped  when  AA/A  s  10"8 .  This  insured  a  reasonably  small 
amount  of  computation  and  in  view  of  the  small  change  in  the 
residuals  this  degree  of  accuracy  was  considered  sufficient. 

The  residuals  were  calculated  directly  from  the  initial 
expression  to  be  minimised.  In  using  such  an  Iterative  pro¬ 
cedure,  the  indirect  estimate  from  the  calculated  coefficients 
and  the  non -homogeneous  terras  of  the  normal  equations  do  not 
provide  a  legitimate  estimate  of  the  sum  of  squares  of  the 


residuals . 
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